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Abstract 

This is a comprehensive review on fault-tolerant topological quantum computation 
with the surface codes. The basic concepts and useful tools underlying fault-tolerant 
quantum computation, such as universal quantum computation, stabilizer formalism, 
and measurement-based quantum computation, are also provided in a pedagogical 
way. Topological quantum computation by brading the defects on the surface code 
is explained in both circuit-based and measurement-based models in such a way that 
their relation is clear. The interdisciplinary connections between quantum error cor¬ 
rection codes and subjects in other fields such as topological order in condensed matter 
physics and spin glass models in statistical physics are also discussed. This manuscript 
will be appeared in SpringerBriefs. 
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Preface 


In 1982, Richard Feynman pointed out that a simulation of quantum systems on classi¬ 
cal computers is generally inefficient because the dimension of the state space increases 
exponentially with the number of particles |Fey82| . Instead, quantum systems could be 
simulated efficiently by other quantum systems. David Deutsch put this idea forward by 
formulating a quantum version of a Turing machine |Deu85| . Quantum computation en¬ 
ables us to solve certain kinds of problems that are thought to be intractable with classical 
computers such as the prime factoring problem and an approximation of the Jones poly¬ 
nomial. It has a great possibility to disprove the extended (strong) Church-Turing thesis, 
i.e., that any computational process on realistic devices can be simulated efficiently on a 
probabilistic Turing machine. 

However, for this statement to make sense, we need to determine whether or not 
quantum computation is a realistic model of computation. Rolf Landauer criticized it 
(encouragingly) by suggesting to put a footnote: “T/izs proposal, like all proposals for 
quantum computation, relies on speculative technology, does not in its current form take 
into account all possible sources of noise, unreliability and manufacturing error, and prob¬ 
ably will not work. ” |Llo99j . Actually, quantum coherence, which is essential for quantum 
computation is quite fragile against noise. If we cannot handle the effect of noise, quantum 
computation is of a limiting interest, like classical analog computers, as a realistic model 
of computation. To solve this, many researchers have investigated the fault-tolerance of 
quantum computation with developing quantum error correction techniques. One of the 
greatest achievements of this approach is topological fault-tolerant quantum computation 
using the surface code proposed by R. Raussendorf et al. |RHG06l IRHG071 IRH07j . It says 
that nearest-neighbor two-qubit gates and single-qubit operations on a two-dimensional 
array of qubits can perform universal quantum computation fault-tolerantly as long as the 
error rate per operation is less than ~ 1%. 

In this book, I present a self-consistent review of topological fault-tolerant quantum 
computation using the surface code. The book covers everything required to understand 
topological fault-tolerant quantum computation, ranging from the definition of the sur¬ 
face code to topological quantum error correction and topological operations on the surface 
code. The basic concepts and powerful tools for understanding topological fault-tolerant 
quantum computation, such as universal quantum computation, quantum algorithms, sta¬ 
bilizer formalism, and measurement-based quantum computation, are also introduced in 
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the first part (Chapter 1 and Chapter 2) of the book. In particular, in Chapter 1, I also 
mention a quantum algorithm for approximating the Jones polynomials, which is also re¬ 
lated to topological quantum computation with braiding non-Abelian anyons. In Chapter. 
3, the definition of the surface code and topological quantum error correction on it is ex¬ 
plained. In Chapter 4, topological quantum computation on the surface code is described 
in the circuit-based model, where topological diagrams are introduced to understand the 
logical operations on the surface code diagrammatically. In Chapter. 5, I explain the 
same thing in the measurement-based model, as done in the original proposal [RHCObj . 
Hopefully, it would be easy to see how these two viewpoints are related. 

Throughout the book, I have tried to explain the quantum operations using circuit 
and topological diagrams so that the readers can get a graphical understanding of the 
operations. The graphical understanding should be helpful to study the subjects more 
efficiently. Topological quantum error correction codes are a nice play ground for studying 
the interdisciplinary connections between quantum information and other fields of physics, 
such as condensed matter physics and statistical physics. Actually, there is a nice corre¬ 
spondence between topological quantum error correction codes and topologically ordered 
systems in condensed matter physics. Furthermore, if we consider a decoding problem of 
a quantum error correction code, a partition function of a random statistical mechanical 
model is naturally appeared as a posterior probability for the decoding. These interdisci¬ 
plinary topics are also included in Chapter 3. 

Almost all topics, except for the basic concepts in the first part, are based on the results 
achieved after the appearance of the standard textbook of quantum information science 
entitled “Quantum Computation and Quantum Information” (Cambridge University Press 
2000) by M. A. Nielsen and I. L. Chuang. In this sense, the present comprehensive review 
on these topics would be helpful to learn and update the recent progress efficiently. In 
this book, I concentrated on the quantum information aspect of topological quantum 
computation. Unfortunately, I cannot cover the more physical and condensed matter 
aspects of topological quantum computation, such as non-Abelian anyons and topological 
quantum field theory. In this sense, this book is complemented by the book “Introduction 
to topological quantum computation” (Cambridge University Press 2012) written by J. K. 
Pachos. Readers who are interested in the more physical aspects of topological quantum 
computation are recommended to read it. 

Hopefully, this review will encourage both theoretical and experimental researchers to 
find a more feasible way of quantum computation. It will also bring me great pleasure 
if this review provides an opportunity to reunify and refine various subdivided fields of 
modern physics in terms of quantum information. 


Kyoto, Japan 


Keisuke Fujii, 
April, 2015 
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Chapter 1 

Introduction to quantum 
computation 


In this chapter, we introduce the basic concepts of quantum computation. We first describe 
the minimum unit of quantum information, the qubit, and define several gate operations for 
it. Then, we explain the Solovay-Kitaev algorithm, which provides a way to decompose an 
arbitrary single-qubit unitary operation into an elementary set of single-qubit gates. Using 
multi-qubit gates, we construct an arbitrary n-qubit unitary operation from an elementary 
universal set of gates, which we call universal quantum computation. Quantum algorithms, 
which run on a universal quantum computer, are also presented. One example of this is 
Shor’s prime factoring algorithm based on the phase estimation algorithm. Another is an 
approximation of the Jones polynomial. Finally, we will introduce quantum noise and see 
how we can describe a quantum system coupled with an environment. 

1.1 Quantum bit and elementary operations 

In classical information science, the minimum unit of information is described by a binary 
digit or bit, which takes the value 0 or 1. Its quantum counterpart is a quantum bit, the 
so-called qubit. The qubit is defined as a linear superposition of two orthogonal quantum 

states 1^) ~ ^ Q ^ ~ 1 

Itjj) = a\0) + (3\1), (1.1) 

where a and j3 are arbitrary complex values satisfying |ap-|-|/3p = 1. The complex 
amplitudes can be expressed as 

6 9 

q; = cos-, /3 = e*‘^sin-, (1.2) 

up to an unimportant global phase. By using the angles 9 and (j), the qubit can be mapped 
onto a point on a three-dimensional (3D) sphere, the so-called Bloch sphere, as shown in 
Fig. [Lll 
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Figure 1.1: The Bloch sphere. 


The time evolution, which maps a quantum state into another, is given as a unitary 
operator in quantum mechanics. The most important operators are the Pauli operators 


I = 


1 0 
0 1 


X = 


0 1 
1 0 


Y = 


0 -i 
i 0 


Z = 


1 0 
0 -1 


(1.3) 


The computational basis states {|0),|1)} are eigenstates of the Pauli Z operator. The 
Pauli X operator flips the computational basis state: 


|1)=X|0), |0)=X|1). 

We define the eigenstates of the Pauli X operator as 


l+) = 


| 0 ) + | 1 ) 
V2 




| 0 ) - | 1 ) 

^/2 ^ 


(1.4) 


(1.5) 


which we call the X-basis states. Similarly, the eigenstates of the Pauli Y operator are 
defined as 


+ i) = 


| 0 )+^| 1 ) 

V2 


- ^ = 


V2 ^ 


( 1 . 6 ) 


which we call the T-basis state. 

The second most important operators are the Hadamard H and phase S operators 


1 


1 0 


H = — (] \ ) and 5 = . ^ . 

y/2\l -I J V 0 z 
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These gates transform between the different Pauli-basis states, i.e., H : {|0),|1)} -H- 
{|+),|—)} and S : {|+),|—)} -H- {| +i))l “ Ol- Equivalently, we may say that these 
operators transform a Pauli operator into another Pauli operator under their conjuga¬ 
tions: 


X = HZH, Y = SXSl (1.7) 

From this property, the H and S gates are called Clifford gates. The Clifford gates are 
depicted as circuit diagrams as follows; 



If we measure the qubit \ip) in the Z-basis {|0),|1)}, we obtain the measurement 
outcomes 0 and 1 with probabilities 


po = |(0|V')l2 = Tr[|0)(0|V’)(V’|], (1.8) 

PI = |(l|V^)p = Tr[|l)(l|V^)(V^|], (1.9) 


respectively. More generally, we may use the measurement operators {Mj}, satisfying 
that Ei = M^Mi is positive semidefinite and that Yli ^■ (If operator A satisfies 

^|u), {u\A\u) > 0, it is said to be positive semidefinite.) The probability of obtaining the 
measurement outcome i is given by 


Pi = (V'IM/M.IV') = TV 




TT[Ei\i;)m. 


( 1 . 10 ) 


Such a measurement and the set of positive operators {Ei = are called a positive- 

operator-valued measure (POVM) measurement and POVM elements, respectively. The 
post-measurement state conditioned on the measurement outcome i is 


V^Tr[M/M,lV')(V'|] 


( 1 . 11 ) 


Suppose we have quantum states |'0) and |(/>) with probability p^ and respectively. 
We perform a measurement with the measurement operator {Mj}. If we assume that the 
measurement outcome i is obtained with probability 


Pi=p^Yi Tp^Tr 


M^Mi 


we can express the classical mixture of and \4>) as an operator 


( 1 . 12 ) 




(1.13) 





which is called the density matrix or density operator. More generally, if we have pure quan¬ 
tum states (IV'fc)} with probability {pk}-, the density matrix is given by p = YlkPk\'4’){i’]- 
A density matrix, which has information of a quantum state including its statistical prop¬ 
erty, has to be a positive hermitian (self-adjoint) operator p satisfying Tr [p] = 1. Note that 
for a given density matrix p, the pure state decomposition of it is not uniquely determined, 
and any decomposition has the same statistical property. 

A mixed state of a qubit p can be represented as a point inside the Bloch sphere 
through the three coordinates calculated from the density matrix 

(r-.,r„r,) = (Tr[Xp], Tr[yp], Tr[Zp]). (1.14) 

If the state is a pure state |V’) = cos f |0) -|- sin ||1), the coordinates can be calculated 
to be 


{rx,ry,rz) = (sin0cos(/), sin0sin()), cos0), (1-15) 

which is consistent with the previous definition. The completely mixed state 1/2 corre¬ 
sponds to the origin of the coordinate system. 

1.2 The Solovay-Kitaev algorithm 

The Pauli operators {I, X,Y, Z} and single-qubit Clifford operators {H,S} form a finite 
group, and hence cannot cover all unitary operations for a qubit. If a non-Clifford oper¬ 
ation exists, e.g., ^ we can generate an arbitrary single-qubit unitary operation 

using the Solovay-Kitaev algorithm |DN05] . The underlying Solovay-Kitaev theorem states 
that if a set of single-qubit operations generates a dense subset of 517(2), then that set is 
guaranteed to fill 51/(2) quickly. 

Suppose we have a basic (0th order) approximation Uq of an arbitrary unitary operator 
U, and that it approximates U with a certain constant error cq: 

\\Uo-U\\<eo, (1.16) 

where ||...|| indicates an operator norm. The Solovay-Kitaev algorithm takes the (n — l)th 
order approximation Un-i with an error e^-i and returns the nth order approximation 
Un with an error as follows. 

First, UU\_i is decomposed in terms of the unitary operators V and IF as a group 
commutator: 


C7C7t_^ ^ FIFFtvi/t^ 

where V and IF are chosen such that 

\\V-I\\<c^f^i, ||IF-/|| <Cv/i;(ZT, 


(1.17) 


(1.18) 
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with a constant c. The above decomposition always exists from the following argument. 
Let V and W be rotations with an angle (p about the x- and y-axes, respectively. Then 
VWV^W^ is a rotation with an angle 9 about some axis, where 


sin(0/2) = 2siv?{(p/2)\J 1 — sin^((/>/2). (1-19) 

We assume that 0, 0 <C 1, and hence 9 ~ cjP'. By inverting the above argument, if 

\\I-VWV^W^\\c^9/2 + 0{9^) ( 1 . 20 ) 

we can find V and W such that 

||/- y ||, ||/- W || ~(/./2 + 0 (</> 3 ). ( 1 . 21 ) 

Because 

||/ - uul_^\\ = \\I - < e„_i, (1.22) 

we can find V and W such that ||/ — y||, ||/ — 1T|| < c^e„_i. 

Second, we calculate the (n — l)th order approximations W-i and ITn-i of V and W, 
respectively. Then the Solovay-Kitaev algorithm returns the nth order approximation 

Un = Vn-lWn-lV^_,Wl_^Un-l, (1.23) 

which satisfies 

\\Un-U\\<en. (1.24) 

Next, we will calculate as a function of e^-i. By using the property of the operator 

norm, we obtain 

\\Un-U\\ < \\Vn-lWn-lV^_,Wl_^ - UUl_MUn-l\\ (1.25) 

< \\Vn-lWn-lV^_,Wl_,-VWV^W^\\\\Un-l\\ (1.26) 

By denoting Vn-i = V + Ay, ITn-i = W + A^y, we obtain 

\\Un-U\\ < \\Vn-lWn-lV}^_^Wl_^-VWV^W^\\ (1.27) 

< WAvWV^W^ + VWAIw^W + WVAwV^W^ + VWV^Aly\\ + 0{A^). 

(1.28) 

Moreover, in terms of Sy and Sw defined hy V = I + Sy and W = W + 5w respectively, 
we obtain 

\\Un-U\\ < llAyFt+ yA[.|| + ||AiyWt +WA[^|| (1.29) 

+ ||Ay(5u/|| + ||A|^(5p^^|| + ||(5yAy/|| + ||(5yA|^|| + O(A^). (1.30) 

< 4ce^/_\+ 0(A2,,52A). (1.31) 


10 



Here we have used that 


AyW + (1.32) 

which can be derived from the unitarity of V and Vn-i = V + Ay. To leading order, the 
error is given by en < with a certain constant d and is calculated to be 

en < (2c'eo)(^/2)7(2c')- (1-33) 

If eo < l/(2c'), the error decreases super-exponentially in n. On the other hand, the nth 
order approximation calls the (n — l)th order approximation three times, i.e., through 
Un-i, 14 .- 1 , and H4_i. Including the resource Rd required for the decomposition (jl.I7l) . 
the overhead Rn for the nth order approximation is given by 

Rn = 3Rn-i + Rd (1-34) 

^ Rn = 0{3^). (1.35) 

Similarly, the number of unitary operations employed in the nth order approximation is 
calculated to be M„ = 0(5"). By using (11.331] and (ll.35|] . the overhead R^ and the number 
Me of gates required to obtain an approximation of U with an error e can be estimated: 

Re = 0(ln'°3/M3/2)(i/g))^ 

Me = 0(ln'°®/^"(3/2)^^/g)^_ (^37^ 

Thus, both Re and Me scale as polylogarithmic functions of 1/e. 

1.3 Multi-qubit gates 

An n-qubit state is given by a superposition of tensor product states 

I'h) = 'y ^ Oqj2...jn 

11,12,...,in 

where = 0,1 and |iii2---4) = |h) <8) 14) <8> ••• <8> |4)- A single qubit gate A acting on the 
/cth qubit is denoted by 

k—1 n—k—1 

' -""-^ ^ 

= (1.38) 

An important two-qubit gate is the controlled-NOT (CNOT) gate, 

Ae,t{X) = |0)(0|e4 + |l)(llcXi. 

For a computational basis input state \i)c\j)t, the CNOT gate acts as kc,t{X)\i)c\j)t = 
\i 0 j). In this sense, the CNOT gate is a quantum generalization of the XOR operation 
in classical computation. If the input state is |+)c|0)t, the output of the CNOT gate is a 
maximally entangled state: 

A,,t(A)|+),|0)i = (|00) + |11))/72. (1.39) 

The CNOT gate is depicted by a circuit diagram as follows: 
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We may also define the controlled-Z (CZ) gate, 

Ae,t(Z) = \0){0\Jt + \l){l\cZt. 

The CZ operation is symmetric; Ac^t{Z) = At^c{Z), and is represented by the circuit 
diagram 



Because the Pauli Z operator is transformed into the Pauli X operator by the Hadamard 
gate, we have the following relation between the CZ and CNOT gates; 



These CNOT and CZ gates are both Clifford gates, i.e., A(yl) {A = X, Z) transforms 
the two-qubit Pauli group onto itself under the conjugation A(yl)[- • •]A(yl)^. For example, 
A{X),^tX,®ItA{X)l^ = X,0Xt, A{X),^tIc^ZtA{X)l^ = Z,0Zt, A{Z),^tXc0ltA{Z)lt = 
Xc® Zt, etc. 

For an arbitrary unitary operator U, the controlled-C/ gate is denoted by 

Ae,i(C/) = |0)(0|,/t + |l)(l|eC/t, (1.40) 

where the qubits c and t are called the control and target qubits, respectively. By de¬ 
composing a single-qubit unitary gate into U = e^°‘AXBXC with the unitary operators 
A,B,C satisfying ABC = I, the controlled-?/ operation A(?7) can be implemented as 
follows; 



12 





(Note that we always have such a decomposition for an arbitrary single-qubit gate |BBC'*~95 
iNdnnj .i 

For example, the controlled-Hadamard gate A{H) can be represented by 



Next, we will discuss one of the most important multi-qubit gates, the Toffoli gate: 

= ilcjc, - |l)(l|ejl)(l|c,)/i + |l)(l|ci|l)(l|c2^i, (1.41) 

which is depicted as a circuit diagram: 



For a computational basis input state |ii)cil* 2 )c 2 |j)t) the Toffoli gate acts as 

^Lc2,i(^)l*i)ciN2)c2lj)i = \k)cAi2)c2\j ® (n • k))t- (1.42) 

The state of the third qubit is equivalent to the output of the NAND operation in classical 
computation. In this sense, the Toffoli operation can be regarded as a quantum extension 
of the NAND operation. The NAND operations are known to be universal in classical 
computation in the sense that any logic gate (boolean function) can be constructed from 
them. This implies that quantum computation trivially includes classical computation. 
More importantly, because unitary operations are reversible, quantum computation can 
simulate classical computation in a reversible way. Suppose we want to calculate a boolean 
function f{x) for an input state x. Then, we can construct a quantum circuit Uf consisting 
of Toffoli and Pauli X gates: 


kf\x) input 10.. .0)ancilla 10) answer — |2:)input |fl'(3:)) anciUa|/(3:')) answer; (1.43) 


where the qubits |•)input, |•)ancilla) and |•)answer indicate the registers for the input state, 
the ancillae for the Toffoli operations, and the answer of the calculation, respectively. The 
output state of the ancilla register |g(x))anciiia is the garbage of the computation. However, 
the garbage can be uncomputed as follows (see also the circuit diagram below): 


(/JA 

answer,out {X)Uf\ x) input 1 0...0)ancilla 1 0) ar 
I input Iff (l^)) ancilla I f (®)) answer I f (®))out 
— |3^)input 1 0...0) ancilla 1 0) answer I / (lZ'))out • 


•| 0 ) 


out 


(1.44) 
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\x) 

I00...0) 

| 0 ) 

l/(^)) 


In this way, we can calculate an arbitrary boolean function in a reversible way. 

Finally, we introduce the multi-controlled unitary gate 

A^(C/) = - |1)(1|®^) (8) / + |1)(1|®^ (8) U, (1.45) 


where U is applied to the target qubit if all k control qubits are |1). The multi-controlled 
gate can be implemented by using the Toffoli gates and k ancilla qubits as follow; 



2< 


|0). 
I |0)- 
|| 0 ). 

o 
c 
CC 

| 0 )- 

target qubit 


■ 



s 

■ 


■ 


■ 

1 

1 



1 

1 

B 


■ 

■ 


^^R 

■ 

■ 



■ 

■ 


^^R 

■ 

m 



m 



m 










The Is?{U) gate can be decomposed into CNOT and single-qubit gates by using an idea 
similar to the decomposition of the A(t/) gate: 



For example, the Toffoli gate can be constructed from the CNOT, Hadamard, and vr/S 
operations as follows: 
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—( 

)— 



J' — 


where we have used the fact that that the controlled-\/X gate is decomposed into the 
CNOT, Hadamard, and tt/S (T = gates as follows: 



It should be noted that a multi-controlled gate A^{U) can also be constructed without any 
ancilla qubit from multi-controlled gates with lower controlled qubits by using 

the trick employed for decomposing h?{U) into A(C/) [BBC'*~95 . 


1.4 Universal quantum computation 


Here we will show that an arbitrary unitary operation can be decomposed into single¬ 
qubit and CNOT gates. To this end, we will first show how to decompose an arbitrary 
unitary operation into a product of two-level unitary gates. Second, the two-level gates are 
decomposed into the multi-controlled gates, which can be constructed from single-qubit 
and CNOT gates as seen in the previous section. 

Let U be an arbitrary n-qubit unitary operator, represented by an mxm unitary matrix 
with m = 2”. Let Tij be a unitary operator such that = 6ki k,l ^ i,j, which we 

call a two-level unitary gate. (The kl element of a matrix A is denoted by {A)ki.) The ii, 
ij, ji, and jj elements define a unitary operation on the two-level subsystem spanned by 
|j) and |i). By choosing Tmm-i appropriately, we have 


UTmm-l — 


Uu 

'^m—l 1 
\ '^m 1 


'^1 m—1 


Ui, 


1 m—1 ^m—1 m 

0 11^ 

^ m 


where {U)ki = Uki- By repeating this procedure, we obtain 


/ 


UTjyi m—lTm m—2 ' ' ' i — 


ni 


^m—1 1 


m—1 ^1 m '' 


V 0 


m—1 m—1 m—1 m 

0 u" 


(1.46) 


(1.47) 
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Due to unitarity, u'{^ = ■■■ = u'^_^ m = 0 and = 1. Defining m-iTm m -2 ■■■Tmi, 

we can decompose U into a product of Rk and a diagonal unitary operator D: 

U = D{Rm---Ri)^. (1.48) 

It is obvious that D can be decomposed into two-level unitary gates. Thus, an arbitrary 
unitary operator U can be decomposed into two-level unitary gates. 

Next, we show that any two-level unitary operator Rj can be implemented by using 
CNOT and single-qubit gates. Let us rewritte i and j {i,j = 0,..., m — 1) by using the n-bit 
strings s = siS 2 ...Sn and t = respectively. It is easy to find a sequence of n-bit 

strings such that s = gi, t = and g^ and gfc+i differ by only one bit. By using 

the Pauli X gate and the multi-controlled-NOT gate controlled by the same n — 1 bits 
and targeting the one different bit, we can transform the basis |gfc) to |gfc_|_i). In this way, 
the basis \i) = |s) = |gi) is transformed into Ig^-i). After this basis transformation, we 
now want to apply a two-level unitary gate between |gd_i) and \j) = |t) = Ig^^). Because 
gd_i and gd differ by only one bit, we can perform such a two-level unitary gate by the 
multi-conditional gate, using the Pauli X gate for bit flips. Finally, the basis |gd-i) is 
returned to |s) by applying the inverse of the basis transformation. 

For example, a two-level unitary operator acting on a subspace spanned by {|000), |111)} 
can be implemented as follows; 



| 000 > | 011 ) | 000 > | 011 ) 


As seen previously, the multi-conditional gate can be decomposed into CNOT and 
single-qubit gates. Moreover, an arbitrary single-qubit unitary operation can be approx¬ 
imated by using the Hadamard and tt/S operations by virtue of the Solovay-Kitaev algo¬ 
rithm. Thus, the CNOT, Hadamard, and vr/S operations form a universal set of operations 
for quantum computations. The Toffoli and Hadamard operations also form a universal 
set as shown in [RZBB941 lBBC~*~95| IDiV95( IBV931INCOOI . 


1.5 Quantum algorithms 

In this section, we explain two representative quantum algorithms; Shor’s prime factoriza¬ 
tion algorithm [Sho94[ ISho97] and the Aharonov-Jones-Landau algorithm for an additive 
approximation of the Jones polynomial |AJL09[|AA06] . We do not go deep into the math¬ 
ematically rigorous details, but we aim to understand how they work. The readers who are 
interested in more details should read Refs. |Sho94( ISho97( IKit95at INCDOt lAJLDht lAAObj . 
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1.5.1 Indirect measurement and the Hadamard test 


An indirect measurement of an observable (hermitian) A with eigenvalues ±1 can be 
performed by using A(yl): 



By denoting the input state IV'), the post-measurement state is given by 

^ + (-i)M)/ 2 |.»)w]. ( 1 . 49 ) 

where s = 0,1 is the measurement outcome. For example, a circuit measuring the eigen¬ 
value of the operator X 1 X 2 X 3 is given by 



According to the measurement outcome s = 0,1, the post-measurement state is projected 
by ^+(~^) XiX 2 Xi ^ indirect measurement will be employed frequently in quantum 

error correction to measure the eigenvalues of the stabilizer operators. 

The Hadamard test of an arbitrary unitary operator U is defined by the following 
circuit: 



The probabilities of the measurement outcomes 0,1 of the JA-basis measurement are cal- 


culated to be 

Po = ^ (1 + Re(0|®”t/|0)®”) , 

(1.50) 


Pi = ^ (1 - Re(0|®”t/|0)®”) . 

(1.51) 
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Similarly, the Hadamard test for the imaginary part is defined: 



Suppose we perform the Hadamard test N times and obtain the measurement outcome 0, 
No times. By using the Chernoff-Hoeffding bound, 


Prob 


No 

-N-P^ 


> e < 2e 


-2t^N 


(1.52) 


we can estimate the matrix element (0|®”'f7|0)®”' with an error e by repeating the Hadamard 
test N = poly(l/e) times. The Hadamard test is employed in various quantum algorithms 
such as approximations of the Jones and Tutte polynomials [AJL091 IAA061 IAAEL07] 
and the partition functions of statistical mechanical models [DICDV dNMD 1 ll IICBB141 
IMFT14] . 

Specifically, if we choose the input state to be a completely mixed n-qubit state. 



then, the Hadamard test provides the trace Tr[C/]/2"' of the unitary operator U. Such a 
restricted type of quantum computation is called a deterministic quantum computation 
with one clean qubit (DQCl) |KL98j . DQCl seems to be less powerful than universal 
quantum computation, because only one qubit is a pure state. However, DQCl can 
evaluate functions, which would be intractable on a classical computer, such as the Jones 
and Homefly polynomials [SJn8( IJWf)9j . spectral density function [KL98| . and fidelity 
decay |PBKLO04] . Recently, classical sampling of the output of DQCl with a few qubits 
measurements has been shown to be intractable unless the polynomial hierarchy collapses 
to the third level [MFF141IMNFT14] . 
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1.5.2 Phase estimation, quantum Fourier transformation, and factoriza¬ 
tion 

If we have an eigenstate \Ei) of a unitary operator [/, for which A(C/) is described by a 
polynomial number of gates, we can estimate the eigenvalue A* of U with a polynomial 
accuracy using the Hadamard test as follows: 



X 


When we do not have the eigenstate, the input state is projected onto one of the eigenstates 
by the Hadamard test. Thus, by repeating the Hadamard test, we can obtain one of the 
eigenvalues with polynomial accuracy. 

Moreover, if a controlled-!/^ gate A(?7^ ) can be described by a polynomial number 
of gates, we can estimate the eigenvalue with exponential accuracy [Kit95b| . Suppose the 
eigenvalue is given by A* = ^ where we employ a decimal of a binary 

number, 

n 

0.jlj2...jn = (1-53) 

k=l 

Then, the Kitaev’s phase estimation algorithm is given by 

1 +) 

1 +) 

1 +) 


E^) 


rtk 

where QFT stands for the quantum Fourier transform. The controlled gate A{U‘‘ ) kicks 
back the phase to the ancilla state |-|-). The phase information 

is transformed into a computational basis state by using the inverse quantum Fourier 
transformation: 

QFrt(|0) +e2-0diA...i.|i))(|0) = \jij2...jn). 

(1.54) 







QFTf 


n\ 

rt\ 

z ] 
z > 



n\ 

_ > 
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Here QFT is given by 


\h) 


bn) 



> + e2’r«-J"-nn|i) 

> +e2’^®-J"|l) 


Then, we can obtain the phase (j) = {2'7Ti)0.jij2---jn through the Z-basis measurements. 
Note that the accuracy of the estimate is improved exponentially compared to the Hadamard 
test (provided that we have a polynomial size description of A{U^ )). 

Let N and x be co-prime integers. The order-finding problem is about finding an order 
r such that x'’ = 1 mod N. This problem can be solved by using the phase estimation 
against a unitary operator 


because we have 


Ux = '^ \xy mod N){y\, 
y 


Ux\Us) 


— g^Triis/r)) 


Us), 


(1.55) 


(1.56) 


for a state 

r-l 

\us) = ^ y e-2^*("/'’)^|x^mod N). (1.57) 

^ k=0 

Note that we can prepare the initial state |ms) randomly by using the fact that |1) = 
Z]I=o !'“«)• Aftsr estimating the phase 27rz(s/r), the continued fraction provides r with 
high probability. The modular exponentiation mod N can be calculated by using the 
square k times. Thus we can implement controlled-, A{U^ ), by a polynomial number 
of gates. 

If we randomly choose s, we obtain an even order r with high probability. Thus, we 
have (x^/^ — l)(x'’/^ -|-1) = 0 mod N. Finally, the euclidean algorithm gives us the greatest 
common divisor of x^/^ — 1 and N or x^/^ -|- 1 and N, and it is a factor of N. This is the 
so-called Shor’s prime factoring algorithm. Kitaev, who formulated the phase estimation 
algorithm, generalized this idea for the more general Abelian stabilizer problem [Kit 95 a] . 
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Figure 1.2: (a) A link diagram, (b) A braid diagram, (c) A tangle diagram. 



1.5.3 A quantum algorithm to approximate Jones polynomial 

Next we will explain Aharonov-Jones-Landau algorithm for approximating the Jones poly¬ 
nomial |AJLn9[ IAAn6| . which is an algorithmic version of the original proposal based on 
topological quantum field theory |FKLWn3| IWit89] . We will show that the approximation 
of the Jones polynomial is BQP-complete; the approximation of the Jones polynomial can 
be done efficiently by using a universal quantum computer, and inversely it is as hard 
as any problems (BQP) solvable by a universal quantum computation. Below we will 
introduce the braid diagram and the braid group, which are related to link diagrams with 
an appropriate closure. A unitary representation of the braid group is constructed by 
embedding the braid group to the Temperley-Lieb (TL) algebra. Using the constructed 
representation, the Jones polynomial and quantum algorithm are connected. 

Let us first define the Jones polynomial. The Jones polynomial is an invariant of a link, 
that is, closed loops in a 3D space, which are tangled in general |Lic97] . It is convenient 
to describe the link on a 2D projected space as shown in Fig. 11.21 (a), which we call a link 
diagram. Specifically two links are equivalent if two link diagrams can be transformed into 
each others under the Reidemeister moves: 



which appropriately reflect the continuous deformations in the 3D space. The Jones poly¬ 
nomial is calculated from a directed link diagram L as follows: (i) Smooth each crossing 
in two ways { ')' f, ^ }. Let s be the resultant diagram consisting of closed loops with 
no crossing, which we call a state, (ii) For each state s, we assign a weight 

IU(s) = (1.58) 

where d = — (A^ -|- A~‘^) is a complex constant called a loop value, |s| is the number of the 
crossing, and s"'' and s~ are the number of smoothings by ^ f and ^ , respectively. By 
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Figure 1.3: (a) A multiplication for two braid diagrams bi and 62 - (b) The generator Uj 
of the braid group, (c) The generator Ei of the TL algebra. 

taking the summation over all states, the Kauffman bracket (L) of the link L is defined: 

{L) = Y,W{s). (1.59) 

S 

The Jones polynomial is defined as a function of t = by multiplying a factor to the 
Kauffman bracket: 

VL{t) = {-A)^^^\L), (1.60) 

where the writhe oj[L) of the link L is defined as the number of typs crossings minus 
that of JX typs- If is easy and a good exercise to confirm that the Jones polynomial is 
invariant under the Reidemeister moves based on the above definition [Lic97] . 

Let Bn be a braid group consisting of the braid diagrams of n strands, where two 
endpoints of each strand are tied at the top and bottom, respectively as shown in Fig. 11.21 

(b) . The multiplication of two braids 61 and 62 are defined by connecting the bottom and 
top endpoints of two braid diagrams as shown in Fig. 11.31 (a). The braid group Bn is 
generated by n — 1 generators {ai} subject to 

iTjiTj = o-jCTi for \i — j\ >2 (1-61) 

— (7 i(T . (1.62) 

The generator corresponds to the braid diagram with only one crossing of fth and z + lth 
strands as shown in Fig. 11.31 fbl. Then, the latter equality is nothing but the Reidemeister 
move III. By introducing an appropriate closure, which closes the endpoints of the strands, 
a braid diagram is related to a link diagram as we will see later. 

In order to construct a unitary representation of Bn, we embed it into the Temperley- 
Lieb (TL) algebra TLn{d) [TL71] on tangle diagrams with no crossing as shown in Fig. 11.21 

(c) . Similarly to the previous case, TLn{d) is generated by a set of generators {Ei ,..., En-i} 
subject to 

EiEj = EjEi for \i — j\ > 2, (1.63) 
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Figure 1.4: (a) The Markov trace closure for the braid b. (b) The plat closure for the 
braid b. (c) The plat closure is deformed into the Markov trace closure. 


EiEi±iEi = Ei, (1.64) 

Ef = dEi. (1.65) 

The generator Ei corresponds to U and H for the zth and (z + l)th strands as shown in 
Fig. 11.31 (c). The braid group is embedded into the TL algebra by 

ai = AEi + A-^l. ( 1 . 66 ) 

In order to relate the braid diagram and the Jones polynomial, we employ a Markov 
trace closure on the tangle diagram as shown in Fig. 11.41 (a). The Markov trance closure 
denoted by mtr has the following property: 

mtr(l) = 1, (1.67) 

mtr(74il) = miT{BA) for A,B G TLn{d), (1.68) 

mtr(j4) = dm.tT{AEn-i) for A E TLn-i{d). (1.69) 


Now the representation of the braid group and the Jones polynomial are related. Eq. 
(I1.66P corresDonds to suDeroosition of two smoothines { of a crossing'^)><^ . More 

precisely, let p and p be a representation of the TL algebra and an induced representation 
of the braid group, respectively. Moreover, the matrix trace has the same property as the 
Markov trace closure. Thus, for a given braid diagram 6, we have 

Ffe.n(^-^) = ATr[p(6)], (1.70) 
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where 6™*“^ is a link diagram generated from the braid diagram b with the Markov trace 
closure, and Amtr = (— 

Let us construct a representation of the TL algebra, the so-called path-model represen¬ 
tation. Suppose G is a one-dimensional graph with k vertices labeled by 1, ...,k from left 
to right and k — 1 edges. We consider an n step walk on G starting from the left endpoint. 
Let p = {vi,...,Vn) be a path of such an n step walk, where ui = 1 and Vj = 0 and = 1 
mean moving to the left and right neighboring vertices, respectively. Then we define a 
Hilbert space 'Hn,k spanned by all possible paths as a basis {|p)}. For each generator Ei, 
we define its representation <l>j = p{Ei) on 'Hn,k as follows: 


^>j|...Uj_i01uj+i...) 


0, 


(1.71) 

^Zi — 1 1 ^-1 

|...Ui_i01ui+i.., 

^ I V^Zi+lXz^-l I 

...Uj_il0ui+i...), 



(1.72) 

■^.22+ 1 1 1 

|...Ui_il0ui+i.., 

\ \/ ^Zi + lXzi-1 , 

■ A^. ' 

...Uj_i01ui+i...), 



(1.73) 

0, 


(1.74) 


where Zi G {1,..., k} is the label of the vertex at the ith step and Xj = sin(j0) with 0 = Trjk. 
The representation p{ai) = A^i + A~^l for the braid group Bn is induced by p. It is easy 
to confirm that, with A = and d = 2cos0, the representation p is hermitian for all 

Ei, and hence the induced representation p is unitary. Since p{<Ji) is a unitary operator 
acting on the neighboring two-qubit, a multiplication of such unitary operators can be 
implemented by using universal quantum computer. Since, by using the Hadamard test in 
Sec. 11.5.11 we can evaluate the Jones polynomial 14mtr(H“^) with an additive error A2”'e 
taking poly(n, m, /c, 1/e) overhead, where m is the number of the crossing. 

The approximation scale is further improved by considering another closure, the so- 
called plat closure as shown in Fig. 11.41 (b). We deform the link diagram obtained by the 
plat closure to another link diagram obtained by the Markov trace closure as shown in 
Fig. 11.41 (c). Then, we have U and H for all strands. In this case, the support of p{b) is 
only on 110101...), which allows us to replace the matrix trace as follows: 


V/pit(H-4) = Apit(10101...|p(6)|10101...), (1.75) 

where 6^^* is a link generated from the braiding b with the plat closure and Apit = 
(_^)3tj(6P'L^n,/2-i^ (Note that for each p{Ei) of re/2 Us and Us, a factor A 2 /A 1 = d is 
substituted.) Again using the Hadamard test, we can estimate 14pit(A“^) with an addi¬ 
tive error Ap^e taking poly(re, rre, fc, 1/e) overhead. 

Finally, we show that the approximation of the Jones polynomial with an additive 
error Apite is BQP-hard. Unfortunately, the computational basis defined by the path 
p = {vi,...,Vn) is inappropriate for this purpose, since the re-step walks cannot span the 
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whole 2”-dimensional Hilbert space. Instead, we employ the four-step encoding [WYOh] , 
where the computational basis states are defined by the following two four-steps: 

|0) = |llOO), |I) = |1010). (1.76) 

The space 7^4n,fe spanned by the 4n-step walks can support the whole 2"'-dimensional space 
for {|0), II)}®”. The dimension of the space spanned by the 8-step walks starting from 
vertex 1 ending at vertex 1 is 14. Thus the unitary representation p of the braid group Hg 
results in unitary operators in SU{1^). If they are dense in S'17(14), we can approximate 
an arbitrary unitary operator in SU (4) spanned by two four-step encoded qubits by using 
the Solovay-Kitaev algorithm |AA06| . which are enough to implement universal quantum 
computation. The density is achieved by fc > 4 and k ^ Q. For such a parameter, the 
approximation of the Jones polynomial with an additive error Apite is enough to solve a 
BQP-complete problem. 

The approximation of the Jones polynomial can be done by universal quantum com¬ 
puter and also is enough to solve the problems solvable by universal quantum computer. 
Thus the approximation of the Jones polynomial is a BQP-complete problem [AJLDQt 
IAA06] . The AJL algorithm for approximation of the Jone polynomial was extended for 
the Tutte polynomial in Ref. [AAELOT] . where the Solovay-Kitaev algorithm for non¬ 
unitary linear operators was developed. 


1.6 Quantum noise 


Quantum coherence, one of the essential properties of quantum systems, is quite fragile 
against noise, due to interactions between the system and the environment. Suppose that 
the system S of interest interacts via a unitary operation U with the environment E, 
where the system and environment are initially uncorrelated. The reduced density matrix 
p'g of the system after the interaction is calculated to be 


p's = Trs 


U{ps ® Pe)U'^ 


(1.77) 


Using a spectral decomposition of the initial state in the environment, pE = YlikPk\^k){^k\-, 
we obtain a map of the system S: 


where 


Ps — ,k)PsK]^i.pky 

{k\k) 


(1.78) 


K{k^k') = V»(efc|t^|efc)- 


(1.79) 


This map satisfies 

^lk',k)^(k',k) 

k 


Pk{ek\U^\ek){e'k\U^\ek) 


Is, 


(1.80) 
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where Is is the identity operator in the system S. 

In general, a map of a quantum state is given as a completely-positive-trace-preserving 
(CPTP) map £, subject to 

• Tt[£p\ = 1 for any density matrix p (preservation of probability) 

• (convex linear map) 

• 1a®£s{pas) > 0 (complete positivity) for an arbitrary ancilla system A and density 
matrix pAs on the composite system 

Such a CPTP map can always be written by Kraus operators {KTfe} |NC00] : 

M 

£p = Y,KkpKl (1.81) 

k=l 


Under the Born-Markov (with rotating-wave) approximation, the time evolution of 
a two-level system coupled with an environment is given by a master equation of the 
Lindblad form |PB02] : 

p{t) = -^[o-_c7+p(t) -F p{t)a-a+ - 2a+p{t)a-] 

— i^[<^+<^-p{t) + p{t)(^+(y- - 2cT_p(t)cj+] 

-'^['^zp{t) + p{t)(Tz - 2azp{t)(Tz\ = Cp{t), (1.82) 

where cr+ = |0)(1|, (t_ = |1)(0| and 7 q (a = 0,-|-, —) are the decay rates of the decay 
channels. One can easily find the eigenoperators of the Lindblad super-operator. These 
eigenoperators form the damping basis |BE93| : 


Cai = = Aio-i, (1.83) 

Ca, = + + (1.84) 

C-Cfz = (7+-F7_)cr3 = A3CJ3, (1.85) 

Cpeq = 0, (1.86) 

where peq = (7+|0)(0|-h7-|l)(l|)/(7-h+7-) = ((To-ha<T3)/2 with a = (7+-7_)/(27+-h27_). 
The solution of this master equation is given by the CPTP map £{t): 


£{t)p = pQ{t)p + ^ pi{t)aipai +f{t){a:ip +pa2,-i(Tipa2 + icr2p(Ti), 

1 = 1 , 2 ,3 

(1.87) 


where 

Po{t) = l(l + e-^i* + e-^=* + e-"3*), 


( 1 . 88 ) 


26 










pi{t) = i(l + e-^i*-e-^=*-e-^3*), (1.89) 

P2it) = ^(l-e-^i' + e-^=*-e-^3*), (1.90) 

p^{t) = l(l-e-^i*-e“'^2t_^g-A3t)^ ^191) 

fit) = ^(l-e-^3i). (1.92) 

If we consider a high temperature case (i.e., a —>■ 0), Eq. (|1.87p can be rewritten as 

3 

£{t)p = [1 -Pi{t) -P2it) - P'iit)] p +'^Piit)aipai. (1.93) 

i=l 


Hence, the CPTP map can be viewed as a stochastic Pauli error with probabilities Piit). 
In general, noise cannot be written by a stochastic Pauli error. However, by performing 
an appropriate operation, one can depolarize the CPTP map into a stochastic Pauli error 
as a standard form, where the noiseless part of the evolution is not altered [DHCBO^ . 
Otherwise, the Pauli basis measurements can collapse the CPTP map into a stochastic 
Pauli error, as we will see later. In the rest of this book, therefore, we will consider 
Markovian stochastic Pauli errors only. 

The fault-tolerance against more general noise has been discussed in Refs. [AB097( 
IABO081IAGP061IAGP081IAP09] . If the decoherence is non-Markovian, the dynamical de¬ 
coupling or quantum Zeno effect can be used to suppress the decoherence |Zur84[ IVGW961 
IDC981 IVL981 IVKL99al lFTP+n5| . Besides, if the decoherence is spatially correlated, 
one can utilize a passive error-prevention scheme, the so-called decoherence free subspace 
(DFS), which is immune to collective noise |ZR97[ IDG971 ILCW98b] . 
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Chapter 2 

Stabilizer formalism and its 
applications 


In general, the description of quantum states is a difficult task because it requires expo¬ 
nentially many parameters in the number of qubits as shown in Eq. (jl.38p . To understand 
these complex quantum systems, it is essential to have efficient tools. The stabilizer for¬ 
malism is one such powerful tool to describe an important class of entangled states. It also 
provides a diagrammatic understanding of quantum states and operations. The stabilizer 
states, described by the stabilizer formalism, play important roles in quantum computa¬ 
tion, such as for quantum error correction codes and resource states in MBQC. In this 
chapter, we introduce the stabilizer formalism, especially focusing on its diagrammatic 
understanding. Based on the stabilizer formalism, we explain quantum error correction, 
magic state distillation, and MBQC. 

2.1 Stabilizer formalism 

We first define an n-qubit Pauli group Vn'- 

iP„:={±l,±z}x{/,X,y,Z}®’^. (2.1) 

An element of the Pauli group is called a Pauli product. For example, the two-qubit Pauli 
group is given by 

V 2 := {±l,±f} 

X {II, IX, lY, IZ, XI, XX, XY, XZ, YI, YX, YY, YZ, ZI, ZX, ZY, ZZ}{2.2) 

where A0 B is denoted by AB for simplicity. (We will frequently use this notation when 
there is no possibility for confusion.) Next, we define an n-qubit stabilizer group S as an 
Abelian (commutative) subgroup of the n-qubit Pauli group: 

5:={5i}s.t. ^ 5 and'^5i,5j G = 0. (2.3) 
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Because —I is not included in the stabilizer group, all elements are hermitian Si = Sj, 
which guarantees that the eigenvalues = ±1. An element of the stabilizer group is called 
a stabilizer operator. The maximum independent subset Sg of the stabilizer group is 
called stabilizer generators. Here, independence means that any element of Sg cannot be 
expressed as a product of other elements in Sg. Any element of the stabilizer group can 
be generated as a product of the stabilizer generators. The stabilizer group S generated 
by the generators Sg is denoted by 5 = {Sg). 

Let us, for example, consider a two-qubit stabilizer group: 

5Beii = {II, XX, ZZ,-YY}. (2.4) 

Because they contain two anticommuting Pauli operators, XX and ZZ commutes. The 
stabilizer group ^Beii is generated by {XX, ZZ}, because —YY can be expressed as a 
product of XX and ZZ. Thus, we can write ^Beii = {{XX, ZZ}). 

For a given stabilizer group S, the stabilizer state is defined as a simultaneous eigenstate 
of all stabilizer elements Si ^ S with the eigenvalue -|-1: 

G 5, Si\ij) = IV^). (2.5) 

It is sufficient that the state is an eigenstate of all stabilizer generators: 

^SiGSg, 5i|V') = |V'). (2.6) 

Let k be the number of elements in the stabilizer generator Sg. Each stabilizer generator 
divides an n-qubit system (Hilbert space) into two orthogonal subspaces associated with 
the eigenvalues ±1. Because all stabilizer operators commute with each other, the k 
stabilizer generators divide the n-qubit system into 2^ orthogonal subspaces. Thus, the 
dimension of the space spanned by the stabilizer states, which we call a stabilizer subspace, 
is 2'^ = When n = k, we can define the quantum state uniquely. The number of 

stabilizer generators is at most n for an n-qubit stabilizer group. In the case of < n, the 
degrees of freedom in the stabilizer subspace can be addressed by using logical operators, 
which commute with all stabilizer generators and also are independent of them. 

Let us consider the stabilizer group ^Beii again. The stabilizer state is the eigen¬ 
state of XX and ZZ with eigenvalue -|-1, and hence given by the Bell state (|00) -|- 
|ll))/\/2 [EPR,35] . If XX is removed from the generators, the two-dimensional subspace 
spanned by |00) and |11) is stabilized. By choosing logical operators Lx = XX and 
Lz = ZI, we can specify the state in the subspace. For example, the eigenstate of Lx 
with the eigenvalue -|-1 is the Bell state. The eigenstate of Lz with the eigenvalue -|-1 is 
100). Another representative example of the stabilizer states is an n-qubit cat state, 

|cat) = ^(|00...0) + 111...1)), (2.7) 

\/n 

whose stabilizer group is given by 

l^iZ2, ..., Zn-lZn, • (2-8) 
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Figure 2.1: The stabilizer formalism: a Heisenberg picture of quantum computation. 
A Clifford operation is represented as a transformation of the stabilizer group by the 
conjugation of U. 

The cat state is a representative example of a macroscopically entangled state. If it is 
determined whether a particle is |0) or |1), the superposition is completely destroyed. If an 
element 0^=1 removed from the stabilizer generator, it defines a stabilizer subspace 
spanned by |00..0) and |ll..l). We can choose Lx = 0?=!"^* as logical 

operators, which anti-commute with each other and behave as logical Pauli operators. 

2.2 Clifford operations 

In the stabilizer formalism, we can describe a restricted class of unitary operations, the 
so-called Clifford operations, acting on the stabilizer states quite efficiently. The Clifford 
operation is defined as an operation U that transforms a Pauli product into another Pauli 
product under its conjugation, [...] —)• U[...]U'^. Let us consider the action of a Clifford 
operation U on the stabilizer state IV’) dehned by a stabilizer group S = 

U\^) = USi\il;) = USiUM^) = S'C/IV'), (2.9) 

where we dehne S[ = USiU"^. The above equality indicates that the state U\il^) is an eigen¬ 
state of the operator S[ with an eigenvalue -|-1 for all S[. Because 17 is a Clifford (unitary) 
operation, the group {5'} is also an Abelian subgroup of the Pauli group. Accordingly, the 
state 171V') is a stabilizer state with respect to the stabilizer group {S''}. In this way, the 
action of U on the stabilizer state can be represented as a transformation of the stabilizer 
groups under the conjugation of U as shown in Fig. 12.11 For example, the stabilizer state 
stabilized by {X 1 I 2 , 11 Z 2 ) is |-|-)i|0)2. The stabilizer group is transformed by A(77)i^2 into 
{X 1 X 2 , Z 1 Z 2 ), whose stabilizer state is (|00) -|- |ll))/\/2. 

The stabilizer formalism corresponds to the Heisenberg picture of quantum computa¬ 
tion, where a minimum number of operators are employed to describe a restricted type of 
quantum states and operations |Got97[ [Got98a] . This representation is powerful because 
it requires us to keep a time evolution of at most n operators, while a straightforward 
state-based approach needs exponentially many states. For example, let us consider the 
following quantum circuit: 
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A straightforward calculation yields the output state |V’), 

1^) = (jOOOOOOO) + IIOIOIOI) + lOllOOll) + IllOOllO) 

+ 10001111) + llOllOlO) + lOllllOO) + IllOlOOl) 

+ 11111111) + lOlOlOlO) + llOOllOO) + lOOllOOl) 

+ 11110000) + lOlOOlOl) + llOOOOll) + |0010110))/4. (2.10) 

It is rather cumbersome to write down the above state. Instead, we can understand the 
output state as a stabilizer state whose stabilizer generators are 

{ZIZIZIZ, IZZUZZ, UIZZZZ, XXXIIII, XXIXXII, IXIXIXI, XIIXIIXl2.ll) 

Equivalently, we may also choose the following stabilizer generators because they generate 
the same stabilizer group: 

{ZIZIZIZ, IZZIIZZ, IIIZZZZ, XXXXXXX, IIIXXXX, XIXIXIX, IXXIXXIp.U) 

Actually, these stabilizer generators are enough to understand the properties of the quan¬ 
tum state IV’). If an explicit description of the state is required, we can systematically 
write it down as follows: 

IV,) = 4 / + 54 I + 53 J + 52 / + 10000000^^ ( 2 . 13 ) 

where 5i = XIXIXIX, S 2 = IXXIIXX, S 3 = IIIXXXX, and S 4 = XXXXXX. 

The above equation means that jOOOOOOO) is an eigenstate for all Z’s stabilizer operators. 

By projecting it into the +1 eigenstate of the stabilizer generator Si by the projection 
we obtain the stabilizer state jV’)- 

In order for the above calculation to work, we have to obtain the stabilizer generators 
of the output state. This can easily be done graphically. We introduce commutation rules 
between the Pauli operators and Clifford operations below. In the case of the Hadamard 
operation, HX = ZH and ZH = HX, and hence we have 
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X.z 

- IL - 

^.^ X 

meaning that the Pauli X operator acting before the Hadamard operation is equivalent 
to the Pauli Z operator acting after the Hadamard operation and so on. Similarly, for the 
phase operation X, we have 


X. 

- A 

z . 


Y 


Z 


The CNOT operation transforms the Pauli operators under its conjugation as follows; 


A,,t{X)X,A,^t{X) = X,Xt, (2.14) 

A,,t{X)XtA,^t{X) = Xt, (2.15) 

A,,i(X)ZeAe,t(X) = (2.16) 

A,,t{X)ZtA,^t{X) = Z,Zt. (2.17) 


The commutation relation between the CNOT operation and the Pauli operators is un¬ 
derstood as follows; 
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In the above circuit diagram, the solid circle commutes with the Pauli Z operator, while 
the Pauli X operator is propagated as the Pauli X operator on the target qubit, making 
a correlation. Similarly, the open circle commutes with the Pauli X operator, while the 
Pauli Z operator is propagated as the Pauli Z operator on the control qubit, making a 
correlation. By recalling that the CNOT operation is transformed into the CZ operation 
by the Hadamard operations on the target qubit, the commutation relation between the 
CZ operation and the Pauli operators are obtained straightforwardly. This is described 
graphically as follows; 
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In this case, note that the Pauli X operation is propagated as the Pauli Z operation. 

This graphical understanding allows us to calculate the stabilizer generators of the 
output of the Clifford circuits. For example, in the following circuit diagram, the first 
qubit is stabilized by X before the Clifford operation. The Pauli X operator is propagated 
toward the right, and we obtain the stabilizer operator ZIZIZIZ for the output: 



The reader should use this graphical technique to calculate the other stabilizer generators 
and verify Eq. (I2.11|] . 

2.3 Pauli basis measurements 

Next, we will see how the Pauli-basis measurements on the stabilizer states are described 
in the stabilizer formalism. Suppose the A-basis {A = X, T, Z) measurement is performed 
on a stabilizer state IV'), whose stabilizer group is given by {Si). (We assume that the 
number of stabilizer generators is equal to the number of qubits, and hence that the 
stabilizer state is uniquely defined.) Depending on the stabilizer group ({S'*}) and A, 
there are two possibilities: 
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(i) The Pauli operator A commutes with all stabilizer generators. In that case, either A 
or —A is an element of the stabilizer group. If A (—A) is an element, the eigenvalue 
+1 (—1) is obtained with probability 1. The post-measurement state is the same as 
the stabilizer state before measurement. 

(ii) At least one stabilizer operator does not commute with A. In this case, we can 
choose another set of generators {S*^} such that S[ anti-commutes with A but all 
other generators commute with A. The measurement outcomes -|-1 and —1 are ob¬ 
tained with an equal probability of 1/2. The post-measurement state is given by 
((—^ 2 ,..., S'j^) depending on the measurement outcomes m = 0,1 correspond¬ 
ing to the eigenvalues (—1)"*. 

For example, suppose we perform the T-basis measurement on the first qubit of the 
Bell state stabilized by ^Beii = {XX, ZZ). We can redefine the stabilizer generators by 
{XX, —YY}. Then the stabilizer group after the measurement is given by {YI, —YY) = 
(yj, —lY). Thus, we obtain | — i) as the post-measurement state on the second qubit. 

2.4 Gottesman-Knill theorem 

Because the stabilizer states and Clifford operations are described efficiently in the sta¬ 
bilizer formalism, it implies that such a restricted type of quantum computation can be 
simulated efficiently on a classical computer. This is stated by the Gottesman-Knill theo¬ 
rem [Got971 IGot98a[ INCOO] . 

Theorem 1 Any Clifford operations, applied to the input state |0)®” followed by the Z 
measurements, can be simulated efficiently in the strong sense. 

Here, an efficient strong classical simulation of a quantum circuit C is a classical polynomial¬ 
time computation that calculates the probability Pc{x) for a given output x of the circuit 
C, including an arbitrary marginal distribution Pc{x). (See, for example. Ref. [B.ISll] 
for the definition of a strong simulation.) Note that this theorem holds true even when the 
initial state is generalized to an arbitrary stabilizer state, and also any Pauli products are 
measured, because they are done in the above setup by modifying the Clifford operations 
appropriately. 

Proof: The stabilizer group of the input state is {{Zi}) {i = 0,1, ...,n — I). By apply¬ 
ing the Clifford operations as mentioned, we obtain the stabilizer generators ({S’i}) of the 
quantum output before the measurements. Suppose the measurement outcome, the clas¬ 
sical output, is given by {mj = 0,1}. Then the probability of obtaining the measurement 
outcome {rui} can be calculated as follows: 

i) Set the stabilizer generators = ({S'i}) and the initial probability = 1. 

ii) For A: = 0,1,...,n — 1, repeat the following procedures. 
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1) If (—G update the probability because the measure¬ 

ment outcome is obtained with probability 1. The stabilizer group after the 
measurement is also updated to 

2) Else, if (—G update the probability = 0, because such a 

measurement outcome does not appear. (You may stop the calculation at this 
stage, and return the probability 0.) 

3) Else, is updated into by removing an anticommuting generator and 

adding (—as a new generator. Because the measurement outcome is ob¬ 
tained randomly with probability 1/2, the probability is taken as = p(*')/2. 

iii) Return as the probability of obtaining the measurement outcome {ruj}. 

Note that, in step (ii), we can efficiently decide which of the three is the case for any k by 
checking the commutability of Zj. with the stabilizer generators of . ■ 

The statement of Theorem [1] can be extended by weakening the notion of the classical 
simulation. 

Theorem 2 Any Clifford operations, applied to any product states of convex mixtures 
of the Pauli basis states, followed by Z measurements can be efficiently simulated in the 
weak sense. 

Here, an efficient weak classical simulation of a quantum circuit C is a classical polynomial¬ 
time randomized computation that samples the output x according to the probability 
distribution Pc{x) of the output of the circuit C. (See, for example, Ref. |BJS11| for the 
definition of weak simulation.) Apparently, a strong simulation includes a weak simulation, 
because we sample the output by using the marginal distributions [TDOdj . 

Proof: Suppose that the Rh input qubit is given by 

Pi = pi*Vl+)(+l +pt]-\-){-\ +4hI + +4*-l “ 

+p«+|0)(0|+pi*L|l)(l|, (2.18) 

where Yla=xy z Yl,v=+ - P^a,u = 1- By using the probability distribution the input 

state of each qubit is randomly sampled. Conditioned by the sampling result, the input 
state is a product of the Pauli basis states, and hence the output probability distribution 
can be calculated as shown in Theorem [TJ Combined with the random sampling of the 
input state, this provides an efficient weak simulation of the Clifford circuit with noisy 
input states (convex mixture of the Pauli basis states). ■ 

The input state can be generalized into a classical mixture of stabilizer states, when 
its polynomial size description of the probability distribution is provided. Similarly, the 
Clifford operations can be extended to stochastic Clifford operations such as the stochastic 
Pauli error. 

The convex mixture of the Pauli basis state lies inside the octahedron of the Bloch 
sphere as shown in Pig. 12.21 It is natural to ask whether or not the Clifford circuit allows 
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Figure 2.2: A convex mixture of the Pauli basis states lies inside the octahedron of the 
Bloch sphere. 


universal quantum computation if the input state lies outside the octahedron. If the input 
state is a pure non-stabilizer state such as we can implement a non-Clifford 

gate by using gate teleportation, explained in Sec. 12.61 Even some mixed states 

can be converted into a pure non-stabilizer state, the so-called magic state, by using only 
Clifford operations. Such a protocol is called magic state distillation [BK05j and will be 
explained in Sec. 12.81 


2.5 Graph states 


In this section, we introduce an important class of stabilizer states, the so-called graph 
states |HDE+n6] , whose stabilizer generators are defined on graphs. The graph states are 
employed as resource states for MBQC as explained in the next section. 

A graph state is defined by a graph G = {V,E). Here, V and E are the sets of the 
vertices and edges, respectively. A qubit is located on each vertex of the graph. The 
stabilizer generator of the graph state |G) is defined as 


Ki = XiY\ Zj for all i G E, 

j&Vi 


(2.19) 


where we define a set of vertices V) := {j\{i,j) £ E}, which are connected to the vertex i 
by an edge on the graph G (see Fig. 12.21) . The graph state |G) is generated from a product 
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G=(V,E) 



Figure 2.3: The graph state |G) associated with a graph G = {V,E). A stabilizer generator 
Ki is also shown. 

state by applying the CZ gate on each of the graphs: 

|G)= n (2.20) 

(iJ)eE 

where |F| indicates the number of vertices of the graph G = {V,E). This can be un¬ 
derstood that the stabilizer generator Xi for the state |-|-) is transformed into Ki by the 
CZ operations U = K[Z)i^j. Especially, when the graphs are regular lattices such 

as one-dimensional (ID), square, hexagonal, and cubic lattices, the corresponding graph 
states tend to be referred to as cluster states |BR01| . Any stabilizer state is equivalent 
to a certain graph state up to local Clifford operations [VdNDDMn4( lHDE~*~n6j . For ex¬ 
ample, the cat state is equivalent to the following graph state by applying the Hadamard 
operation on the /cth qubit: 


k 



Unfortunately, the graph associated with a stabilizer state is not uniquely defined, because 
there are local Clifford operations that change the underlying graph. This property is 
called the local complementarity of the graph states [VdNDDM04[ |HDE'*~0^ . 

Next, we will see how the Pauli basis measurements transform the graph states. For 
simplicity, we assume that the state is projected into an eigenstate with eigenvalue -|-1. 
Let us consider a ID graph state as follows: 

<KXXK> 

i-1 i i+1 
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whose stabilizer generator is given by 


Ki = Zi-iXiZi+i. ( 2 - 21 ) 

We first consider the Z basis measurement (projective measurement of the observable Z) 
on the zth qubit. Following the procedure seen in Sec. 12.31 Ki is removed from the stabilizer 
generator. By adding Zi instead, we obtain the stabilizer group for the post-measurement 
state 


{...,Ki_i,Zi,Ki+i,...). ( 2 . 22 ) 

After the projection, the fth qubit is |0) and hence decoupled from the other qubits. By 
rewriting the stabilizer generators, we obtain three decoupled stabilizer groups 

(..., Zj_2Ai_i), {Zi), (Aj+iZj+2,...). (2.23) 

This means that the graph is divided into two parts as follows: 

-► O O 

1-1 I I+1 i-i i+i 

For any graph, this property of the Z-basis measurement holds; the post-measurement 
state is defined by a modified graph, where the vertex corresponding to the measured 
qubit and the edges incident to it are removed from the original graph. 

Next, we consider the W-basis measurement. The observable Xi does not commute 
with ATj-i and ATi+i, but does commute with iFj_iiFj+i = Zj_ 2 A'j_iX 4 +iZj_|_ 2 . Following 
the procedure in Sec. 12.31 the stabilizer group for the post-measurement state is calculated 
to be 


(...,Zi_2W_iW+i Zi+2^ Zi-iZi^i, ...), {Xi). (2.24) 

By performing the Hadamard operation H on the (i —l)th qubit, we obtain a new stabilizer 
group 


(..., A'j_iZj+i,...), {Xi), (2.25) 

which indicates that the graph is transformed into the following graph with the Hadamard 
operation: 

0 ^ 2)-0 -►-(> = 

i-1 i i+1 i+1 
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Instead of the (i — l)th qubit, we can obtain a similar result by performing the Hadamard 
operation on the (i + l)th qubit as shown above. 

Suppose the ith and (i + l)th qubits are measured in the X-basis on a ID graph state. 
This is equivalent to measuring the (i + l)th qubit of the above post-measurement graph 
state in the Z basis, because the Hadamard operation is applied on it as a byproduct. 
From the previous argument, the Z-basis measurement remove the measured qubit from 
the graph. Thus, two neighboring X-basis measurements remove the measured qubits and 
connect the left and right hand sides directly: 

-(XE>-(E>-0->■-O -O 

which we call a contraction. 

Finally, we consider the T-basis measurement. The observable Yi does not commute 
with either Ki-i, Ki, or Xj+i, but does commute with Ki-iKi = Zi- 2 Yi-iYiZi^i and 
XjXj_|_i = Zi-iYiYi^iZi^ 2 - The stabilizer group for the post-measurement state is calcu¬ 
lated to be 

(..., Zi—2Yi—iZi^i, Zi—iYi^iZi^2 ) •• •)) (Yi). (2.26) 

By performing the phase gates S on the (i — l)th and (i -|- l)th qubits, we obtain a new 
stabilizer group 

(..., Zj_2Xi_iZj+i, Zj_iXj+iZj+2, ■■■), {Yi). (2.27) 

This indicates that the graph is directly connected up to the phase operation 5 as a 
byproduct: 


- 0 -( 2>-0 -©- 

i-1 i i+1 i-1 i+1 

Suppose three neighboring qubits (i — 1), i, and (i -|- 1) are measured in the T-basis. 
This is equivalent to measuring the ith qubit in the T-basis, and then measuring the 
(i — l)th and (i -|- l)th qubits of the post-measurement graph state in the X-basis, because 
there is a phase operation S acting on them as a product. As seen previously, the X-basis 
measurements on two neighboring qubits result in a contraction of the two qubits on the 
graph. Thus, the X-basis measurements on three neighboring qubits contract them from 
the ID graph state. 

TMEHrHrKT-►-O-O 

This property is useful to change even and odd of the length of the ID graph state. 

While we have considered the Pauli-basis measurements only on the ID graph state, 
we can generalize these arguments into graph states of general structures. A graph state 
is still mapped into another graph state up to some single-qubit Clifford operations as 
byproducts. 
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2.6 Measurement-based quantum computation 


Measurement-based quantum computation (MBQC) is a model of quantum computation, 
where quantum gates are implemented by adoptive measurements on a highly entangled 
resource state [RBOll IRBB03t IBauPd] . Specifically, certain graph states, the so-called 
cluster states, are employed as resource states in MBQC. Below we will first demonstrate 
quantum teleportation, a building block of MBQC. Then, we explain how adoptive mea¬ 
surements on a graph state enable us to emulate universal quantum computation via 
quantum teleportation. 

Quantum teleportation is a quantum communication protocol, in which Alice sends 
a quantum state to Bob by using a shared entangled state and classical communica¬ 
tion BBC'*~93] . Suppose Alice and Bob share a maximally entangled state, the Bell state, 

|0)„|0)6 + |l)a|l)fe 




(2.28) 


For an unknown input state and the half of the Bell state, Alice performs a Bell 
measurement, which is a projection onto the Bell basis states 


(2.29) 

where mi, m 2 = 0,1 correspond to the measurement outcomes. A straightforward calcu¬ 
lation provides 


{^{mi,m2)ka = ZpX^^\i;),/2. (2.30) 


Hence, the unknown input state is teleported to Bob with a byproduct operator . 

If Bob does not know the measurement outcomes (mi, m 2 ), the teleported state is a 
completely mixed state for Bob. However, if Alice sends the measurement outcome as a 
classical message. Bob can undo the byproduct and obtain the unknown quantum state 
at Bob’s side. The circuit diagram of quantum teleportation is: 


Alice — 


Bob 


Bell measurement 



maximally entangled state 


By using the following circuit equivalence, we can decompose the teleportation circuit into 
two elementary teleportations, the so-called one-bit teleportations: 
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One-bit teleportation is useful as a building block of the teleportation-based gates em¬ 
ployed in MBQC. A single-qubit Z rotation can be implemented in a teleportation- 
based way. Its action can be understood from the following circuit equivalence: 


lO 

1 +) 

W 

1 +) 



where we utilized the fact that and A(Z) commute. The controlled-Z operation is 
also implemented in a teleportation-based way as follows: 
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{X2Zir-{XxZ2r^K2{Z)HiH2\i^)i,2 


in 1 

l+> 



mi 


in 2 

I+) 



m 2 



{X2Z^r^{XiZ2r^Xl.2{Z)H,H2\^)l.2 


That is, instead of performing the A(Z) gate after one-bit teleportations, we can pre¬ 
pare a special resource state, on which the X{Z) gate is pre-implemented, and the A(Z) 
gate is then performed via teleportation. These quantum operations based on quantum 
teleportation are called gate teleportation |GC99] . 

Now we are ready to formulate MBQC. An arbitrary single-qubit unitary operation U 
can be decomposed, up to an unimportant global phase, into 

U = He^^^e^^^e^^^ (2.31) 

= (2.32) 

This indicates that we can perform an arbitrary single-qubit unitary operation by a se¬ 

quence of one-bit teleportations. The resource state for the sequential one-bit teleporta¬ 
tions is a ID cluster state: 

where the stabilizer generator for the left-most qubit is removed and the input state 
is encoded. We have to take care of the byproduct Pauli operators depending on the 
measurement outcomes. Fortunately, we can propagate the Pauli byproduct operators 
forward as follows: 

JJ ^ ^ 2 . 33 ) 

_ j^mi+2tSrni ^rrii+i ^ ( 2 . 34 ) 
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IV’) 

1 +) 

1 +) 

1 +) 



j^mi+2©mi 2;mi+i It/^) 


By choosing O' = (—and (j)' = (—l)™'*+ii;/) adaptively, depending on the previous mea¬ 
surement outcomes, the random nature of the measurements can be managed. This proce¬ 
dure is called feedforward. The Pauli byproduct is propagated and updated throughout the 
computation. Note that the classical processing required to determine the measurement 
angle has only XOR (addition modulo two) operations [RBBOSj. 

Next, we will investigate the measurement-based two-qubit gate operation. The re¬ 
source state for the gate teleportation is the following cluster state: 


/• /+! 


(H>0 

im outi 


®Kj- 

ing out2 


ini 

I+) 


rt\ 

c 

= mi 

Uti 

l+^ „ 

(A2Zi)™RAiZ2)™^Ai,2(Z)//i//2iV')l,2 

ina 


c 

r/\ 

Ut2 


To adjust the timing of the two-qubit operation, we can insert identity operations depend¬ 
ing on the even and odd lengths as follows: 


inlSKEKfe, in(2HYKY>Oou, 


Without loss of generality, we can assume that all input states of the quantum computation 
are given by |-|-), which are automatically encoded by preparing the graph state. At the 
end of the computation (on the right-most qubits), measurements are performed to read 
out the result as follows: 
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In this way, universal quantum computation is simulated solely by measurements on a 
brickwork-like cluster state. This state can be generated from a cluster state on a square 
lattice by using the Pauli basis measurements as shown above. Accordingly, the square 
lattice cluster states are universal resources for MBQC. 

The above circuit-based explanation of MBQC [Nie06| is very intuitive and straight¬ 
forward. However, for a complicated resource state, as will be shown, an operator-based 
understanding of MBQC [RBB03j is quite useful. Let us reformulate MBQC from an 
operator viewpoint. Suppose again an MBQC on a ID cluster state. The measurements 
are executed up to the {i — l)th qubits, and hence the operator Ki {I < i) is removed from 
the stabilizer generators. The logical degree of freedom on the remaining resource state 
can be specified by the ith logical operators 


tW 


(2.35) 

T (*) 


(2.36) 

T (b 

= 

(2.37) 


These logical operators commute with all remaining stabilizer generators Ki {I > i + \). 
Moreover, they anticommute with each other, satisfying the commutation relations for the 
Pauli operators. Thus, they specify the state encoded in the graph state. As seen above, 
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a Z-rotation e jg applied before the X-basis measurement. Because Zi = this 

rotation induces a unitary transformation U of the logical operator 

cos + sin 0Ly\ (2.38) 

L^y COS 0Ly^ — sin (2.39) 

(2) (2+1) 

Because Lr^ = XiL^^ , the logical X operator after the X-basis measurements is given 
by (—depending on the measurement outcome rrii = 0,1. On the other hand, 
the logical operators Ly ^ do not commute with the X-basis measurement; they are not 
relevant logical operators after the measurement. If two operators are equivalent up to 
multiplications of the stabilizer operators, their action on the stabilizer state is also the 


same. By using this fact, we can replace the logical operators in (|2.39p with 

L® ~ L«Xi+i = Xi+iZi+2 = (2.40) 

4^ ~ Ki+i = XiYi+iZi+2 = XiLp^\ (2.41) 

where ~ indicates that two operators are equivalent up to stabilizer operators. After the 
X-basis measurement, Xj can be replaced by its eigenvalue (—l)™'^ Then the logical 
operator of the post-measurement state is given by 

4 ^ (-1)™*(cos0l 44 sin (2.42) 

4 ^ = (2.43) 

if 4"^^^ = (2.44) 


We now realize that the logical operators for the ith step are transformed into those for 
the {i -|- l)th step rotated hy U = 

Similarly, a two-qubit gate in MBQC can also be regarded as a propagation of a 
correlation by a projection on the stabilizer state. Consider the following graph state. 



stabilizer 


logical operators 
after the measurement 




Z 



The logical operators for the rth step are given by and mul¬ 

tiplying the stabilizer operator, we obtain 


L 

L 


(i) 

Zl 

(i) 

Z2 




Xi,,+iZi,i+iZ2,i+l = 4t'^42'^ 

X2,i+lZ2,i+lZ^,i+l = 
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(2.45) 

(2.46) 




The logical operators for the {i + l)th step after the projections are calculated to be 

{lSuS} ^ = {CLjt'Vt.yig+'Vt}, (2.47) 

T (b 
^X 2 ^^Z 2 


{Z,® .i®} ^ = {V'Lj^Vt.V'Lg'Vt}. (2.48) 


Again, we realize that the logical operators for the zth step are transformed into those for 
the (z + l)th step with a two-qubit unitary operation 


V = {XiZ2r^{X2Zir^ki^2{Z)HiH2. (2.49) 

By combining single-qubit rotations and the two-qubit operation {XiZ 2 )'^^ {X 2 Zi)^^ Ai^ 2 iZ)Hi. 

as seen above, we can perform a universal quantum computation. In this way, MBQC can 
be understood in the Heisenberg picture. 

Suppose the logical Pauli operators of the kth. input and output qubits are related by 
the measurements as follows: 


r r (In) 


r{In)4 
^Z,k S 


{UL^x!k^U^ 


Z,k 


(2.50) 


The unitary operator U is performed on the input qubits. Here a Pauli byproduct, de¬ 
pending on the measurement outcomes, is also included in U. Moreover, if two graph 
states, which perform U and V, are concatenated with the appropriate feedforwarding of 
the Pauli byproducts, then VU is performed: 


input 

output 


Of 

OH" 

OH 


Off 

of 

oJ.. 

oW.. 

OH 

jj- input 

output 


Let us consider the example shown in Fig. 12.41 lai. The logical operators on the inputs 
are replaced by multiplying stabilizer generators so that they commute with the X-basis 
measurements. Then the X operators on the measured qubits are replaced by ±1. The 
measurements transform the input logical operators as follows: 


rr(In) r(In) ^ (In) T'lln)-, 
\^X,l ’ -^Z,l ’ ^X,2 ’ ^Z,2 i 


(Bme ji^iOut) 2^7774©ms ^(Out) | 


(2.51) 


Thus, the ki^ 2 {Z) gate is implemented up to a Pauli byproduct. 

Using this fact and concatenation of the input-output relations, we can construct a 
measurement-based CNOT gate between the separated two-qubit as follows |RBB03j : 
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The input-output relation of the above graph state is equivalent to that for the following 
circuit: 


iritarget 


IHcontrol 

































OUttarget OUtcontrol 




In this way, a CNOT gate between two arbitrary separated qubits can be implemented 
using a constant depth (constant width) resource state. 

Topologically protected quantum computation was first proposed as a MBQC on a 3D 
cluster state |RHG06l IRHG07| as we will see in Ghapter [5l It was further applied to the 
circuit-based model in two dimensions (2D) |RH07( IFSGOQj . which will be explained in 
detail in Ghapter 01 Specifically, in the formulations, the operator-based understanding is 
quite useful for tracking how correlation is propagated in a topologically protected way. In 
Ghapter O we will formulate a topologically protected MBQG from the operator viewpoint 
and see how it is related to the circuit-based understanding explained in Ghapter 01 

To conclude, we summarize the properties of MBQG and recent progress in this area. 
A unique feature of MBQG is that the resource state for universal quantum computa¬ 
tion is prepared offline. Entangling operations, which would be one of the most difficult 
tasks in experiments, are employed only in this stage. Quantum computation is executed 
solely by adaptive measurements. This property is useful for experimental realization in 
certain physical systems. For example, a deterministic entangling operation is difficult to 
achieve in an optical quantum computation. In such a case, we can utilize linear optics 
and measurement-induced nonlinearity to generate a cluster state [YR.n31 INieOdl IBR,n5j . 
Importantly, the entangling operation can be nondeterministic, as long as the successful 
or non-successful outcome is heralded. By using such a probabilistic entangling operation, 
we can gradually expand the cluster state. After successful cluster state generation, we 
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can start the measurements for quantum computation. Note that the probability of suc¬ 
cessful cluster state generation is not exponentially small by using a divide and conquer 
approach |Nie04l IBR051IDR051IMBFIOI ILBSBlOi IFTlOj . 

The clear separation between the quantum stage requiring entangling operations and 
the measurement stage is useful, not only for the physical implementation, but also in 
a quantum cryptographic scenario. Suppose that Bob (server) possesses a fully fledged 
quantum computer and that Alice (client), who has a less advanced quantum device, such 
as a single-qubit state generator, wants to delegate quantum computation to Bob. By using 
the idea of MBQC, such a delegated quantum computation can be made unconditionally 
secure. This is called a blind quantum computation and was proposed by Broadbent, 
Fitzsimons, and Kashefi (BFK) [BFROQ] (see also the related earlier works |Chin5(lASn6j ). 
In the BFK protocol, Alice sends randomly rotated qubits to Bob, where the 

angle is chosen to be 9j = hjir/A {kj = 0,1,..., 7). Bob generates a cluster state by 
using the randomly rotated qubits. In the computation phase, Alice sends a classical 
message 6j = cj)j 9j + r^vr. Here, cj)j is the measurement angle with which Alice want to 
perform a measurement. The angle 9j is for the randomly rotated state (which is secret 
to Bob). The random bit r G {0,1} makes the measurement angle completely random 
for Bob. Then Bob performs the measurement in the {e*'^-’ |±)} basis. Because the initial 
state is pre-rotated by 9j (from Alice’s viewpoint), Bob performs the measurement in the 
|g*( 0 j+rj 7 r)Ij which is what Alice wants to do. However, from Bob’s viewpoint, 

the state is a completely mixed state with no information about {(pj}- Thus, Bob is 
blind to any information about the input, the algorithm, and the output. Instead of the 
state generation, Alice, who has a measurement device, can also perform a blind quantum 
computation, whose security is guaranteed by the no-signaling principle |MF13| . A fault- 
tolerant blind quantum computation has been proposed, based on topologically protected 
MBQC |MF12j . 

2.7 Quantum error correction codes 

In this section, we introduce stabilizer codes, which are a class of quantum error correction 
(QEC) codes. 

Three-qubit bit-flip code 

The QEC codes can be described elegantly in the stabilizer formalism. Let us first consider 
the simplest one, the three-qubit bit flip code, whose stabilizer generators are given by 

5i = Z 1 Z 2 , S 2 = Z 2 Z 3 . (2.52) 

The stabilizer subspace is spanned by the following two logical states: 

\0l) = |000), \1l) = |111). (2.53) 

The logical Pauli-A operator is given by Lx = A 1 A 2 A 3 . The logical Pauli-Z operator 
is defined as Lz = Zi. We may, equivalently, choose the logical Pauli Z operator to be 
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or Z 3 , because their actions on the code space are equivalent. The present code is a 
quantum analogue of the classical three-bit repetition code. Consider a bit flip error with 
an error probability p: 


£iP = (1 - p)p + pXipXi. (2.54) 

If the initial state = ckIOl) + /3|1l) undergoes the bit flip error independently, the 
output state is transformed in leading order as 

£1 0^2 o^3|V'L)(t^L| = (1 +P{^ + 0{p^). 

i 

(2.55) 

The error Xi maps the code space to an orthogonal space. We perform a projective 
measurement onto the orthogonal subspaces, = (/ ± Sk)/2, which we call a syndrome 
measurement, to know in which orthogonal space the state lies. Note that the encoded 
quantum information is not destroyed by the syndrome measurement, because it commutes 
with the logical operators. According to the measurement outcomes, the logical state can 
recover from the error as follows: 

7^o£:l o£:2oT3|V^l)(V’l| = [( 1 -p)^+ 3p(l- p)^]|V’l)(V’l|+ 0(p^), (2.56) 

where the recovery operator is given by 

Tip = PtP^pP^Pt+ X^P{P+pP+P{X^+X2P{P^pP^P{X2 

+X^P+P^pP^P+X^. (2.57) 

The four terms in TZp correspond to the measurement outcomes (eigenvalues) (-1-1, 4-1), 
(—1,4-1), (—1, —1), and (4-1,—1) of the stabilizer generators, respectively. By comparing 
Eqs. (j2.55l) and (I2.56p . one can understand that if p is sufficiently small, the fidelity of the 
logical state is improved. 

Similarly, we can construct a three-qubit phase flip code, which can correct a phase 
flip error, by changing the basis with the Hadamard transformation: 


(ZiZ2,Z2Z3)^(XiX2,X2A3). 


(2.58) 


9-qubit Shor code 

The three-qubit bit-flip code cannot correct Z errors, which commute with the stabilizer 
generators. A QEC code that can correct all A, T, and Z errors was developed by Shor 
based on a concatenation of three-qubit bit-flip and phase-flip codes |Sho95| . The stabilizer 
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generators of the 9-qubit Shor code are given as follows: 
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The code space is spanned by 

|0l) + |1l) ^ (loop) + 1111 )) (loop) + 1111 )) (loop) + | 111 )) 
^ 2\/2 

|0l)-|1l) ^ (|PPP) - |111)) (|PPP) - |111)) (|PPP) - |111)) 
^ 2\/2 


(2.59) 


(2.6P) 

(2.61) 


The logical Pauli operators are given hy Xl = and = Z®®, which are bitwise 
tensor products of physical Pauli operators. If the logical A operator is given by a bitwise 
tensor product of the physical A operators on the QEC code, we say that the operation A 
has transversality. The 9-qubit code is capable of correcting all X, Y, and Z errors for each 

qubit, which can be understood because the three-qubit phase flip code {|+++), |-)} is 

constructed by using the three logical qubits of the three-qubit bit flip codes {|PPP), |111)}. 

Note that any single-qubit noise £ can be described by using the Kraus operators {Kj}: 

£p = Y,K,pK]. (2.62) 

3 


Any operator Kj can be decomposed into the Pauli operators ao = I, ai = X, a 2 = Y, 
and = Z: 


Kj = '^Cjiai. (2.63) 

i 

Thus, if the X and Z errors on a single qubit are both corrected appropriately, we can 
correct any single-qubit noise automatically. Specifically, because noise contains a super¬ 
position of the Pauli errors, it can be collapsed by the syndrome measurements. 


Stabilizer codes 

To summarize the above examples, let us formalize the stabilizer quantum error correc¬ 
tion codes and their properties. The code space of a stabilizer QEC code is defined by 
a stabilizer group ({5*}). The encoded degree of freedom is specified by the mutually 
independent logical operators {T®}, which commute with all stabilizer generators and are 
independent of the stabilizer generators. The computational basis state of the code state 
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is completely determined by the stabilizer group ({Sj}, {(—1 )"*jL^}). We can always find 
another set of the logical operators {L^} being subject to 

LfLf = (2.64) 

where 6ij is the Kronecker delta. Hence, the pair of logical operators Lf and Lf represents 
the ith logical qubit. In terms of the numbers n and (n — k) of qubits and stabilizer 
generators, respectively, the number of pairs of logical operators is k. 

Let us define the weight wt(S'i) of a Pauli product Si as the number of qubits on which 
a Pauli operator (except for the identity 1) is acting. The minimum weight of the logical 
operator over all possible logical operators is called the code distance d. This implies that 
all Pauli products whose weights are smaller than d are elements of the stabilizer group 
or anticommute with the stabilizer generators. Thus, they act trivially on the code state 
or map the code state into an orthogonal subspace. If the weight of a Pauli product as 
an error is less than {d — l)/2, we can find a unique recovery operator that returns the 
erroneous state into the code space. Thus, we can correct weight-[(d — l)/2j errors. Such 
a stabilizer QEC code is called a [[n, k, d]] stabilizer code. For example, the code distance 
of the 9-qubit code is a [[9,1,3]] stabilizer code correcting weight-one errors. 

The nine-qubit code is not the smallest QEC code that can correct all weight-one X,Y, 
and Z errors. The smallest code is the five-qubit code, found independently by Laflamme 
et al. [LMPZ96] and Bennett et al. [BDSW9^ . The stabilizer generators and the logical 
Pauli operators are given as follows: 
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We see that the code distance is three, and hence an arbitrary single-qubit error can be 
corrected. 

Calderbank-Shor-Steane codes 

The readers who are familiar with classical coding theory might already be aware of the 
correspondence between stabilizer codes and classical linear codes. Let us recall the 9- 
qubit code. The X and Z errors are detected independently through the Z-type and 
X-type stabilizer generators, respectively. This implies that X and Z error corrections 
are described by classical coding theory, where two classical error corrections are subject 
to a certain constraint to appropriately form a stabilizer group. 

To formulate this, we briefly review classical linear codes. A [[n, A:]] classical linear 
code C is defined as a A-dimensional space Vc over GF{2^) by using an n x k generator 
matrix 

G= (bi,...,bfc), (2.66) 
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where the column vectors {bj are the basis vectors of V^- A fc-bit classical information 
y is encoded into the code c as 

a = Gc. (2.67) 

To detect and analyze the errors, we define an (n — k) x n parity check matrix H such 
that Hhk = 0 for all basis vectors {b^}. Suppose an error e occurs on the code state, 
a' = a 0 e, where 0 indicates a bitwise addition modulo two. By using the parity check 
matrix H, we can detect the error 

Ha' = H{a (B e) = He = s, (2.68) 


where s is called an error syndrome. 

For example, the three-bit repetition code is defined by the generator 


G = 



(2.69) 


A classical bit 0 and 1 is encoded into (0,0,0)"'" and (I,!,!)"*", respectively. The parity 
check matrix is defined to be 


H = 


1 1 0 \ 

0 1 1 y ■ 


(2.70) 


Now, we realize that the positions of the Is of the parity check matrix are exactly the 
same as those of the Zs in the stabilizer generators of the three-qubit bit flip code. This 
suggests to use the parity check matrices H^ and H^ of the two classical linear codes Gx 
and Gz, respectively, in the definition of the X-type and Z-type stabilizer generators: 

5^^ = , Sf = . (2.71) 

j j 

For these operators to commute with each other, the two parity check matrices have to 
satisfy 


HxHj = 0, 


(2.72) 


where 0 indicates a matrix with all elements = 0. To define the logical Z operators, we 
define a quotient space Keic{Hx)/VHz '■= {c\HxC = 0, c G where Vh^ is the 

space spanned by the column vectors of HJ. Denoting the basis vectors of the quotient 
space Kei{Hx)/VHz by {[b|]}, we define the logical Z operators 

L (2.73) 

i 

Similarly, we can define the logical X operators 

(2.74) 
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using the basis vectors {[b^]} of a quotient space /Vh^, where is chosen such 

that = (—Note that dimensions of these kernel subspaces are the 

same, and we can easily find such pairs of anticommuting logical operators. The above 
stabilizer code constructed from two classical linear codes is called a Calderbank-Shor- 
Steane (CSS) code. 

Let us see an important example of CSS codes, the 7-qubit code introduced by Steane |Ste96| . 
Specifically, we utilize a classical linear code, the [[7,4, 3]] Hamming code, whose generator 
and parity check matrices are given by 


G 


/ 1 0 0 0 \ 
0 10 0 
0 0 10 
0 0 0 1 
0 111 
10 11 
V 1 1 0 1 / 


H = 


1 0 1 0 1 0 1 \ 

0 110 0 11 
0 0 0 1 1 1 1 / 


(2.75) 


Because HH^ = 0, we can employ the Hamming code to dehne both X- and Z-type 
stabilizer generators: 
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There is an element (1,1,1,1,1,1,1)"*" in the quotient space Ker{H)/VH- The logical 
operators are given by 


Lx = XXXXXXX, 
Lz = ZZZZZZZ. 


The 7-qubit code is quite useful for fault-tolerant quantum computation. Both the X- 
and Z-type stabilizer generators are defined from the Hamming code, and the stabilizer 
group is invariant under the transversal Hadamard operation H = . Moreover, the 

logical X operator is mapped into the logical Z operator, HLxH = Lz- Thus, the 
transversal Hadamard operation acts as the logical Hadamard operation for the encoded 
degree of freedom. Similarly, a transversal phase operation S = (ZS)®'^ acts as a logical 
phase operation, SLxS^ = LxLy- Furthermore, a transversal CNOT operation A{X) = 
A(X)®^ keeps the stabilizer group of two logical qubits invariant: 

({5^^ (8) I®^}, {Sf (g) I®^}, {I®^ (g) }, {I®^ (g) Sf}) (2.78) 

= ({5« ® 5«}, {5« ® I®^}, {I®^ ® 5f}, ® 5«}). (2.79) 


The logical Pauli operators are subject to the transformation rule of the CNOT gate. Ac¬ 
cordingly, the transversal CNOT operation A(X) acts as a logical CNOT operation for the 
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encoded degree of freedom. Because the Hadamard, phase, and CNOT operations are im¬ 
plemented transversally, whole Clifford group elements can be implemented by transversal 
operations. 

The transversal implementation is fault-tolerant because the operations apparently do 
not increase the number of errors on a code block and there is no internal interaction 
between the qubits in the same code block. Combined with a fault-tolerant gadget for 
measuring the error syndrome as explained in Appendix I A. 11 we can implement Clifford 
operations fault-tolerantly. 

Unfortunately, the non-Clifford operation does not transform a Pauli operator into 
another Pauli operator. For example, the vr/S operation transforms the Pauli X 

operator into a Clifford operator: 

^-i{-Kl%)z = (X + y) / V2. (2.80) 

This implies that a transversal non-Clifford operation hardly results in a logical non- 
Clifford operation. Thus, a fault-tolerant non-Clifford gate operation is not so straightfor¬ 
ward. To settle this, we can utilize magic state distillation, consisting of noisy non-Clifford 
resource states and ideal (fault-tolerant) Clifford operations, as explained in the next sec¬ 
tion. 

2.8 Magic state distillation 

2.8.1 Knill-Laflamme-Zurek protocol 

A fault-tolerant implementation of a non-Clifford gate was first proposed in an earlier 
paper by Knill, Laflamme, and Zurek |KLZ98bl IKLZ98a] . Instead of implementing the 
non-Clifford gate directly, we consider a fault-tolerant preparation of the non-stabilizer 
state, the so-called magic state, 

|A) = (2.81) 

The magic state can be utilized to implement a non-Clifford gate A = by using 

one-bit teleportation consisting of Clifford gates and the Pauli basis state preparations and 
measurements. Thus, if we can prepare a clean magic state, we can create a fault-tolerant 
non-Clifford gate by using fault-tolerant Clifford gates. 

The Knill-Laflamme-Zurek construction of the fault-tolerant preparation of the magic 
state was based on the fact that |A) is an eigenstate of H. The Hadamard operation has 
transversality, e.g., on the 7-qubit code. Hence, if we perform a projective measurement 
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of we obtain a clean magic state. The circuit is given as follow 
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(2.82) 


where A = and we used the fact that = Ac^t{H). The above 

circuit consists of three parts, encoding of the logical magic state into the 7-qubit code, 
indirect measurement of and decoding by one-bit teleportation. Note that all Clifford 
operations are assumed to be ideal, because they are easily made fault-tolerant by using 
a stabilizer code. With an appropriate randomization operation made up by Clifford 
operations, a noisy magic state can be transformed into a mixture of |A) and lA.-*-) = Y\A): 


PA = il-p)\A){A\ +p\A-^){A\ 


(2.83) 


Hence, the noise on the magic state can be expressed as a T error. The Y error can 
be detected by transversal Z measurements for the decoding (see also Knill’s gadget for 
a fault-tolerant syndrome measurement in Appendix lA.ip . Assuming an ideal Clifford 
operation, the error probability p decreases as 0{p^), when we employ the 7-qubit code 
with distance 3. Including the A gate for the initial encoding, we need 15 noisy magic 
states I A) to obtain a clean magic state. This distillation protocol works for all self-dual 
CSS codes, which are symmetric between the X- and Z-type stabilizer generators and has 
transversality for the Hadamard gate. Recently, improved protocols have been proposed 
based on this approach |MEK12( l,Tonl3b| . 


2.8.2 Bravyi-Kitaev protocol 

Bravyi and Kitaev proposed another magic state distillation protocol based on a 15-qubit 
code [BK05] . While their and the Knill-Laffamme-Zurek protocols seem quite different, 
they are, interestingly, known to be equivalent [Rei05] . In the Bravyi-Kitaev protocol, a 
[[15,1,3]] quantum code is defined by the [[15,7,3]] classical Reed-Muller code, whose parity 
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check matrix is given by 
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so that HxHJ = 0 and the Z-type stabilizer generators are defined similarly. The logical 
operators are given by Lx = and Lz = The logical states are written down 

explicitly as 


|0l) = n(^ +4^)100-0)’ 

i=l 
4 

lu) = 

i=l 

The number of Is in each term of |0 l) and |1 l) is 8 and 7, respectively. By applying the 
T = gate transversally, we obtain 

( 2 . 88 ) 

= e-*"/®|0L). (2.89) 

Thus, the transversal T gate acts as a logical gate. Note that this transversality 
does not hold in the orthogonal (erroneous) subspace, e.g., spanned by {X^IOl),X^ ll^,)}. 
However, we can show that this is enough to perform a fault-tolerant logical T gate. 
Instead of applying the T gate directly, we implement it using a one-bit teleportation: 


( 2 . 86 ) 

(2.87) 



\T) 



(2.90) 


56 





The Z-basis measurement and CNOT operation are both implemented transversally on 
a CSS code. Thus, if the preparation of the non-Clifford ancilla state |T) = (IOl) -|- 
e-*^/4| 1 1 ))/ \/2, called a magic state, is done fault-tolerantly, one can ensure fault-tolerance 
of the logical T gate. 

By using an appropriate randomization process, we can prepare a noisy magic state as 
follows: 


Pt = {1-p)\T){T\+pZ\T){T\Z. (2.91) 

Thus, a phase error Z is located on the ideal magic state with probability p. This phase 
error causes a Z error after the T gate by one-bit teleportation. Because the code space 
is invariant under the transversal T gate, we can detect such a Z error by the following 
circuit. 



The first part of the above circuit consisting of CNOT operations is an encoding circuit 
for the quantum Reed-Muller code. The transversal T gate is applied by using one-bit 
teleportation. The logical |r) state is measured in the X-basis transversally, which detects 
Z errors on the code state, projecting the code state on the local X-basis. The input state 
in the second lowest wire is entangled with the ancilla qubit in the lowest wire, where the 
distilled magic state is teleported. 

Let c be a 15-bit string specifying the location of the Z errors, E{c) = Yli ■ 
If E(c) commutes with the X-type stabilizer generators, the state passes through the 
distillation circuit. To calculate this probability, we define a weight enumerator of a 
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subspace V G GF{2^) 


Wv{x,y) = Y, (2.92) 

c&V 

The probability of passing the distillation circuit is calculated to be 

Ppass = Wy^ (1 -p,p) = (1,1 - 2p) = l + 15(l-2p)^ ^ ^2.93) 

"a: \Vx\ iO 

where the orthogonal subspace is equivalent to the kernel of Hx, Ker(//x)- We also 
used the MacWilliams identity |MS77] : 

Wv{x,y) = ^^Wv^{x + y,x-y). (2.94) 

Similarly, the error probability of the output can be calculated to be 

Wvh, (p, 1 - p) = 


Accordingly, the error probability, under the condition of passing the distillation circuit, 
is given by 

If p' > 0.141, we can reduce the error probability on the magic state via the distillation 
circuit. After I rounds of distillation, the error probability decreases to {\/^p)^‘/^/^. At 
each round, we need 15 noisy magic states. Because the probability of successfully passing 
the distillation circuit converges rapidly to 1, the average number of noisy magic states 
consumed after I rounds becomes 15^ Accordingly, the average number of noisy magic 
states required to achieve an error probability e of the magic state scales like 

[log(V^e)/log(V^p)]'°§(^®)/'°§® =o(iog2-5g). (2.98) 

In Sec. a we saw that if an input state is a convex mixture of the Pauli basis 
states followed by Clifford operations and Pauli basis measurements, the measurement 
outcomes can be simulated classically in the weak sense. The noisy magic state px = 
(1 — p)|r)(r| + p|r)(r| lies on a line in the x-y plane, as shown in Fig. 12.51 If p > 
(I —\/2/2)/2 = 0.146, lies inside the octahedron, and the Gottesman-Knill theorem is 
applicable. On the other hand, if p < 0.141, magic state distillation allows us to implement 
universal quantum computation with an arbitrary accuracy as seen above. Unfortunately, 
magic state distillation based on the Reed-Muller code does not provide a tight distillation 
threshold against the classically simulatable region. In Ref. [Rei05] . a distillation protocol 
using the 7-qubit code was proposed and achieved a tight threshold p = (1 — \/2/2)/2. In 
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this sense, the classically simulatable and quantum universal regions are divided sharply 
on the x-y plane. 

By combining the magic state distillation and fault-tolerant Clifford operations on 
the CSS code, we can perform universal quantum computation fault-tolerantly. In order 
to make the error probability arbitrarily small, we can employ concatenated quantum 
computation, in which logical qubits of a lower concatenation level are utilized as the 
physical qubits at a higher level. At the higher level, all operations, including logical qubit 
preparations and syndrome measurements for QEC, have to be done fault-tolerantly. If 
the error probability is smaller than a certain constant value, which we call the noise 
threshold, the logical error probability at the highest concatenation level decreases super- 
exponentially. On the other hand, the overhead increases exponentially. Thus, we can 
make the logical error probability small enough to maintain a quantum computation of 
size N with a polylogarithmic overhead polylog(A^). This implies that we can obtain a 
quantum benefit even for a quantum algorithm with a quadratic speedup, such as the 
Grover algorithm |Gro96| . In Appendix lA.ll we briefly review fault-tolerant syndrome 
measurements, concatenated quantum computation, and the threshold theorem. 

While the resource increment for protecting quantum computation scales polylogarith- 
mically, its constant factor is quite huge. Almost all overheads for the fault-tolerant quan¬ 
tum computation are employed for the magic state distillation [VMLFYlOt I.IVMF~^T^ . 
Thus, much effort has been spent recently on developing resource-efficient magic state 
distillation |MFK12l IBH121 IFasl3l IJonl3a] . 
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measurements 



logical operators multiplied by 
stabilizer generators 


Figure 2.4; (a) A graph state and a measurement pattern, (b) The logical X operator of 
input 1 is multiplied by the stabilizer generators and we obtain a correlated operator on 
outputs 1 and 2 (left). The logical Z operator of input 1 is multiplied by the stabilizer 
generators and we obtain the logical Z operator of output 1. The gray colored X operators 
are replaced by ±1 depending on the measurement outcomes. 



Figure 2.5: The noisy magic state in the Bloch sphere. 
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Chapter 3 


Topological stabilizer codes 


Protecting quantum information from decoherence is of prime importance to realize quan¬ 
tum information processing. Several approaches have been proposed toward reliable 
quantum information processing, ranging from passive to active protections, such as 
decoherence-free subspaces |LCW98^ , dynamical decoupling [VKL99b] , and quantum er¬ 
ror correction |Sho95] . Among these, the most comprehensive approach is fault-tolerant 
quantum computation based on quantum error correction |DS96] . Quantum fault-tolerance 
theory ensures scalable quantum computation with noisy quantum devices as long as the 
error probability of such devices is smaller than a threshold value (see Appendix lA.ip . 
The first threshold values were obtained, at almost the same time (1996), independently 
by Aharonov et al ., who used a polynomial code of distance five [AB0971 IABO08| ; Knill 
et al., who used the Steane seven-qubit code [KLZ98bl IKLZ98a] : and Kitaev, who used 
toric (surface) codes [Kit97] . All of these studies were based on a concatenated quan¬ 
tum computation and achieved similar threshold values ~ 10“®. Since then, fault-tolerant 
quantum computing has been studied as one of the most important issues in quantum 
information science. In 2005, Knill proposed a novel scheme based on the C^/Cq error¬ 
detecting code, namely the Fibonacci scheme [Kni05] . and achieved a considerably higher 
threshold, a few %. Recently, Fuji! et al. have further improved this approach by making 
use of measurement-based quantum computation (MBQC) [RHG061 IRHG07j on logical 
cluster states, which gives the highest threshold value so far, ~ 5% [FYlOallFYIOb] . 

All these approaches rely on the availability of two-qubit gates between arbitrarily 
separated qubits. However, all physically available interactions are local in space. More¬ 
over, if we consider manufacturing convenience and the selective addressability of in¬ 
dividual qubits, 2D nearest-neighbor architectures are necessary. It was thought that, 
if we restrict the two-qubit gates to nearest-neighbor ones in 2D, then the threshold 
value would decrease significantly |STD07| . The situation changed completely by the 
proposal made by Raussendorf et al., i.e., topologically protected quantum computa¬ 
tion |RHn7[ iRHGOTl IRHGOB] . They utilized a surface code, a kind of topological quan¬ 
tum code, originally proposed by Kitaev as a toric code on a torus |Kit03] and inves¬ 
tigated by Dennis et al. in detail |DKLP02] . All the stabilizer generators of the sur- 
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face code are spatially local, as will be seen later, and hence syndrome measurements 
are done by using nearest-neighbor two-qubit gates. More importantly, Raussendorf 
et al. developed a novel way to implement a fault-tolerant logical gate operation by 
braiding defects on the surface, which can be done by using only single-qubit measure¬ 
ments and nearest-neighbor two-qubit gates. Nevertheless, its threshold value is very high 
~ 1% [RH07[ [RHG071 rRHGOGi IWFHllj . Based on these analyses, physical implementa¬ 
tions and quantum architecture designs have recently been suggested |VMLFYT0llLBSB10l 
IFTml l■TVMF+12[ IFYKT1 21 ITTbT^ INT.B1.3l l(TF(T12[ IFMM(112) . which clarifies the experi¬ 


mental requirements for building a fault-tolerant quantum computer. On the other side, 
extensive experimental resources have been used to achieve these requirements, and very 
encouraging results have already been obtained in several experiments in superconducting 


systems [OOC+12[ ICOM+14[ IBKM+Uai ICMS+14| . 


Another important motivation for studying topological quantum codes are their con¬ 
nection with topologically ordered systems in condensed matter physics [Kit03] . Topo¬ 
logical stabilizer codes, whose stabilizer generators are local and translation invariant, 
provide toy models for topologically ordered systems. Using quantum coding theory with 
a geometrical constraint such as locality and translation invariance, we can understand 
the nature of topologically ordered condensed matter systems. One of the most important 
issues in this direction is to answer whether or not thermally stable topological order exists 
in three or lower dimensions [BTOQl lYosllb] . 

When we consider the decoding problems of topological quantum codes, we encounter 
classical statistical physics |DKLP02] . Specifically, the posterior probabilities for the error 
correction correspond to partition functions of classical spin glass models |DKLP0^ . such 
as the random-bond Ising model. If the corresponding random-bond Ising model lies inside 
a ferromagnetic phase, the error correction on a surface code succeeds, and the logical error 
probability decreases exponentially. 

In this chapter, we introduce a representative example of the topological stabilizer 
codes, the surface codes. Before that, we also briefly review the Z 2 chain complex, which 
is a useful tool to define the surface codes. We also apply it to a classical repetition 
code as an exercise. After introducing the surface codes, we explain how topological error 
correction is performed. We also mention the relation between topological error correction 
and the classical spin glass model. Another example of topological stabilizer codes, the 
topological color codes, is also presented. Finally, we explain the connection between 
topological stabilizer codes and topologically ordered systems studied in condensed matter 
physics. 


3.1 Z 2 chain complex 

Before defining the surface codes, we introduce a mathematical tool, the Z 2 chain complex, 
which is useful for describing the surface codes. 

Consider a surface G = {V, E, F) that consists of the vertices V = {vk}, edges E = {e;}, 
and faces F = {fm}- (Here we consider a graph G = (U, E) embedded on a 2D manifold. 
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the surface, on which the faces F are defined.) We define a vector space Cq over Z 2 , using 
each vertex Vk G V as a basis B{Cq) = {ufc}. A vector in Cq 

Co = '^ZkVk, (3.1) 

k 

where Zk = {0,1}, is referred to as a 0-chain. The vector space Cq is an Abelian group 
under component-wise addition modulo 2. Similarly, we can define the Abelian groups Ci 
and C 2 using the edges E = {c;} and faces F = {fm} as base vectors B{Ci) = {c;} and 
B{C 2 ) = {fm}, respectively. The elements 

Cl = '^ziei, (3.2) 

i 

C2 = '^Zmfm, (3.3) 

m 

where zi,Zm = {0,1}, are called the 1-chain and 2-chain, respectively. With a mild abuse 
of the notation, a set of i-dimensional elements (i.e., vertices, edges, and faces) is also 
specified by a i-chain Cj G C* as a set of elements having Zj = \. 

We can define a homomorphism di : Ci ^ Ci-i such that 

di o di-i = 0. (3.4) 

Specifically, diCi is defined as an (i — l)-chain that is the boundary of c*. For example, 
dei is a 1-chain for which = 1 if and only if (iff) Vk is an endpoint of c;. Thus, the 
homomorphism di is called a boundary operator. When there is no risk for confusion, we 
will denote di simply by d. A chain c* is called a cycle, if it is in a kernel of the boundary 
operator di, i.e., dci = 0. For example, the boundary d 2 fm of a face fm is a cycle, because 
^2 0 01 = 0. Such a sequence of Abelian groups Ci, connected by the boundary operators di 
with diodi-i = 0, is called a Z 2 chain complex. We define a trivial i-cycle c* if there exists 
an (i -|- l)-chain such that Cj = dci+i, i.e., Cj G Img(5i+i). By regarding the trivial cycles 
as generators of the continuous deformation, we may define an equivalent class having the 
same topology. More precisely, a homology group hi is defined as a quotient group formed 
by the cycles and the trivial cycles: 

hi = ker(5j)/Img(9j+i). (3.5) 

An element of the quotient group hi is called a homology class. If two i-chains q and 
c' belong to the same homology class, there exists an (i -|- l)-chain Cj+i such that c* = 
c' -|- dci+i- Two such i-chains Cj and c' are said to be homologically equivalent. 

We also define the dual surface G = {V,E,F), where we identify the elements of the 
dual and original (primal) lattices such that V = F, E = E and F = V. Specifically, the 
dual lattice is constructed such that the two vertices v,v' G V are connected by an edge 
e, if the corresponding two faces / and /' share the same edge e. We can also define a Z 2 
chain complex on G using the dual bases B{Ci), a dual i-chain q G Ci and a boundary 
operator di : Gi ^ Q-i. When there is no risk of confusion, di will be denoted by d. 
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In the construction of the surface codes, a qubit is defined on each edge ei £ E of the 
surface G (or equivalently each dual edge ei of the dual surface G). The Pauli product is 
defined using a 1-chain ci = ziei (or dual 1-chain ci) such that 

lT(ci) = n<S (3.6) 

i 

where Wi G {X;, Y), Z;} is a Pauli operator acting on the qubit on edge e;. Specifically, we 
have 


W{ci)Wic[) = W{ci + c'l). (3.7) 

Consider the two operators X{ci) and Z{c'i), defined by the two 1-chains ci and 
respectively. The commutability of these two operators is determined by the inner product 
of the two chains (vectors) ci • c'^ = Yhi where the addition is taken modulo 2: 


Cl • = 0, iff X{c\)Z{c!i) = Z{c!i)X{ci) (commute), (3.8) 

Cl • Cl = 1, iff X{ci)Z{ci) =—Z{(^i)X{ci) (anticommute). (3.9) 

Let M{di) be a matrix representation of di with respect to the basis vectors B{Gi) and 
B{Gi-i). We have 

{M{di)ci) ■ Ci-i = Ci ■ {M{difci-i). (3.10) 

By identifying the primal and dual bases B{Cq) = B{G2), B{Ci) = B{Gi), and B{C2) = 
B{Go), the duality relation between primal and dual lattices can be expressed by 

M{di) = M{d 2 f, M{d 2 ) = (3.11) 

where indicates the matrix transpose. 


Using Eq. (13.4p di o 82 = 0, we obtain M{d 2 )'^M{d 2 ) = 0. Thus, for any primal and 
dual 2-chains, C 2 and C 2 , we have 

dc 2 ■ dc 2 = 0, (3.12) 

i.e., X{dc 2 ) and Z{dc 2 ) always commute. This property is useful to construct a stabilizer 
group, because stabilizer generators have to commute with each other. Using Eqs. (|3.10p 
and (13.lip , we also have 

Cl • d 2 C 2 = dci ■ C2. (3.13) 


3.2 A bit-flip code: exercise 

As an exercise, we define a bit-flip code (classical repetition code) using the Z 2 chain 
complex. (Readers who are familiar with this topic may skip this section.) Consider 
a regular polygon G{V,E,F) on a sphere consisting of n = |E| edges and two faces 
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Figure 3.1: (a) A regular polygon embedded on a sphere, (b) A stabilizer operator defined 
as Ak = Z{dfk) = 0(51)^ The logical X operator is given by Lx = X{dfi). Because 
primal and dual 1-chains share an even number of qubits, Lx and Ak commute with each 
other. 

F = {/i,/ 2 } (corresponding to the top and bottom hemispheres) as shown in Fig. 13.11 
The number of qubits is n. 

We define a stabilizer generator for each dual face fk = as follows: 

Ak = Z{dfk)= n (3T4) 

eiGdvf, 

where, for convenience, 6 vk = df is defined as a set of edges that are incident to the vertex 
Vk- (Note that both the dual and primal objects are identihed.) 

Because A^ = /, there are n — 1 independent stabilizer generators. The di¬ 

mension of the stabilizer subspace is two. The code subspace is described by a one-cycle 
Cl = dfi = 9/2 surrounding the sphere. We define the logical X operator 

Lx =X{dfi) = X{ci). (3.15) 

Note that Lx = X{dfi) and A^ = Z{dfk) commute, due to Eq. (13.121) . The logical Pauli 
Z operator, which commutes with {A^} and anticommutes with Lx, is dehned as Lz = Zi. 
(We may choose any of the Z;s on the edge ei, because they act equivalently on the code 
space.) The code is a bit-flip code, which protects quantum information against bit flip 
errors. 

A string of bit errors X{ci) is defined using a 1-chain ci, which we call an error chain. 
The error is detected by measuring the eigenvalues of the stabilizer generators A^ = 
Z{6vk), i.e., through syndrome measurements. From Eq. (13.131) . we have ci -dvk = dci -Vk- 
Thus, X{ci) anticommutes with the stabilizer generator Ak if dci ■ Vk = 1. Hence, X(ci) 
anticommutes with Ak on the vertex Vk that is a boundary Vk G dci of the error chain 
Cl. The eigenvalue of the stabilizer generator becomes -1. Error correction is the task of 
finding a recovery chain c[ such that d{ci -|- c^) = 0. In the case of the bit-flip code, there 
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error syndrome: dci 


Figure 3.2: An X error chain X{ci) and the corresponding error syndrome dci. The 
eigenvalue of the stabilizer generator becomes —1 on vertices in dci. 

are only two possibilities ci + = 0 or ci + = dfi. In the former case, error correction 

succeeds, while the latter case results in a logical error Lx- 

The bit-flip code cannot protect quantum information against phase errors. A natural 
extension of the classical repetition code to handle both bit and phase errors is the surface 
codes introduced in the next section. 


3.3 Definition of surface codes 

3.3.1 Surface code on a torus: toric code 

Now we are ready to define the surface code. Let us consider a square lattice G = {V, E, F) 
on a torus with a periodic boundary condition as shown in Fig. 13.31 A dual square lattice 
G = (V, E, F) is also defined on the torus. We define stabilizer generators of the Z- and 
A-types for each face fm and vertex Vk as follows: 

Am = Z{dfm), Bk = X{5vk)=X{dfk). (3.16) 

Because dfm ■ dfk = 0 from Eq. (j3.12l) . Am and Bj. commute. The stabilizer generators 
Am and B^ are called plaquette and star operators , respectively. By the definition of the 
stabilizer code, the code state I'L) satisfies 

Arnm = I'L), = m, (3.17) 

for all fm & F and Vk E V. 

If two 1-chains, ci and d^, are homologically equivalent, the actions of Z{ci) and Z{di) 
on the code state are the same, because there exists a 2-chain C 2 such that 

Z{d,) = Z{ci)Z{dc 2 ) = Z{ci) I n I • (3-18) 

\/mGC2 / 
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Figure 3.3: A square lattice on a torus with a periodic boundary condition. The Z and 
X stabilizer generators are defined by being associated with each face and vertex (dual 
face). The logical Z and X operators are defined by non-trivial cycles. 


We will denote this simply by Z{d^ ~ Z{ci). This is also the case for homologically 
equivalent dual 1-chains, ci and c'l, and the actions of X{ci) and A(c'^) on the code state. 
Hence, the homology classes, which are equivalent classes over trivial cycles, of the primal 
and dual 1-chains correspond to the actions of Z- and A-type operators on the code state, 
because their actions are equal up to stabilizer operators. 

Let us define the logical operators Z{ci) and A(ci). The logical operators have to com¬ 
mute with all stabilizer generators and also be independent of them. The former condition 
implies that dci = 0 and dci = 0 from Eq. (13.131) . This is because the commutability im¬ 
poses that Cl • 5vk = dci ■ Vk = 0 for all vertices k and ci • dfm = dci ■ Vm = 0 for all 
faces (dual vertices) m. The latter condition implies that ci and ci are nontrivial cycles, 
because for two homologically equivalent cycles ci and c[, the actions of Z{ci) and Z{c[) 
on the code state are the same as seen in Eg. (13.181) . We can find two non-trivial cycles, 
which belong to different homology classes, for each primal and dual 1-chain, as shown in 
Fig. 13.31 (This is a natural consequence, because the homology group hi on the torus is 
Z2 X Z2.) Then, we define two pairs of logical Pauli operators 


{4') = Z(cS')),4) = ^(44 and {4^ = Z{4 


,( 2 ) 


(i) 

Note that and anticommute with each other if i = 


),4^=A(44. (3.19) 

j. Otherwise, they commute: 


LfL^ = (-l)'5-4)4). 


(3.20) 


That is, they satisfy a commutation relation equivalent to that of the Pauli operators for 
two qubits. The logical Pauli basis states are defined as follows: 

44z(s1,S2)) = (-1)^'|'&z(5i,S2)), 

44x(s1,S2)) = {-ir\^x{si,S2)). 
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(3.21) 

(3.22) 






The number of stabilizer generators of the n x n square lattice on the torus is given by 


|F| + \F\-2 = |F| + \V\-2 = 2r? - 2, (3.23) 

where —2 comes from the fact that O/^eF -^rn = I and OufeeF = I, and hence there 
are two non-independent operators. On the other hand, the number of qubits is given by 

\E\ = 2r?. (3.24) 

Thus, we have a = 2^-dimensional stabilizer subspace. The above two 

pairs of logical operators appropriately describe the degrees of freedom in the code space. 
The code distance, the minimum weight of the logical operators, is determined by the 
linear length of the square lattice. Thus, the surface code on the torus is an [[re^,2,n]] 
stabilizer code. 

In short, an operator that commutes with the stabilizer generators corresponds to a 
cycle, i.e., ker(5i). The stabilizer operators correspond to the boundaries of the 2-chains 
C2 G 6 * 2 . The logical operators, which commute with and are independent of the stabilizer 
generators, correspond to the homology class hi = ker(cli)/Img(cl 2 ). We summarize the 
correspondence between the surface code and the Z 2 chain complex in Tab. 13.11 

Table 3.1: The correspondence between Z 2 chain complex and stabilizer codes. 


stabilizer code 

chain complex 

Z-type stabilizer generator 

faC6 fm 

A-type stabilizer generator 

vertex Vk (dual face fk) 

Z- and A-type stabilizer operators 

boundaries of 2 -chains, 

Img( 52 ) and Img( 92 ) 

commutability 

cycle conditions: dci = 0 and dci = 0 

the operators commuting with stabilizer generators 

ker(5i) and ker(9i) 

the logical Z and A operators 

nontrivial cycles ci and ci 
(homology classes [ci] = ker(( 9 i)/Img( 52 ) 
and [ci] = ker( 9 i)/Img(( 92 )) 

Z and A errors 

1 -chains ci and ci 

Z and A error syndromes 

dci and dci 


All properties so far hold for a general tilling G = {V,E,F), and hence we can define 
a surface code on general, e.g., triangular and hexagonal, lattices [FT12] . Moreover, the 
numbers of edges, faces, and vertices are subject to the Euler characteristic formula 

\F\ + \V\-\E\=2-2g, (3.25) 

where g is the genus, the number of “handles” of the surface. Thus, the dimension of 
the stabilizer subspace is calculated to be = 2^^. These are the degrees of 

freedom equivalent to 2 g qubits in the stabilizer subspace. 
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Figure 3.4: A planar surface code. Top and bottom are smooth boundaries consisting of 
three-qubit star operators. Left and right are rough boundaries consisting of three-qubit 
plaquette operators. At the boundary, we can define sets of 1-chains Fi = { 71 } G Ci 
and Fi = { 71 } G Ci for the primal and dual lattices, respectively, by which the relative 
homology is defined. 


Because almost all arguments are made from the operator viewpoint (Heisenberg pic¬ 
ture), there is no need to write down S 2 )) explicitly. However, for the readers who 

worry about such things, we will give an explicit description from the state viewpoint; 


|^z(si,S2)) 


|^x(si,S2)) 




I + Bk 




Afx{c^^^y^x{cf^y^ I II 


I + Ar, 


| 0 )®n 2 

(3.26) 

1 1 +)®"' 

(3.27) 


L /m€F 


where M is a normalization factor. The code state can be viewed as an equal weight 
superposition of all cycles belonging to the same homology class. 


3.3.2 Planar surface code 

The periodic boundary condition might be hard to implement in experiments. We can 
also define a surface code on a planar n x (re — 1 ) square lattice with an appropriate 
boundary condition as shown in Fig. 13.41 The top and bottom boundaries consist of 
three-qubit star operators, which we call smooth boundaries, because they are complete 
plaquette operators. On the other hand, the left and right boundaries consist of three- 
qubit plaquette operators, which we call rough boundaries. The X operators on a dual 
1-chain can terminate at the smooth boundary, while the Z operators on a 1-chain can 
terminate at the rough boundary. Thus, we should use a relative homology to define logical 
operators, where two chains Cj and c( are said to be (relative) homologically equivalent iff 
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Figure 3.5: A planar surface code with a defect. The stabilizer generators are not defined 
inside a defect, which provides a nontrivial cycle wrapping around it. 

= Cj+9ci+i+7i with 7 j G Fj C Cj. In this case, Fi is chosen, specifically, to be the vector 
space spanned by the set of 1 -chains each of which consists of three edges corresponding to 
the three-qubit plaquette operator at the top and bottom smooth boundaries. Similarly, 
we also define Fi at the left and right rough boundaries of the dual lattice. Because 
Z("fi) and A( 71 ) are both stabilizer operators, if the shapes of two logical operators are 
homologically equivalent, their actions on the code state are the same. The number of 
edges, faces, and vertices are now \E\ = In? — 2n + l, |F| = n? — n, and |F| = — n. Thus, 

the stabilizer subspace is a 2D subspace. We can define the logical operators Lz = Z{ci) 
and Lx = A(ci) using horizontal and vertical 1-chains ci and ci, as shown in Fig. 13.41 
If one chooses all boundaries smooth, the stabilizer state is uniquely defined, and 
hence there is no logical degree of freedom. By punching a hole, we can introduce a 
defect on the surface, whereby we can define a nontrivial cycle for a logical Z operator, as 
shown in Fig. 13.51 The logical X operator can be chosen as a dual 1-chain that connects 
the defect and the smooth boundary, as shown in Fig. 13.51 As mentioned previously, 
the properties of the logical operators are the same if the corresponding 1 -chains are 
homologically equivalent. 

3.4 Topological quantum error correction 

Let us return to the surface code on the torus to explain how errors are corrected. There 
are orthogonal subspaces, the so-called syndrome subspaces, each of which is an 

eigenspace of the stabilizer generators and has the same structure as the code space. These 
orthogonal subspaces are utilized to identify the location of errors and to infer a recovery 
operation. Suppose the X and Z errors A(cf) and .Z’(cf), defined by error 1-chains cf 
and cf, respectively, occur on the surface code, as shown in Fig. 13.61 The code state is 
now mapped into one of the orthogonal subspaces. From Eq. (13.1311 . A(c|) and .Z’(cf) 
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Figure 3.6: (a) A Z error chain Z(ci) is detected at the boundary dci. The eigenvalue of 
the stabilizer generator on the face /fc £ C 2 = 5ci becomes -1. (b) An X error chain 
X(ci) is detected at the boundary dci. The eigenvalue of the stabilizer generator Am on 
the face /m £ C 2 = dci becomes -1. 

anticommute with Am and on the face fm £ 9cf and vertex £ dcf. Thus, the 
orthogonal subspace is specified by error syndromes 

= c|, del = Cq. (3.28) 

More precisely, the eigenvalues with respect to the stabilizer generators Am and B^ are 
given by (-1)^- and (-l)^fc, where c| = and = Yjk^k'^k (see Fig. [3T]), 

respectively. Error correction is the task of finding recovery 1-chains cl and cl such that 
9(cf + c^) = 0 and 9(cf + cl) = 0 , meaning that the state is returned into the code space 
by applying X(cl) and Z(cl), respectively. Below we will, for simplicity, explain only how 
to correct the Z errors, but the extension to the X errors is straightforward. 

Suppose each Z error occurs with an independent and identical probability p. Con¬ 
ditioned on the error syndrome Cq = def, the posterior probability of an error Z{ci) 
occurring with ci = zici is written as 

PMc‘o) = vn(^)"b„.,. (3.29) 

where AA is a normalization factor. One way to find a recovery chain is by maximizing the 
posterior probability 

cl = argm^axP(ci|cg) = argmin l 9 ci=cg, (3.30) 
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Figure 3.7: From left to right, error chain Z(cf), recovery chain Z(ci} with = Cq, and 
net actions Z(ci + cf). If + cf is a trivial cycle (top), error correction succeeds. If + cf 
is a non-trivial cycle (bottom), the error correction results in a logical error. 


where ci = 'Y^iZiei. As seen in the l.h.s., this corresponds to minimizing the number of 
errors such that dc^ = Cq. Hence, this is called minimum distance decoding. If 

Cl + cf is a trivial cycle, the error correction succeeds, because the net action Z(ci + cf) 
after the recovery operation is identity on the code state as shown in Fig l3.7[ If + cf is 
a nontrivial cycle, the recovery operation results in the logical operation Z(ci + cf) ~ Lz, 
i.e., a logical error. 

Minimum distance decoding is hard in general, because it can be mapped into an inte¬ 
ger programing problem, which is NP-hard |HLGll[|KL12l[iP13j . However, in the present 
case, there is a nice geometrical property that makes the decoding problem feasible. The 
condition dci = Cq and min^ read that it is sufficient to find a 1-chain that connects 
pairs of two vertices in Cq with a minimum Manhattan length. There is a classical poly¬ 
nomial time algorithm to do such a task, the so-called minimum-weight perfect matching 
(MWPM) algorithm of Edmonds |Edm65l lBar82[ [CR99] . The algorithm scales like 0(n®), 
with n being the linear length of the lattice and with the fixed error probability p. Typ¬ 
ical examples of the error chain cf and the error plus recovery chain cf -|- cf are shown 
in Fig. 13.81 for each p = 0.05 (b), p = 0.10 (c), and p = 0.15 (d). Here, we employ an 
implementation of MWPM, blossom V [Kol09| . The higher the physical error probability 
is, the longer the error plus recovery chain becomes. For a high physical error probability 
p = 0.15, the error plus recovery chain becomes a large cycle, and unfortunately results 
in a logical error. Such a logical error probability is plotted as a function of the physical 
error probability p in Fig. 13.81 (a) for each n = 10 (solid line), n = 20 (dashed line), and 
n = 30 (dotted line). If the error probability is sufficiently smaller than a threshold value, 
the logical error probability decreases for increasing n. The threshold value for decoding 
by the MWPM algorithm has been estimated to be 10.3% (MWPM) [DKLPO^ IWHPO^ . 

The minimum distance decoding with MWPM is not optimal for our purpose, i.e.. 
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(c) p=10% 
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(b) p=5% 


(d) p=15% 



Figure 3.8: (a) The logical error probability is plotted as a function of the physical error 
probability p. Error chains cf and error plus recovery chains cf + cf for p = 0.05 (b), 
p = 0.10 (c), and p = 0.15 (d). 

making the logical error probability as small as possible. An error correction with another 
recovery chain , for which c\ + Ci belongs to the same homology class as cf+ C]|, provides 
exactly the same result. We may use such a recovery chain Ci to correct the error. Thus, 
we should maximize, not each posterior probability p{ci\cq), but a summation of it over 
the same homology class. This problem originated from the degeneracy of the surface 
code, where each syndrome is assigned not uniquely, but for many error instances. A 
prototypical example of an error syndrome, for which we should consider not olny the 
weight of errors, but also combinatorics (an entropic effect) of the error configurations, is 
shown in Fig. 13.91 

Denoting the homology class by hi, the posterior probability for a homology class hi 
is given by 

B= F(c;VS). (3.31) 

cl'\ci+cl'ehi 

where c[ is a recovery chain satisfying dc[ = Cg and chosen arbitrarily as a reference frame, 
and the summation is taken over all 1-chains such that -|- c[ belongs to the homology 
class hi. 

The posterior probability may be rephrased by using the stabilizer language as follows 
(see also Appendix [B|). Let Q and £ be the stabilizer and logical operator groups, respec¬ 
tively. For a given error syndrome Cg, we define the recovery operator Z(ci) a priori, such 
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Figure 3.9: Suppose the top and bottom boundaries are connected by a periodic boundary 
condition. The minimum-weight error is shown in the left panel. The error shown in the 
right panel is not minimum-weight but has a fourfold degeneracy. The recovery chain 
should be chosen by comparing the total error probabilities and 4 x p^. 


that the erroneous state is returned into the code space. We can decompose an arbitrary 
error operator Z(ci), providing the syndrome Cq, into 

Z{cl) = Z{c\)GLi, (3.32) 


where G £ Q and Li £ C are a stabilizer and logical operator, respectively. The posterior 
probability of the logical operator Li is calculated by summing over all stabilizer operators 

G£g 


Pi = P{Li\cl) 


1 j;P[Z(c5)GL,]. 
^ Gag 


(3.33) 


We choose the most likely homology class hi or, equivalently, the most likely logical oper¬ 
ator Li that maximizes the probability pi = P{Li): 


Lj = argmaxP(Lj|co). (3.34) 

Li 

The error correction is completed by applying LiZ{c\), and the logical error probability 
is given by 1 — pi- 

In the next section, we relate the posterior probability summed over the same homology 
class to a partition function of the random-bond Ising model. The partition function of the 
Ising model on a planar graph with general coupling strengths (without magnetic fields) 
can be calculated in polynomial time using the Kasteleyn-Barahona algorithm with the 
Pfaffian method [Kasbll rBar82( IBROT] . Thus, the optimal decoding is also implemented 
by a polynomial time classical processing, though it takes more overhead than MWPM. 
The threshold value for optimal decoding is also discussed in the next section. 

Efficient decoding has been one of the most important issues for the realization of fault- 
tolerant quantum computing, because the coherence time of quantum information would 
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stabilizer operators and pure errors 



edge operators 


(b) 6 

6 6 6 



Figure 3.10: (a) A unit cell of the level-1 logical qubit. Two pairs of logical Pauli operators 
(top). The 6 stabilizer operators and the pure error operators, each of which anticommutes 
with an stabilizer operator (middle). The 4 pairs of edge operators. These 12 pairs of 
mutually anti-commuting operators generate the Pauli group of the 12 qubits, (b) Using 
the level-1 logical operators, the level-2 stabilizer generators are defined, (c) Using 12 
level-1 logical qubits, a unit cell of the level-2 logical qubit is defined. 

be very short; a fast classical processing is essential. Fowler et al. proposed an efficient 
decoding method based on MWPM |FWH12] . Because a long error chain is exponentially 
suppressed, we can assign weights between each pairs of vertices having —1 eigenvalues 
according to their Manhattan length. This allows us to reduce the number of edges from 
0{n^) to O(n^). Moreover, the exponential suppression of longer error chains allows us to 
search the pairs within a local small region; matching a long-distance pair is exponentially 
rare. Because the matching process employs almost exclusively local information, this 
algorithm can be parallelized to an 0(1) average time per round. 

Another decoding method is using a renormalization technique [DCPIO^IDCPIOI^ . As 
explained in Appendix[Bl we can efficiently execute an optimal decoding on a concatenated 
quantum code by using a brief-propagation on a tree factor graph |Pou06| . Although the 
surface code itself does not have such a hierarchal structure, a renormalization technique 
is employed to introduce a hierarchal structure on the surface code [DCPlOa] IDCPlOb] . 

Using 12 qubits as a unit cell, we define a couple of level-1 logical qubits, as shown 
in Fig. 13.101 (a), which include 2 pairs of logical operators 6 stabilizer generators 

We define 6 pure error operators each of which anticommutes with a stabilizer 
generator as shown in Fig. 13.101 fa) (middle). Moreover, we define 4 pairs of anticommuting 
operators as shown in Fig. 13.101 (a) (bottom), which we call edge operators to generate 
the Pauli group of the 12 qubits. Any Pauli operator A on the unit cell can be decomposed 
in terms of these operator, 

A = (3.35) 
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where E for B = L,G, G, E. The pure error operator is chosen uniquely 
according to the error syndrome to return the state into the code space. (If a stabilizer 
generator Gi has an eigenvalue —1, then we employ Gi, which anticommutes with Gi.) 
From the error distribution P{A), the posterior probability of the level-1 logical operator 
is calculated by taking a marginal over and 

P{A = (3.36) 

which are utilized to model the error distribution at the level 2. 

Similarly, a level-fe unit cell is defined by using 8 level-(fe — 1) unit cells and the 12 pairs 
of the level-(fe — 1) logical operators on them. Similarly to the level-1 case, we define the 
level-/c logical, stabilizer, pure error, edge operators. At the highest level k = I, we obtain 
the logical operators of the surface code. For example, the level-2 stabilizer generators are 
shown in Fig. 13.101 (b). A unit cell of the level-2 logical qubit is shown in Fig. 13.101 (c). 

The posterior probability of the level-A; logical operator is calculated us¬ 

ing the posterior probabilities at the level {k — 1) by assuming that they are independent 
for each level-(fc — 1) unit cell. Under this assumption, we can calculate the posterior 
probability conditioned on all error syndrome by using the belief propa¬ 

gation [PouOOl IDCPTTR IDCPI 0b| (see Appendix [B] for decoding the concatenated codes 
by using the belief propagation). Then the maximization procedure over all logical oper¬ 
ators provide us a most likely logical operator, 

L* = argmaxP(L*^*^|5^^^). (3.37) 

L(0 

The level-(A; — 1) unit cells share physical qubits with each other, and hence the conditional 
logical error probability is not independent for level-(fc — 1) unit cells. If we employ a 
message passing in order to reweigh the level-(/c — 1) logical error probability as the error 
model for the level-A: unit cell, we can further improve the approximation of P(L^^'^\S^^'^). 
If the decoding process is done in parallel, the belief propagation takes 0(log2(n)) time, 
similar to the case of the concatenated code as explained in Appendix [Bj Because the 
posterior probability is approximated, the decoding based on the renormalization is not 
optimal, but results in a reasonable threshold value, ~ 9% for the independent X and 
Z error and ~ 15.2% for the depolarizing error (against 10.3% and 15.5% obtained by 
MWPM, respectively) |DCP10^ IDCPlOb] . The renormalization decoder can also be 
applied for an arbitrary topological code [BDCPl^ . 

3.5 Error correction and spin glass model 

The behavior of the logical error probability in Fig. 13.81 suggests the existence of a critical 
phenomenon behind the error correction problem. Indeed, there is a beautiful correspon¬ 
dence between quantum error correction on the surface code and a spin glass model, the 
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so-called random-bond Ising model (RBIM) [DKLP02]. More precisely, the posterior prob¬ 
ability Eq. (j3.31l) of the logical operator is mapped into a partition function of the RBIM 
as seen below. 

To solve the condition dci = Cq in Eq. (j3.29j) . we rewrite the recovery chain as ci = 
c^'“-|-c^, where c\ € Img(92) is atrivial cycle, is further decomposed into = cf-|-c^^\ 
where cf determines the actual location of errors, and is a (nontrivial) cycle belonging 
to the homology class hk- Note that this decomposition corresponds to Eq. ()3.32p . The 
posterior probability is rewritten as 




(3.38) 


where c“ = zf^ei, with z" E {0,1}. In order to take the condition c\ E Img(52) auto¬ 
matically, we introduce a gauge valuable Zm £ {0,1} on each dual vertex Vm and a dual 
0-chain Cg = Ylm ^mVm- Any trivial cycle c\ is replaced by the gauge valuables using the 
following relation 


4 = 0 4 

Vm&dSl 


(3.39) 


For example, if we choose Cq = hm, we obtain c\ = dfm- (This corresponds to a multi¬ 
plication of the face stabilizer generator defined on the face fm-) There is a one-to-one 
correspondence between a trivial cycle and a gauge dual 0-chain. Using Cq we can formally 
solve the condition c\ E Img(92) in Eq. p.29p . 


^’(cilcS) 


vn 


I 



'vm^dei ' 


(3.40) 


The binary valuables Zi E {0,1} are transformed into spin variables Uj E {-|-1,—1} by 
Uj = (—1)^\ Moreover, we define a coupling constant e~‘^ = y^p /(1 — p). Then Eq. ()3.40p 
is rewritten as 


Pici\4) 


M'e 


JT.I 


I m{l) m'{l) 


(3.41) 


where Af is a normalization factor and m(Z) and m'{l) are end points of the edge e;. 
By changing the notation I —>• ij and i,j, the posterior probability is 

reformulated as a Boltzmann factor of the ± J RBIM: 

P(ci|cg) = (3.42) 

where J-4 = JcrR- edges where the errors are located, anti-ferromagnetic in¬ 

teractions are assigned. The posterior probability of the logical operator Eq. (I3.3ip is 
calculated by taking summation over all gauge spin configurations (this corresponds to 
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the summation over all stabilizer operators in Eq. (j.S..Sljl l: 


E ^(LlcS) 

(3.43) 

ci'\c['-+ci'ehk 


E p(c;*+c‘iic5) 

(3.44) 

c\\c\elmg{d 2 ) 



(3.45) 

wn 



(3.46) 


where is the partition function of the ±J RBIM. 

Let us consider the performance under this decoding by taking an average of the 
logarithm of pk over the error distribution, which corresponds to a sample average with 
respect to quenched randomness: 


P'k = MPk)] 


'^P{c\)Pk 

^1 

{•^ij } 

-Ffe + InAA', 


(3.47) 

(3.48) 

(3.49) 


where [•] indicates the sample average and is the sample average of the logarithm of 
the logical error probability p^. The quenched randomness is determined by the error 
distribution 


F(4) = {l-p')5{Jt^ - 1) +p'5(4 + 1). (3.50) 

Note that p is a parameter in the posterior probability and could be different from the 
actual error probability p'. Of course, optimal decoding, in the sense of Bayesian inference, 
is achieved by maximizing the posterior probability with a true error probability p = p' ■ 
Fq = —[lnZ({J®})] is the free-energy of the RBIM, while F^ = {k = 

1,2,3) is the free-energy with respect to the interactions where a domain wall 

(k) 

of anti-ferromagnetic interactions corresponding to a cycle c) of a homology class is 
inserted. 

From Eqs. (j3.46jl and (I3.49h . the relation between quantum error correction and a 
spin glass model becomes apparent; the posterior probability of the logical operator is 
proportional to the partition function of RBIM, whose Hamiltonian is given by F = 
— Ylij location of the Z errors, represented by = — 1, corresponds to 

the anti-ferromagnetic interaction due to disorder as shown in Fig. 13.111 (see also Tab. l3.2|l . 
The error probability and the coupling strength (inverse temperature) are related by e~'^ = 
y/p /(1 — p). The error syndrome Cq = dcf corresponds to the distribution of frustrations 
of the Ising interactions and also the end points (Ising vortex) of the excited domain walls. 
The probability of the homology class is expressed by the domain-wall free energy: 
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Figure 3.11: A gauge spin aj located at the face-center of a plaquette. The distribution 
of anti-ferromagnetic interactions corresponds to the error chain cf. The error syndrome 
Cq = 9cf corresponds to the distribution of frustrations. The ground state configuration 
is determined by a domain-wall consisting of excited domain-wall and anti-ferromagnetic 
bonds, both of which have frustrations at their end points. 


-{P^^-p^^) = Fk-Fo. (3.51) 

If the physical error probability p is smaller than the threshold value, —{k = 1,2,3) 
diverges in the large n limit. On the other hand, if p is higher than the threshold value, 
~P^k ~ converges to 2 In 2 meaning the stored information becomes completely 

destroyed. This non-analytical behavior also appeared in the r.h.s. of Eq. (I3.49p . i.e., the 
free-energy of RBIM. Indeed, the difference between the free energies A = F^ — Fq is 
an order parameter, the so-called domain-wall free energy |Huk99| . for a ferromagnetic 
ordered phase. The phase diagram of RBIM with respect to J = (l/2)ln[p/(l — p)]) and 
p' is shown in Fig. 13.131 As mentioned before, optimal decoding is achieved with p = p', 
i.e., e-'^ = vWO" — p'), which is called the Nishimori line |Nis81| . The critical point on 
the Nishimori line, which is called a multi-critical point, has been numerically calculated 
to be 10.94 ±0.02% by Honecker et al. and Merz et al. |RDQDS^ IHPPOlt IMC02] . The 
optimal threshold value of the surface code is good in the following sense: The existence 
of CSS codes with asymptotic rate R = k/n is guaranteed if 

R=l- 2H2 {p), (3.52) 

by the quantum Gilbert-Varshamov bound under independent X and Z errors with prob¬ 
ability p [CS96| . The rate becomes zero with p = 11.00%. The optimal threshold of 
the surface code, which consists only of local stabilizer generators, achieves a value very 
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AF = excited domain wall 




Figure 3.12: Examples of ground states in a ferromagnetic phase are shown for 
with k = 0,1. In the ferromagnetic case, the insertion of the anti-ferromagnetic bonds 
results in an excited domain-wall and hence increases the free energy (ground state energy 
at zero temperature). 


close to this. Indeed, Nishimori conjectured that the multi-critical point of the RBIM is 
determined by 


H2{p) = 1/2, (3.53) 

arguing from the self-duality in RBIM with a replica method |NNn2| INis79l IKW41] (if 
the reader is interested this derivation, please see a review in Ref. [Ohzl3| ). Ohzeki has 
evaluated the multi-critical point precisely to be 10.92% by using a real-space renormal¬ 
ization [Ohz09b| . which is in a good agreement with the numerical result (HTPVOS]. 

The minimum distance decoding with MWPM is achieved in the limit p —)• 0, which is 
the low temperature limit J —?■ oo where the entropic effect is suppressed. The threshold 
of MWPM corresponds to the critical point at zero temperature [Bar82j . which has been 
investigated numerically and determined to be 10.4 ±0.1 by Kawashima et al. [KR971 
IWHP031IDKLP02] . 


3.6 Other topological codes 


The surface codes have been also studied on general lattice tillings, such as, triangle, 
hexagonal, and random lattices |Ohz09bl IRWH+1^ IFT121IOF121 [XSWP131 [Kayl4| . Sup¬ 
pose a surface code is defined on a lattice G = {V,E,F). As seen in the previous section, 
the X error correction is mapped into a RBIM on a lattice G = {V,E). On the other 
hand, the Z error correction is mapped into a RBIM on its dual lattice G = {y,E). The 
mutual duality relation [NNr)2] between G and G allows us to predict the relation between 


80 































































probability of AF bond 

Figure 3.13: A phase diagram of RBIM with respect to the coupling strength J = 
(l/2)ln[p/(l — p)] and the probability of anti-ferromagnetic interaction p'. The proba¬ 
bility p in j corresponds to the parameter in the posterior probability. When p = p', 
called the Nishimori line, an optimal decoding is achieved. 


the optimal threshold values for X and Z errors: 


H{p^) + H{p:,) = 1 , 


(3.54) 


where px and pz are the X and Z error probabilities, respectively. This equality is the same 
as the quantum Gilbert-Varshamov bound evaluated for the independent X and Z errors 
with probabilities px and pz, respectively. The precise locations of the optimal thresholds 
for regular lattices have been investigated by Ohzeki using a real-space renormalization 
technique [Ohz09b| . The thresholds have been investigated using the MWPM algorithm 
by Fujii et al. [FT12| . The thresholds with MWPM (i.e., the critical points of RBIMs at 
zero temperature) even approaches Eq. (I3.54jl . as shown in Fig. 13.141 These codes with 
an asymmetry between the X and Z error tolerances would be useful to correct a biased 
error |ABD~*~09 . In Refs RWH'*~12 , IOF12| , the asymmetry is continuously controlled by 
changing a lattice parameter. 

A leakage process or qubit loss is an important source of noise. Unlike the X and Z 
errors discussed so far, the qubit loss is detectable (heralded), and hence we can tolerate 
much more loss rate than for the (undetectable) error rate. Stace et al. proposed to cope 
with the qubit loss on the surface code [SBDn9l ISBIO] . Suppose a qubit is lost on the 
surface as shown in Fig. 13.151 (a). The plaquette operators containing the lost qubit are 
undetermined. However, we can construct a super-plaquette multiplying two neighboring 
plaquette so that the super-plaquette does not have the lost qubit. By using the super- 
plaquette as a stabilizer generator, we can perform MWPM. If two super-plaquettes are 
neighbors, the weight between these super-plaquettes has to be modihed appropriately 
in MWPM, because the plaquettes share two physical qubits and the error probability is 
effectively doubled. The logical operator is defined by avoiding the lost qubits. Unless 


81 




























Table 3.2: The correspondence among random-bond Ising model, Z 2 chain complex and 
stabilizer codes. 


RBIM 

stabilizer code 

chain complex 

dual lattice 

Z error correction 

primal lattice 

primal lattice 

X error correction 

dual lattice 

Ising interactions 

qubits 

edge ei 

gauge spin 

Z-type stabilizer generators 

faC6 fm 

gauge spin 

A-type stabilizer generators 

vertex Vk (dual face fk) 

domain wall 

logical X and Z operators 

nontrivial cycles ci and ci 

anti-ferromagnetic interactions 

X and Z errors 

1-chains ci and ci 

frustration 

X and Z error syndromes 

dci and dci 


the lost qubits, which are located on edges, are percolated throughout the lattice, we can 
find such a logical operator. Thus, the threshold for qubit loss in the large lattice limit 
without any (undetected) error is determined by the bond percolation threshold. The 
trade-off curves between qubit loss and (unheralded) error rates with MWPM for various 
lattice tillings are shown in Fig. 13.141 (b). 

On the other hand, when the error correction problem is mapped into the RBIM, the 
qubit loss corresponds to a bond-dilution. The duality relation in the presence of the 
bond-dilution on the square lattice is given by [Qhzl2| 

{I-q)h{p)+q= 1/2. (3.55) 

From the numerical data in Refs. [FT12] . we could expect a more general equality, 

(1 - qx)h{px) + (1 - qz)h{pz) + qx + qz = 'i-, (3.56) 

where qx and qz are the probabilities of the heralded X and Z errors, respectively. With 
Px = Pz = 0, this is reduced to qx + qz = ^, which corresponds to Kesten’s duality relation 
of the bond percolation thresholds between mutually dual lattices. With qx = qz = 0, Eq. 
(j3.56jl is reduced to Eq. (I3.54jl . With qx = qz = q and Pz = Px = P, Eq. (I3.56p is reduced 
to Eq. (|3.55l) . 

Another important class of local stabilizer codes is the topological color code proposed 
by Bombin and Martin-Delgado [BMPObl IBMDOT] . The topological color codes are defined 
on trivalent graphs with faces that can be colored with three colors, such as hexagonal 
(6.6.6) and (4.8.8) lattices, as shown in Fig. 13.171 (a) and (b), respectively. A qubit is 
defined on each vertex v and stabilizer generators are defined on each face /: 

=n =n 

vef vef 

Specifically, the color code on a (4.8.8) lattice, shown in Fig. 13.161 (c), allows all single¬ 
qubit Clifford gates transversally. The distance-3 topological color code on the (4.8.8) 


82 






























p\OSS 


Figure 3.14: (a) The thresholds {px,Pz) of the surface codes on square, Kagome, hexagonal, 
and triangle-hexagonal lattices, as well as their duals. The curve is the quantum Gilbert- 
Varshamov bound hitting zero asymptotic rate under a bit and phase flip channel with 
probabilities px and pz, respectively, (b) The trade-off curves between qubit loss rate 
pioss unheralded error rates px and pz are shown for square, Kagome, hexagonal, and 
triangle-hexagonal lattices. If the loss and error rates are below these curves, the logical 
information is protected. 


lattice is equivalent to Steane’s 7-qubit code |Ste96j . The extension to a 3D lattice also 
enables a transversal non-Clifford gate on the code space |BMD07] . whose distance-3 
version corresponds to the Reed-Mullar 15-qubit code |BK05| . 

The topological color codes are also described by a Z 2 chain complex on hyper¬ 
graphs |Dell4] . Consider a trivalent graph G = {V,E), on which a topological color 
code is defined. We define a hyper-graph G = {V,£,T) consisting of the sets of hyper¬ 
vertices V, hyper-edges £, and hyper-faces E as follows (see Fig. 13.17]) : The hyper-vertices 
V are vertices V of the dual graph G. A hyper-edge e G £ is defined as a triplet of 
(hyper-)vertices v G V on a face / of the dual graph G. A hyper-face f G T is defined as 
a set of hyper-edges e that are incident to the vertex v of the dual graph G. As defined 
in Sec. EH we can define both Z 2 valued vector spaces and Abelian groups C' 1 ^ 2,3 with 
the hyper-graph elements as their bases. Because the hyper-edge and hyper-face are de¬ 
fined as sets of hyper-vertices and hyper-edges, respectively, we can define boundary maps 
di : Gi ^ Gi-i naturally by such sets. The plaquette and star stabilizer generators of the 
surface code defined on such a hyper-graph are as follows: 

Am=X{ Z,, Bfc = n (3.58) 

eea/ egiSC 

This definition is equivalent to the previous definition, Eq. (I3.57h . In Ref. |Dell4| . the 
authors dehned the projections from the Z 2 chain complex on the hyper-graph into chain 
complexes on a dual graph G. Specifically, a certain subset of vertices is removed in each 
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Figure 3.15: (a) A super-plaquette defined as a product of two plaquette operators. The 
lost qubit is not contained in the super-plaquette operator, (b) The logical operator 
is chosen appropriately by avoiding the lost qubits. The lost qubits (bonds) are not 
percolated, so we can find such a logical operator. 

projection. This procedure corresponds to the removal of stabilizer generators colored by 
one of the three colors. This allows us to utilize MWPM to decode the topological color 
codes efficiently on the projected surface codes [DelUj . 

The decoding problem of topological color codes is mapped into random three-body 
Ising models using the mapping between quantum error correction and spin glass models. 
This can be understood as follows. For each face center a gauge spin valuable is located 
associated with each stabilizer generator; three gauge spin valuables (stabilizer genera¬ 
tors) that share the same qubit (vertex) interact with each other. The locations of the 
optimal thresholds have also been investigated via the spin glass theory and Monte Carlo 
simulations [OhzOQat . 

3.7 Connection to topological order in condensed matter 
physics 

Topological order is an exotic quantum phase of matter, which cannot be characterized 
by the Landau-Ginzburg theory of symmetry breaking, in which an ordered phase can 
be characterized by local order parameters. The ground state of a topologically ordered 
system is degenerate, but its degeneracy cannot be destroyed by any local perturbation. 
Thus, no local order parameter can succeed to capture the topologically ordered phase. 
Understanding the nature of topological order is one of the most important goals of modern 
condensed matter physics. Moreover, the ground state of a topologically ordered system is, 
by definition, robust against any local perturbations, and hence it is also useful for storing 
quantum information. There is a beautiful correspondence between topological quantum 
codes and topologically ordered systems, which provides a promising way to understand 
condensed matter physics via quantum information. 

Let us first consider the bit-flip code defined in Sec. 13.21 The Hamiltonian, which we 
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Figure 3.16: (a) and (b) Trivalent lattices that can be colored with three distinct colors, 
(c) Color codes defined on (4.8.8) lattices with a triangle open boundary condition. The 
logical operators are shown by solid lines. Any string at the boundaries can be employed 
as a logical operator. 

call a stabilizer Hamiltonian, is defined as follows: 

n—1 

.falsing = —J = — J (3.59) 

k i=l 

Here, the summation is taken over all stabilizer generators, but one stabilizer operator, 
which is not independent, is removed. The Hamiltonian Eq. (j3.59p corresponds to the 
Ising model in ID with an open boundary condition. By its construction, the Hamiltonian 
is diagonalizable, and the stabilizer state becomes the ground state. The bit-flip code has 
a 2D stabilizer subspace spanned by {|00...0), 111...!)}. This means that the ground state 
is degenerate. The bit-flip errors occurring on the code space excite the ground state to an 
excited state. Thus, the states in the orthogonal subspaces of the bit-flip code correspond 
to excited states. 

To address topological order, let us consider the robustness of the ground state against 
perturbations from transversal fields hxYli^i- By using standard perturbation the¬ 
ory [Kat95] . we can show that the degeneracy of the ground states is not lifted up to 
the {n — l)th order of perturbation. This happens because the code distance of the n- 
qubit bit flip code against X errors is n, and any X errors of weight up to n — 1 map the 
code state into an orthogonal space. Accordingly, the ground state degeneracy is robust 
against the transversal fields. 
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Figure 3.17: (left) A hexagonal lattice with three-coloring, (right) A triangular lattice, the 
dual of the hexagonal lattice. A hyper-graph is defined on a triangular lattice. The color 
code on the hexagonal lattice corresponds to a surface code defined on the hyper-graph. 

On the other hand, longitudinal fields spoil the ground state degeneracy in 

the large n limit, even if hz is small. More precisely, the energy between |00..0) and 
|ll..l) is shifted by nhz- Thus, a superposition Q!|00..0) -|-/3|11..1) in the ground subspace 
is easily destroyed by the longitudinal fields. In this sense, the stabilizer Hamiltonian 
constructed by the bit-flip code is not topologically ordered. However, if any perturbation 
with respect to the longitudinal fields is prohibited due to some symmetry of nature, the 
ground state degeneracy is robust under that symmetry. This type of robustness of ground- 
state degeneracies is called symmetry protected topological order [GWOQirPTBOlO] . (Note 
that in this case the ground-state degeneracy is not related to the geometrical property of 
the underlying manifold in contrast to the genuine topological order in 2D.) 

The symmetry prohibiting the longitudinal fields seems to be somewhat artificial. We 
can, however, impose the symmetry by transforming the Ising Hamiltonian to a free- 
fermion model using the following Jordan-Wigner transformation |JW28] : 


i—l 

021-1 = XkZi, 

(3.60) 

k=l 


i—l 


^2i 

k=l 

(3.61) 


The operators, called Majorana fermion operators, are hermitian = a|, and satisfy 
the fermion commutation relation { 0 ^, 0 ^/} = The Ising stabilizer Hamiltonian is 

reformulated in terms of a^: 

n—1 

-^^Isng J ^ • (3.62) 

i=l 
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Figure 3.18: (left top) (1) A pair of A-type excitations are created, moved around the 
torus, and annihilated. (2) A pair of A-type excitations are created, moved around the 
torus, and annihilated. (3) Do the process (1) again. (4) Do the process (2) again, 
(middle) The creation, movement, and annihilation process is continuously deformed, 
(right bottom) A braiding operation of an A-type excitation around a Z-type excitation. 


The logical operators acting on the degenerated ground states are given by 

/n-l \ 

Lz = ai = Zi, Ly = aia2n = ( IT 


(3.63) 


\k=2 


The degree of freedom in the degenerated ground states is called the Major ana zero mode 
or the unpaired Major ana fermion [KitOl] . If the parity of the number of fermions is 
preserved, the fermion operators would appear with a quadratic form a^a^. Under such 
a symmetry, there is no perturbation that lifts the ground state degeneracy. Thus, the 
ground state of the unpaired Majorana fermion is symmetry protected. 

Next, we will provide a genuine topologically ordered system based on the surface 
code (Kitaev’s toric code). The stabilizer Hamiltonian, the so-called Kitaev’s toric code 
Hamiltonian |Kit97j . is given as a summation of all plaquette and star operators: 


TI^Kitaev — J ^ ^ J ^ ^ . 


(3.64) 


The ground state has a fourfold degeneracy corresponding to the code space. Errors on 
the code state correspond to excitations. Specifically, there are two types of excitations, 
corresponding to the Z error Z(ci) and A error A(ci). Excitations appear at the bound¬ 
aries of the error chains dci and dc\, because the local energy changed from —J to +J 
there. Such excitations, i.e., at the end points of the error chains, are always created as 
pairs, can be viewed as a pair creation process on the ground state. 

Suppose these two types of the excitations were created on the system and moved as 
shown in Eig. 13.181 (left top). This process can be described by 

Z(cf))A(cf))Z(cS'))A(cW)|T) = -|T). (3.65) 
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On the other hand, by continuously changing the trajectory of the particles, as shown in 
Fig. 13.181 this process can also be regarded as a braiding process of an X-type excitation 
around a Z-type excitation. After the braiding operations, a phase factor is applied 
to the state as shown in the r.h.s. of Eq. (j3.65j) . Thus, the excitations are neither 
bosonic nor fermionic, which are invariant under the braiding operation, i.e., the swapping 
operation twice. In this sense, the excitations on the surface code are referred to as anyons. 
Specifically, since the phase factors is Z 2 , they are called Z 2 Abelian anyons. By using 
the generalized Pauli operators on for the qudit, we can also define Kitaev’s toric 
code |Kit03| . on which excitations are Z^ Abelian anyons. More generally, using a finite 
group G, we can define a quantum state \g) {g E G) in a |G|-dimensional Hilbert space. 
Then we define four types of operators for each g ^ G: 

L\ = Y. \gh){hlL^_ = ^ \hg-^){h\,Tl = \h){hlT^ = \h-^){h-\ (3.66) 

h£G heG 

The non-Abelian Kitaev’s toric code model, which is called the quantum double model [KitOdj . 
is defined as 


H = -jY,Mfm)-jY,Bivk), (3.67) 

m k 

in terms of the following plaquette and star operators 

MM = (3.68) 

91929394=1 

BM) = (3.69) 

' ' 96G 

where the four edges E dfm and £ ^Vk = dfk are labeled clock wise. The 

quantum double model supports non-Abelian anyonic excitations |Kit03l IPacl2] . which 
allows us to implement universal quantum computation solely by braiding them. 

Let us return to the Z 2 Kitaev’s toric code Hamiltonian. The code distance of the 
surface code on the n x n torus is n. Thus, neither X-type (transverse) nor Z-type 
(longitudinal) fields can lift the ground state degeneracy up to the (n — l)th order of 
perturbation. More generally, no local perturbation can lift the ground state degeneracy 
in the large n limit. Thus, the ground state of the Kitaev’s toric code Hamiltonian is 
topologically ordered. For any properly defined stabilizer code, we can define a stabilizer 
Hamiltonian, whose ground state is topologically ordered. However, in condensed matter 
physics, the local interactions are of central importance. Thus, it is natural to restrict the 
stabilizer generators to be spatially local, i.e., topological stabilizer codes. 

Only the nearest-neighbor two-body interactions are attained in physically natural 
systems. It has been known that the Kitaev’s toric code Hamiltonian is obtained as an 
effective low energy model of another model consisting only of two-body nearest-neighbor 
interactions |Kit06j . Let us consider the following two-body nearest-neighbor model, called 


















Figure 3.19: A hexagonal lattice on which Kitaev’s compass model is defined. In the large 
Jz limit, the two spins on each vertical edge are confined into a 2D subspace forming a 
dimmer. The dimmer is located on each edge of a square lattice shown by the dotted lines. 


Kitaev’s compass model: 


Hcomp =-Jx ^ XiXj - Jy ^ YiYj — Jz ^ ZiZj, (3.70) 

{i,j)^Ey ii,j)eEz 


where Ex, Ey, and E^ are sets of right-up, left-up, and vertical bonds, respectively, of a 
hexagonal lattice (see Fig. 13.191 (left)). If we take the large Jz limit, each vertical bond 
favors the two-dimensional subspace spanned by {|00),|11)}, because it is stabilized by 
ZiZj. Thus, in the large Jz limit, we can derive an effective low energy Hamiltonian, which 
commutes with the ZiZj interactions, by using perturbation theory: 


J2j2 

Tj _ ^ y 

-ti pff — 


16|J^ 


Vx fX fYfY 

z e{ e§ e 


f 

3 ^4 


Jz Y. 

e£Ez 


(3.71) 


where Ag = AiAj with e = (i,j) and A = X,Y,Z, depending on e G Ex, Ey, Ez- The 
edges {e{} are left-up and right-up edges on a hexagonal face / and the summation Ylf 
is taken over all faces. 

Let us define a qubit {|0) = |00), |1) = 111)} for each dimerized edge and Pauli opera¬ 
tors X = XX OT = YY and Z = Z ® I or = I ® Z. Now a qubit is assigned to each edge 
of a square lattice, which is an vertical edge in the hexagonal lattice as shown in Fig. 13.191 
(right). Using this definition, the effective Hamiltonian can be reformulated as 


^efr = - 


16|J.|3 


/ 


iZ^jZjZj, 


(3.72) 


f f f f 

where {ej, ei, ei, ey} are the left, right, top, and bottom edges, respectively, on a square 
face /. If we apply the Hadamard transformation on all horizontal edges, we obtain the 
Kitaev’s toric code Hamiltonian. In this way, the stabilizer Hamiltonian can be obtained 
































as a low-energy effective model of a two-body nearest-neighbor system [KitOB] . This 
shows the validity of employing topological stabilizer codes and quantum coding theory 
to understand the quantum phase of matter in condensed matter physics. 

A complete classification of the topological stabilizer codes in 2D has been obtained 
by Yoshida |Yoslla| (see also the classification by Bombin et al. |BDCP12j ). Specifically, 
the quantum phases of the 2D stabilizer Hamiltonians are classified by the geometric 
shapes of the logical operators. The thermal stability of the topological order in stabilizer 
Hamiltonian systems at finite temperatures has also been studied via quantum coding 
theory by Bravyi and Terhal |BT09] for 2D, and Yoshida for 3D |Yosllb| . 

If there is a thermally stable topological order, we can store quantum information re¬ 
liably even at a finite temperature without any active error correction, i.e., we have a self- 
correcting quantum memory. Of course, if a fault-tolerant quantum computer were real¬ 
ized, we could store quantum information reliably with a repetitively performing quantum 
error correction, which, however, requires selective addressing of each individual qubits. 
There are also several intermediate approaches for a reliable quantum storage without 
selective addressing using global dissipative dynamics [FNIK14] . an interaction with an 
engineered environment |HCC09l IPHWLl^ IKCS14] , and decoding by cellular automata 
with local update rules [HCEKl^ IHar04| . However, a genuine topologically ordered self- 
correcting quantum memory seems to be hard to achieve in 2D even in the presence of 
effective long-range interactions [LCYPPTK] . 
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Chapter 4 


Topological quantum computation 


In this chapter, we explain how to perform topological fault-tolerant quantum computation 
on the surface code. All operations employed are nearest-neighbor (at most) two-qubit 
gates and signle-qubit measurements on a 2D array of qubits. This property is quite favor¬ 
able for the fabrication of qubits and their control lines on a chip. Fault-tolerance ensures 
that no local noise during any sort of operations ever spoil the quantum computation, if 
the noise strength is smaller than a certain threshold value, as seen below. 

We first introduce the defect pair qubit, a logical qubit using a pair of defects on 
the surface, which allows us to encode many logical qubits on the surface. Then, we 
explain the elementary operations of defects by local quantum information processing. A 
braiding operation, based on the elementary operations, of the defects is further utilized 
to implement the logical CNOT gate between the defect pair qubits. Combining with the 
topologically protected operations, the magic states are injected and distilled for fault- 
tolerant universal quantum computation. Based on these understandings, we introduce a 
topological diagram and the topological calculus on it, as a set of the transformations that 
preserve the logical actions on the code space. These provide us with a deep understanding 
of topological quantum computation on the surface. Finally, we return to the microscopic 
viewpoint to explain how topological quantum error correction is executed and analyzed. 

4.1 Defect pair qubits 

In order to perform quantum computation of N qubits, we have to arrange N logical qubits 
in the code subspace while keeping their code distance long enough. A way of doing this 
is to use a surface of higher genus. Specifically, we use a large planar surface and punch 
holes on it, which we call defects on the surface. 

In Sec. 13.3.21 we have seen that we can define a logical qubit by introducing a defect on 
a planar surface, where the stabilizer generators inside the defect region are removed from 
the stabilizer group. The logical operators are defined by a nontrivial cycle surrounding 
the defect and a chain connecting the defect to the boundary. If we introduced many 
degrees of freedom, i.e., many defects based on this strategy, the logical operators defined 
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Figure 4.1: Primal (left) and dual (right) defect pair qubits. The primal defect pair qubit 
is defined by removing the plaquette operators inside the defect regions D and D'. Then 
the cycle dD around the defect becomes a nontrivial cycle, by which the logical Z operator 
Z{dD) is defined. The logical X operator X{ci) is defined as a dual 1-chain ci connecting 
two defects. Any homologically equivalent logical operators acts the same way on the code 
space. The dual defect pair is defined similarly, but the basis is changed by the Hadamard 
transformation. 

by chains connecting each defect to the boundary would become complicated. To avoid 
this, we use a pair of defects to define a logical qubit, which we call a defect pair qubit. 

We define two defect regions D £ C 2 and D' & C 2 , as shown in Fig. 14.11 We remove all 
Z-type stabilizer generators inside these regions, i.e., {Am}/^g_DuD'- Thus, the operators 
Z{dD) and Z{dD') are not stabilizer operators. Instead, we append a Z-type stabilizer 
operator Z{dD + dD) as a stabilizer generator. Moreover, we also append the Pauli X 
operator on all edge qubits inside (not including the boundary) the regions D U D', i.e., 

By this definition, the star (A-type) stabilizer generators in the 
defect region are still in the stabilizer group. We choose Z{dD) as a logical operator, 
because it commutes with all stabilizer operators, but does not belong to the stabilizer 
group. Because Z{dD + dD) is a stabilizer operator, Z{dD') also acts the same way, i.e., 
Z{dD) ~ Z{dD'). Moreover, the actions of operators represented by the homologically 
equivalent cycles are the same. The logical X operator is given by A(ci), with a dual 
1-chain ci connecting two defects, as shown in Fig. 14.11 The actions of the operators 
represented by any homologically equivalent chains in the sense of relative homology are 
the same. That is, we may choose any dual 1-chain ci which connects two defect regions 
D and D'. The code distance is given by the circumference of the defect or the distance 
between two defects. Because the plaquette operators on the primal lattice are removed, 
we call these defects and the logical qubit primal defects and primal defect pair qubit, 
respectively. 
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Figure 4.2: (a) A defect creation. The qubits inside the defect regions (excluding the 
qubits on the boundary), denoted by circles, are measured in the A-basis. (b) The defect 
region D is expanded into D' by measuring the qubits inside the region D' in the A-basis, 
where the encoded quantum information is preserved, (c) A defect region D' is contracted 
into D by measuring the plaquette operators in D" = D'\D. The logical operator Z after 
the contraction is defined depending on the measurement outcome m. 


Similarly, we can also define a logical qubit by removing the star (A-type) stabilizer 
operators on the dual defect regions D and D' defined on the dual lattice as shown in 
Fig. 14.21 We call such defects and the logical qubit dual defects and dual defect pair qubit, 
respectively. Hereafter, the planar surface code state, on which the defects are introduced, 
is referred to as vacuum. Below we will explain how defect qubits are created, deformed, 
and moved in the vacuum. 

4.2 Defect creation, annihilation, and movement 

The defect creation is accomplished by measuring the qubits inside the defect region D, not 
including the qubits on its boundary, in the A-basis [see Fig. l4.2l (al]. These measurements 
remove the plaquette operators inside D from the stabilizer group, because these A- 
basis measurements do not commute with them. On the other hand, the measurements 
do commute with the star operators, and hence a parity of four measurement outcomes 
corresponding to a star operator is always even (if there is no measurement error). While 
the measurement outcomes are random, we can prepare all measured qubits to be in the 
1+) state by applying Z operators according to the measurement outcomes. (In practice. 
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there is no need to apply the Z operations. It is enough to keep the information of the 
byproduct Pauli operator dependent on the measurement outcomes.) 

As the qubits on the boundary are not measured, the post-measurement state is sta¬ 
bilized by Z{dD). In order to create a pair of defects, we do the same thing for another 
defect D'. Apparently, the resultant state is stabilized by Z{dD + dD') and satisfies the 
definition of the defect pair qubit. Moreover, the logical qubit is stabilized by Z{dD) and 
hence a logical Z-basis state is prepared. 

The defect region can be expanded by creating a larger defect D\ which includes the 
original defect region D, i.e., D C D', as shown in Fig. 14.21 (b). We can choose Z[dD') 
as a logical Z operator of the defect pair qubit, so that the information with respect to 
Z{dD') is preserved during this operation. (Recall that logical operators act the same if 
the corresponding cycles are homologically equivalent. Thus we can choose a large enough 
cycle surrounding the defect in advance.) We can expand the defect region step-by-step 
such that, at each step, a logical X operator is untouched as shown in Fig. 14.21 (b). (Or 
equivalently, we can compensate for the difference in the logical X operator according 
the A-basis measurement outcomes inside the defect region.) Thus, the information with 
respect to the logical X operator is also preserved. Accordingly, the logical information 
of the defect pair qubit is stored, but now the defect region D is expanded into D + D'. 

The defect annihilation is executed by measuring the plaquette operators inside the 
region D to restore them to the stabilizer group. The surface with a defect D is rewritten 
by 



where \v) indicates the surface code state without defect, i.e., vacuum. Thus, \D) is a 
superposition of all possible applications of X operators on the vacuum |u). The mea¬ 
surements of the plaquette operators collapse the superposition. By applying a recovery 
operation inside the region D such that all eigenvalues becomes -|-1, the defect is annihi¬ 
lated. (Note that there is no need to actually apply the recovery operation. It is enough 
to keep the record of the eigenvalues.) The parity of all measurement outcomes of the 
plaquette operators inside D corresponds to the eigenvalue of Z{dD), because we have 

z{dD) = n z{dfm) = n ( 4 . 2 ) 

fm^D fmSD 

This indicates that a logical Z-basis measurement can be done by annihilating the defect 
completely. 

Suppose a defect region D inside the defect region D' (i.e., D G D') is annihilated by 
measuring the qubits inside D, as shown in Fig. 14.21 (c). We can then obtain the eigenvalue 
(_l)"i 2) of Z{dD). Let D” be the complement of D in D'. Because we have 

Z{dD') = Z{dD")Z{dD), (4.3) 
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Figure 4.3: A defect movement by expansion and contraction. 

depending on the eigenvalue (—1)™', the operator {—l)'^Z{dD") acts as the logical operator 
of the defect D". Thus, the defect D', consisting of the defect pair qubit, is contracted 
into D" without changing the stored quantum information (up to the logical Pauli X flip). 

The defect movement on the surface is implemented by combining the defect expansion 
and contraction, as shown in Fig. 14.31 At first, we expand a defect D into D' {D G 
D') by the previously mentioned procedure. Second, the defect region D is annihilated, 
and we obtain the eigenvalue (—1)™' of Z{dD). As we already pointed out, the stored 
logical information is unchanged under these procedures up to the logical Pauli X operator 
depending on the measurement outcome m. 

The logical A-basis state preparation is executed by combining the defect creation and 
annihilation, as shown in Fig. 14.41 We first create a defect region Dq = Di + D 2 + D 3 , 
which consists of three adjacent defect regions Di, D 2 , and D^. The two defects Di and 
T >3 are utilized as a defect pair. We can define a logical X operator by employing the 
qubits inside the defect region D 2 . Because all qubits inside the defect region are in the 
1+) state, the eigenvalue of the logical X operator at this time is +1. Then the defect 
region D 2 in-between Di and is annihilated by measuring the star operators. These 
measurements commute with the logical X operator, and hence its eigenvalue is still +1 
after the annihilation. Now we have a defect pair qubit, which is the eigenstate of the 
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Figure 4.4: A preparation of the logical A-basis state by creation of a defect region Dq 
and annihilation of a defect region D 2 between two defect regions Di and D^. Note that 
after the first defect creation, the state is stabilized by the logical X operator A(ci), which 
is untouched during the following defect annihilation. 

logical X operator. A measurement in the logical X basis can be implemented by doing 
the A-basis state preparation in an inverse way. More precisely, two defects Di and D 3 
are connected by making a larger defect Dq = Di + D 2 + D^. We can choose a logical A 
operator such that all constituent qubits belong to the defect region Dq. Then, we obtain 
the eigenvalue of the logical A operator. 

The dual defect creation, annihilation, expansion, contraction, and propagation can be 
done in the same way on the dual lattice with the Hadamard transformation. In Sec. 14.51 
these elementary operations of the primal and dual defects are depicted by a topological 
diagram. 

4.3 Logical CNOT gate by braiding 

Using the elementary operations explained in the previous section, we can perform logical 
gate operations on the defect pair qubits. Indeed, the defects created from the vacuum 
behave like anyonic particles, i.e., braiding a defect around another defect results in a 
nontrivial operation. (It is better to say that we change the manifold continuously to 
deform the ground state, in contrast to the anyons appeared as an excitation on the toric 
code Hamiltonian.) All operations employed in the elementary operations are single-qubit 
or local stabilizer measurements, which can be done by nearest-neighbor two-qubit gates 
and single-qubit measurements. This contrasts with concatenated quantum computation, 
where the logical operations are implemented as transversal operations and hence, essen¬ 
tially, non-local two-qubit gates are employed. 
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Figure 4.5: A logical CNOT gate created by braiding, (a) A primal defect is braided around 
a dual defect. (b,c) The logical X operator is deformed continuously during the braiding 
operations, (d) After the braiding, the logical X operator is represented by the dual chain 
c'l winding around the primal defect. By applying star operators inside the rectangle 
denoted by dotted lines, c[ is decoupled into a chain connecting two primal defects and 
a chain surrounding the dual defect. The corresponding operators are equivalent to the 
products of the logical X operators on the primal and dual qubits. 


We first consider the CNOT gate between primal (control) and dual (target) defect 
pair qubits. Suppose there are primal and dual qubits on the surface code as shown in 
Fig. 14.51 fa), where logical operators are specified by the chains {dD,ci} and {dD,ci} for 
the primal and dual qubits, respectively. We can freely move the defect everywhere we 
want by the defect expansion and contraction. Let us braid the primal defect around the 
dual defect, as shown in Fig. 14.51 (a)-(d). After the braiding operation, the operator X{^^) 
in Fig. 14.51 (d) has the same information as the operator X(ci) before the braiding. Using 
the equivalence relation, 

A(c'i) ~ X{ci)X{dD), (4.4) 

we understand that a correlation is made by the braiding operation between the logical X 
operators X{c\) and X{dD) of the primal and dual defect pair qubits. A similar observa- 
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Figure 4.6: The time evolution of the logical Z operator caused by braiding. The logical 
Z operator on the dual qubit is transformed into the product of the logical Z operators 
on the primal and dual qubits. 

tion holds for the logical operator Z{c\)^ which is transformed into ■Z’(c^) ~ Z{ci)Z{dD) 
(see Fig. 14.61) . On the other hand, Z{dD) and X{dD) are invariant under this operation. 
In short, the braiding operation transforms the logical Pauli operators as follows: 


Z{dD) Z{dD), (4.5) 

Z(ci) ^ Z{ci)Z{dD), (4.6) 

X{dD) X{dD), (4.7) 

X{ci) X{ci)X{dD). (4.8) 


This transformation is equivalent to that for the CNOT gate, Eqs. (|2.14112.171) . Thus, 
braiding the primal defect around the dual defect results in a logical CNOT gate between 
the primal and dual defect pair qubits. 

Unfortunately, in the above CNOT gate, the primal (dual) defect is always a control 
(target) qubit. Such CNOT gates always commute with each other, which is a natural 
consequence of the fact that the defect qubits on the surface code are Abelian. In order to 
realize a genuine CNOT gate between the same type of qubits (and hence noncommuting 
gates), we utilize a teleportation-based gate |GC99] . as shown in Fig. 14.71 Only CNOT 
gates between primal (control) and dual (target) qubits are employed. The Pauli basis 
measurements are also done as mentioned in the previous section. Accordingly, the CNOT 
gate between primal qubits is realized by braiding the primal defects around virtual dual 
defects, which are created and annihilated as ancillae. 

The above braiding operations are implemented with topological quantum error cor¬ 
rection at each elementary step, which will be explained later in detail. All operations 
considered so far can be executed keeping the defect size (circumference) and defect dis- 
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Figure 4.7: The circuit diagram of the teleportation-based CNOT gate between primal 
qubits using a dual qubit as ancilla (left). The control and target qubits of the CNOT gates 
are always primal and dual qubits, respectively. The corresponding braiding operation 
(right). 

tance larger than a length d, which provides a code distance of the logical qubits. Thus, 
the logical error probability for the logical CNOT gates decreases exponentially by in¬ 
creasing the characteristic length d of the system. In other words, the CNOT gates are 
topologically protected. 

We should mention that the braiding operations explained above are not the only 
way to perform fault-tolerant operations on the surface code. There are another ap¬ 
proaches to perform logical operations fault-tolerantly for the encoded degrees of freedom 
of the surface code. One is the lattice surgery scheme |HFDVM12] . where the bound¬ 
ary conditions of two planar surface codes are engineered to perform a logical operation. 
Another is to employ twists, which are topological defects introduced by point defects 
on lattices [BomlOl IBomll| . All Clifford gates can be implemented by the twist cre¬ 
ation, braiding, and annihilation similarly to the defect pair qubits explained above. All 
these different approaches can be view as the logical operations by the code deforma¬ 
tions |DKLP02l IRHG061 IBMD09| and seem to be a unique feature for the topological 
codes contrasting to the transversal logical gate for the concatenated codes. 

4.4 Magic state injections and distillation 

Unfortunately, the topologically protected CNOT gates do not allow universal quantum 
computation, because Clifford circuits can efficiently be simulated classically due to the 
Gottesman-Knill theorem |Got98b] . Here, we explain how to perform single-qubit rota¬ 
tions on the defect pair qubit, while, unfortunately, they are not topologically protected. 

Suppose we have a defect pair qubit, whose logical operators are given by Z{dD) and 
X{ci). We first consider a logical Z rotation defect pair qubit. To 
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Figure 4.8; (a) A magic state injection by contracting the defect into an elementary cell, 
(b) A circuit diagram for the logical Z rotation on the defect of an elementary cell. 

this end, the defect region D is contracted to a single face fi by using the annihilation 
process, as shown in Fig. 14.81 (a). Now the logical operator Z{fi) is a four-body operator. 
(This implies that the code distance becomes 4 at this stage.) A rotation with respect to 
the logical operator Z{fi), g^n be implemented indirectly by using an acilla 

qubit located at the center of the face shown in Fig. 14.81 (b). First, the four CZ gates are 
applied between the ancilla and the four edge qubits on the face. Second, after applying a 
A-rotation Rx{0) = the ancilla qubit is measured in the Z-basis. According to 

the measurement outcome m = 0,1, a logical rotation is implemented. 

After the above procedure, the defect, as a single face /;, is expanded again into the defect 
region D to restore the code distance of the defect pair qubit. 

Next we consider a logical X rotation, In this case, the two defects are 

moved and deformed near each other such that X{di) becomes a one-body operator as 
shown in Fig. 14.91 The rotation is now easily implemented by a single-qubit 

gate. After that, the distance between two defects is restored. 

In this way, we can perform an arbitrary single-qubit unitary operation on the defect 
pair qubit by continuously deforming the defect pair qubit. These deformations and the 
operations for the logical rotations are depicted in Fig. 14.101 Unfortunately, during the 
deformation, the code distance inevitably becomes relatively small, so that we can perform 
logical rotations directly with nearest-neighbor operations. Thus, these processes are not 
topologically protected. 

These non-topological operations can be utilized to inject noisy magic states on the 
surface. Specifically, the logical state is injected as shown in Fig. 14.101 (c). 

Then, these noisy magic states are distilled by topologically protected operations to clean 
magic states, which have a fidelity high enough for reliable quantum computation. Note 
that we are allowed to use only the CNOT gates, and there is no topologically protected 
single-qubit Clifford gate. To manage this, we distill two types of magic states. One is 
the eigenstate of the Pauli-U operator, which is used to implement the S gate via gate 
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Figure 4.9: A magic state injection by connecting two defects. If two defects regions are 
adjacent, the logical X operator is a physical Pauli X operator between two defects. 

teleportation [GC991 IZLC00| . as shown in Fig. 14.111 (a). Another is the eigenstate of 
the {X + Y)ly/2 operator, which is used to implement the vr/S gate, a non-Clifford gate 
necessary for universal quantum computation. The magic state distillation for the Y -basis 
state is executed solely by the topologically protected CNOT gates using the Steane 7- 
qubit code, similarly to the method introduced in Sec. 12.81 The topologically protected 
CNOT gates and the S gates with the distilled T-basis states are further employed to 
distill the {X + y)/\/2-states through the Reed-Muller 15-qubit code, as explained in 
Sec. ESI 

Using these distilled magic states and the CNOT gates, we can implement the single¬ 
qubit gates shown in Fig. 14.111 (a) and (b), which together with the CNOT gate constitute 
a universal set of gates. Accordingly, universal quantum computation is performed reliably 
on the surface code. Note that all operations employed are single-qubit gates, two-qubit 
nearest-neighbor gates, and single-qubit measurements on the 2D array of qubits. 

4.5 Topological calculus 

Based on the microscopic understanding of the topological operations, we can introduce a 
topological diagram and a topological calculus, which allows us a diagrammatic description 
of topological quantum computation on the surface codes. 

In this diagram, the 3D space-time trajectory of the defects is depicted by projecting 
it onto a 2D plane like link diagrams mentioned in Sec. 11.5.31 The trajectory of the primal 
and dual defects as tubes in 3D space-time are denoted by pairs of solid and gray-colored 
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Figure 4.10: (a) The logical Z rotation, (b) The logical X rotation, (c) The logical 
e“*(®/^)^|+) state injection. 
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Figure 4.11: (a,b) Single-qubit rotations by one-bit teleportations [ZLCOOj . 


lines: 


primal defect pairs dual defect pairs 

The braiding operation is denoted by a double crossing as follows: 



The logical X state preparation and X-basis measurement of the primal defect are denoted 
by closures: 
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The logical Z state preparation and Z-basis measurement of the primal defect are denoted 
by endpoints of the tubes: 


The state preparations and measurements for the dual defect pair qubit are depicted in 
a similar way with gray (curved) lines, except for the basis change by the Hadamard 
transformation. 

The trajectories of the logical operators are represented by surfaces in 3D space-time 
like a Seifert surface, which we call correlation surfaces. Specifically, the logical Z operator 
for the primal defect corresponds to a surface wrapping around the primal defect. The 
logical X operator for the primal defect corresponds to a surface whose boundary is the 
primal defect tubes. (This is also the case for the dual defect, except for the Hadamard 
transformation.) For example, the time evolution of the logical X and Z operators by the 
braiding operation can be represented by the following surfaces 


The logical operators before and after the 
boundaries of these surfaces, respectively. We can confirm Eqs. ()4.5l) - (j4.8p from the left 
and right boundaries the correlation surfaces of the above diagrams. 

The CNOT gate between the primal defect pair qubits is depicted as follows: 



operation correspond to the left and right 
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In this diagram, we can easily confirm that the braiding operations transform the logical 
operators according to the rule for the CNOT gate, i.e., X®I —)■ X®X and I®Z —)• Z®Z'. 



As seen above, the action on the code space is defined by the left and right boundaries 
of the correlation surfaces. Thus, the logical action of the topological operation is invariant 
under transformations that do not change the topology of the left and right boundaries 
of the correlation surfaces. Similarly to link diagrams, the boundary topology is invariant 
under the Reidemeister moves: 



In the present case, we have further transformations under which the logical action on 
the code space is invariant due to the properties of the defect pair qubits. At first, the 
following two crossings of tubes of the same type are equivalent: 


(4.9) 

Second, if two defect tubes are a defect pair qubit, a dual loop wrapping around a defect 
pair can be removed: 




4 


(4.10) 
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This is because the primal defect pair qubit is stabilized by a Z-type loop operator sur¬ 
rounding the defect pair, and hence the dual loop is nothing but a trivial measurement of 
the stabilizer operator. Third, because the defect pair qubit is a Z 2 Abelian anyon, if we 
braid a primal defect around a dual defect twice, it results in an identity operation: 



(4.11) 

Fourth, and most importantly, two tubes can be connected or disconnected by the following 
procedure, which we call a 4>-transformation; 



where we note that the dual ring in the middle serves to stop the Z-type correlation surface 
from propagating toward the bottom and also serves to mediate the A-type correlation 
surface toward the right. Specifically, if two defect pair qubits are connected by the <1>- 
transformation between their closures, we can remove the dual ring by using the rule 

dm]): 


DC-y<‘DCZ‘DC 

By using the <h-transformation and rules (j4.1ip and (I4.10p . we can transform the defect 
pair qubit as follows: 
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We can also easily confirm that the logical operators are invariant under this transforma¬ 
tion as follows: 



The <l>-transformation is very useful to simplify a complicated topological operation on 
the surface by joining the tubes as follows: 



where we frequently used the ^>-transformation and rule (jd.lOjl . Here, the dual triple-ring 
wrapping around each of the two defect tubes serves to reflect the logical Z operator from 
the upper to the lower tubes as follows: 
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For example, if we apply this transformation to the CNOT gate between the primal 
defect pair qubits, we obtain a much simpler diagram as follows: 



where the dual ring in the middle serves to keep the code space of the defect pair qubit as 
shown in the left above. It is straightforward to check the transformations of the logical 
operators. In general, by joining the defect tubes directly, we can transform the logical 
operators under multiple CNOT gates. However, if two tubes are joined, the definition 
of the defect pair is broken. In order to keep the defect pair qubit encoding, we need the 
dual ring wrapping around the two tubes at the joint. 

The magic state injection is denoted by the following diagram: 



which can be viewed as a superposition of two correlation surfaces for the two anti¬ 
commuting logical operators. The non-Clifford gate by one-bit teleportation can be de¬ 
scribed in the topological diagram as follows: 
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Figure 4.12: Syndrome measurements for plaquette (left) and star (right) operators, where 
qubits on a face center and on a vertex are employed as ancillae. The alternative syndrome 
measurements for the plaquette and star operators are done repeatedly. 



The right diagram, which is topologically equivalent to the left one, indicates that the 
non-Clifford gate by one-bit teleportation with the magic state injection is equivalent to 
simply performing the non-Clifford gate by the method mentioned in the previous section. 

The present diagrammatic description of the topological operations and correlation 
surfaces on them provide us an intuitive understanding of how topological quantum com¬ 
putation on the surface code is performed. These transformation rules will be useful to 
optimize the complexity (space-time volume required) of the braiding operations. 

4.6 Faulty syndrome measurements and noise thresholds 

We have seen how universal quantum computation is executed on the surface code. At 
each step, we have to perform topological quantum error correction to protect the quan¬ 
tum information encoded by the defects. In Sec. 13.41 we analyzed topological quantum 
error correction on the surface code. However, at that time, we assumed ideal syndrome 
measurements. However, in a fault-tolerant quantum computation, we have to take into 
account all reasonable sources of noise, including faulty syndrome measurements, and the 
quantum computation has to tolerate them as well. Here, we explain how faulty syndrome 
measurements are handled during topological quantum computation, which completes the 
big picture of fault-tolerant topological quantum computation. 
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During the topological operations, the star and plaquette operators are measured in 
the vacuum region to obtain the error syndrome. These measurements are implemented 
using ancillae located on each vertex and face center for the star and plaquette operators, 
respectively (see Fig. I4.12|l . The ancilla state is prepared to be |+) and the A(X) or 
A(Z) gates are performed from the ancilla qubit as a control to the four qubits on dfi 
or 6vk- By measuring the ancilla qubit in the X-basis, we obtain the eigenvalue of the 
star or plaquette operator. If an error is introduced during the measurement, we obtain 
an incorrect syndrome value, which has to be dealt with in the topological quantum error 
correction. 

Again, we assume that the error is given as a Pauli error for simplicity. We also 
assume that the X and Z errors, which might be correlated in general, are corrected 
independently. The former assumption could be justified as follows. Any Kraus operator 
of the noise map can be decomposed into a superposition of Pauli operators. Such a 
superposition is collapsed into Pauli errors by the syndrome measurements, which map 
different Pauli errors (of low weight) into orthogonal subspaces. Note that, in this case, 
we have to model the error per gate carefully. The latter assumption makes the analysis 
very simple, but only results in an underestimation of the noise threshold. Below we only 
consider correction of the Z errors, but it can be applied straightforwardly to the X errors. 

Let cf(t) = be a 1-chain specifying the space-time Z error location, where 

if zi{t) = 1, an Z error occurs on qubit e; at time step t. The Z operators on the code 
state at time t are denoted by ci{t) = zi{t) (i.e., Z[ci(t)]).The state at time t satishes 
the following equation of the motion: 


Cl (t +1) — Ci(t) -|- Ci(t -|- 1). 


(4.13) 


At each time step t, we measure the star operator Bk{t) and obtain the measurement 
outcome 


^k{t) = I 0 Zi{t) 1 ®zl{t), 
Keie5vk 


(4.14) 


where the bit z^{t) E {0,1} indicates an error on the measured syndrome, which we call 
a measurement error. To obtain an equation consisting only of errors, we calculate the 
parity of the syndrome at time steps t and t — 1, 


Sk{t) = mk{t) © mk{t - 1 ) = 


© 2f(‘) 

yeieSvk 


© zl{t) © zl{t — 1). 


(4.15) 


Hence, errors on the measured syndrome can be detected by the parity (difference) of the 
syndrome values at time t and t — 1. We redefine {sfc(t)} as an error syndrome of the 
space-time errors including the measurement error. 

In the case of the perfect syndrome measurement, the error syndrome is given by 
the boundary of errors 5cf, which allows us to use the MWPM algorithm for the error 
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//„ (measurement error) 



Figure 4.13: A mapping from the chain complex in 2D with time evolution to another 
chain complex in 3D. 

correction. What about the case with the faulty syndrome measurement? Using a chain 
complex in a 3D manifold, we can again reformulate the error syndrome {sfc(t)} as the 
boundary of a space-time error chain as shown in Fig. 14.131 We consider a chain complex 
{Co, Cl, ( 72 , 6 * 3 } on a cubic lattice. The basis for the 3-chain is given by a set of cubes 
{qr}- The 3-chain C 3 = € C 3 {zr G {0,1}) is a linear combination over Z 2 . The 

dual chain complex {Cq, Ci, ( 72 , C 3 } is also defined. Specifically, Ci and Cs-j are identified 
via the duality transformation. 

Let us explain how to embed the space-time error location to a 1-chain in the chain 
complex in 3D. The edge e/(f) and its coefficient Zi{t) at time t are mapped into a vertical 
face fm and its coefficient in the 3D chain complex, respectively. The measurement 
error z^it) is mapped to the coefficient z^, of a horizontal face fm'- Then, the space-time 
error location including the measurement errors is described by a 2 -chain or a dual 1 -chain 
in the 3D chain complex: 

C2 = '^Z^fm, ^ ^i='^Z^m&m- (4.16) 

m m 

A syndrome Sk{t) is assigned on each cube qr-, or equivalently each dual vertex in the 
3D chain complex, and denoted by Sr- Now we realize that 

^ ® (4-17) 

fm 

By using the same argument for the 2D case, = 1 iff Vr belongs to the boundary clef 


no 

























of the dual 1-chain cf, specifying the space-time error location. Thus, we can apply the 
MWPM algorithm on the cubic lattice to perform topological quantum error correction in 
space-time. 

To calculate the noise thresholds, we have to model the noise distribution P{c\). First, 
we consider the simplest case, where the errors are located on each dual edge e; of the 
dual cubic lattice with probability p. This means that Z errors occur on each qubit 
independently with probability p at each time step. Moreover, the measured syndrome is 
flipped with probability p. Such a noise model is called a phenomenological noise model. 
Using the MWPM algorithm, the noise threshold has been estimated to be 2.93% |WHP03] . 
The error correction problem under the phenomenological noise model can be mapped into 
a phase transition of the random-plaquette Z 2 gauge model by using the same argument 
as in Sec. 13.51 Specifically, the loop condition for the dual 1-chain (primal 2-chain) in the 
3D model can be solved by introducing a gauge spin on each primal edge ai and defining 
the dual trivial 1-cycle ci = Ylm^rnem G Img(52) via (-1)^"* = Hejea/m where Cm and 
fm are related through the duality relation. In this way, the variable on the face center is 
provided as a product of the gauge spins on the boundary of the face, which leads to the 
random-plaquette Z 2 gauge model. The ordered Higgs and disordered confinement phases 
correspond to fault-tolerant and non-fault-tolerant regions, respectively [DKLP02] . The 
threshold value 2.93% corresponds to the critical point at zero temperature. The optimal 
threshold is provided on the Nishimori line and has been estimated to be 3.3% using the 
Monte Carlo simulation |OAIM04] . 

In fault-tolerant quantum computation, we have to consider any source of noise, in¬ 
cluding the gate operations, during the syndrome measurement. As a standard way to 
model a realistic situation, suppose that each elementary gate is followed by a depolarizing 
channel. This is called a circuit-based noise model. More precisely, an ideal single-qubit 
gate is followed by single-qubit depolarizing noise, 

{l-pi)p+ 

Ae{X,Y,Z} 

An ideal two-qubit gate is followed by two-qubit depolarizing noise, 

{1-P2)p+ Y B)p{A0 B). (4.19) 

A,B£{I,X,Y,Z}\I'S>I 

A faulty Pauli-basis state preparation is modeled by the state 

pi{pp)^{i-ppd-^+Pr^. U20) 

where A = X,Y,Z. A faulty Pauli-basis measurement is modeled by a POVM measure¬ 
ment with POVM elements 

(Pm))-^4 (4-21) 

where 

Pa^P^) ~ ~ Pm) 2 ^ Pm 2 ■ (4.22) 
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Thus, the measurement outcome is flipped with probability pm- 

In Refs. [RHG061IRHG071 IRH07| . the error probabilities were parameterized as pi = 
P 2 = (3/2)pp = {2>/2)pm- The noise threshold with the MWPM algorithm was obtained 
by numerical simulations to be p 2 = 0.75%. In Ref. [WFHll] . the threshold value was 
further improved to ~ 1% by assigning a weight for each edge in MWPM appropriately, 
according to the amount of possible errors. Roughly speaking, the threshold is located 
where the expectation value of each syndrome becomes 0.7: 

((-1)^’-) = 0.7. (4.23) 

In the phenomenological noise model, the expectation value is provided by 

((_1)^.) = (1 _ 2p)\ (4.24) 

which hits 0.7 for p = 2.89%; thus being in good agreement with 2.93%. Moreover, in 
the case of the circuit-based noise model, p is roughly given by Pp + pm + 4:^2 for the 
measurement error, which means that one state preparation, one measurement, and four 
two-qubit gates are employed. The error probability p for the code state is given by 4p2, 
meaning that four two-qubit gates are employed in the syndrome measurements of the 
star and plaquette operators. This yields 

((-1)^’-) = (1 - 2pp - 2pm - 8P2)"(1 - 8P2)^ (4.25) 

which hits 0.70 for p 2 = 0.63% with pi = P 2 = (3/2)pp = {3/2)pm- This is again in a 
good agreement with the numerical result 0.75%. This simple calculation provides rough 
estimates of the threshold values, but there is no validity. If we need a more accurate 
threshold value, we should perform a numerical simulation by taking noise propagation 
and correlation into account [P,HG07( IWFH11| . 


4.7 Experimental progress 


Topologically protected quantum computation in 2D has been utilized as a platform to 
design fault-tolerant architectures for quantum computation. One promising approach is 
on-chip monolithic architectures, such as quantum dots |VMLFYin( l■IVMF'*~12) . silicon- 
based nuclear spins jONR+14) , and superconducting qubits |GFG12l IFMMCl^ lGGM+141 
IGMS"*"!!] . where huge number of qubits are integrated on a single chip and each individ¬ 
ual qubit and the interactions between the qubits are manipulated by multiple lasers or 
electrodes. Among them, the superconducting qubits system is one of the most promis¬ 
ing candidates for implementing topological quantum computation on the surface code, 
because they can be fabricated on a 2D chip, and all elementary operations have al¬ 
ready been demonstrated experimentally |NPT99( lYPA'*~n3l IGNHMO^ IWSB'*~n4] . The 
coherence time and gate fidelity of the superconducting systems have improved rapidly. 
An important breakthrough was made by Martinis’s group at University of Galifornia 
Santa Barbara in 2014, where single-qubit gates with a fidelity of 99.92% and two-qubit 
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gates with a fidelity of 99.4% were demonstrated on a ID array of five superconduct¬ 
ing transmon qubits [BKM"*" fdb IGibl4] . In addition, repetitive quantum non-demolition 
measurements were demonstrated on a ID array of 9 qubits, which improved the fidelity of 
a state preparation by a factor of 8.5, even when using 9 qubits and faulty two-qubit gates 
and measurements |KBF~*~14] . This is an important building block of the fault-tolerant 
quantum computation on the surface code. 

Another approach is the distributed modular architecture, where few-qubit local mod¬ 
ules are connected with quantum channels mediating interactions between separate mod¬ 
ules [LBSBlOl IFTlOl IFYKI121ILB121 INLB13[ lMRR+14] . The few-qubit quantum module 
has already been experimentally realized in various physical systems, such as nitrogen- 
vacancy centers in diamond and trapped ions. Furthermore, entangling operations be¬ 
tween separate local modules have been experimentally demonstrated. This experimental 
and theoretical progress will gradually lead us to large-scale quantum computation. 
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Chapter 5 

Topologically protected MBQC 


In this chapter, we reformulate topological fault-tolerant quantum computation explained 
in the previous chapter in terms of meausrement-based quantum computation. 


5.1 Topological cluster state in 3D 


Consider a (primal) cubic lattice £ and Z 2 chain complex on it, {Cq, Ci, £* 2 , Cs}, where 


Co = ^ ZkVk € Co, Cl = ^ Ziei G Ci, 

k 1 

(5.1) 

^2 — ^ ^ ^mfm ^ C3 — ^ ^ ^nQn ^ 

(5.2) 


m n 


with Zk, zi, Zm, Zn G Z 2 . We also consider a dual cubic lattice C through the relations 
Vk ^ Qk, ei ^ fi, fm ^ em, and qn o Vn- 

Qubits are defined on the edges and faces of the primal lattice £ (or equivalently primal 
and dual edges), as shown in Fig. 15.11 We define an operator A{ci) {i = 1, 2) in terms of 
a 1-chain ci = zjCj or a 2-chain C 2 = Ylj ^jfj ^ 


A{ci) = \{A^K (5.3) 

j 

The stabilizer generators of a 3D cluster state for topologically protected MBQC are 
defined on the primal and dual elementary faces fm, fm' (see Fig. 15.21 (a)): 


Kf^ = Xf^Zidfm), (5.4) 

Xj, = Xf- Zidfm'). (5.5) 

A unit cell of the 3D cluster state is shown in Fig. 15.11 This notion of stabilizer generators 
of the cluster state is quite useful; it provides a connection between the operators and the 
chain complex as follows. By multiplying the two stabilizer operators and Kf^,, we 
have 


= X{fm + fm')Z{d{fm + fm')), 


(5.6) 
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primal cubic lattice 



Figure 5.1: A unit cell of the cluster state for topological MBQC. The primal and dual 
cubes, faces, and edges are also shown. 

(see also Fig. 15.21 (b)). By using this property, we can define a stabilizer operator on a 
2 -chain C 2 , 

Kic2) = n ^/: = Xic2)Z{dc2), (5.7) 

m 

(see also Fig. 15.21 (c)). Furthermore, for the two 2-chains, C 2 and C 2 , we have 

K{c2 + c'^) = K{c2)K{c!^). (5.8) 

Let us see how the 3D cluster state is related to topological quantum computation on 
the surface code explained in Chapter 01 Recall the circuit diagrams for the syndrome 
measurements of the plaqette and star operators in Fig. 14.121 The measurement for the 
plaquette operator is done by applying the CZ gates between the ancilla qubit on the face 
center and the four qubits on the edges. This operation is the same as generation of the 
cluster state stabilized by K{fm) with a horizontal face fm- The syndrome measurement 
for the star operator can be done by the CZ gates with the basis change by the Hadamard 
gates. This corresponds to generation of the cluster state stabilized by K{fi) with a 
horizontal dual face /;. Moreover, the horizontal edge qubits, which constitute the surface 
code, are connected by applying the CZ gates vertically in order to perform the Hadamard 
gates for the basis change. In this way, we recover the 3D cluster state stabilized by 
K{fm) and K{fi) for all primal and dual faces fm and /;. Two of three dimensions are 
employed for the spatial degrees of freedom, constituting the surface code. One is for 
the time evolution of measurement-based quantum computation. The measurements are 
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Figure 5.2: (a) A stabilizer operator defined on a primal face, (b) , = X(/m + 

fm')Z{d{U + fm')). (C) K{c 2) = X{c2)Z{dc2). 


done along the time-like axis, where even and odd layers, corresponding to the syndrome 
measurements of the plaquette and star operators respectively, together with constitute 
an elementary time step of the topologically protected MBQC. Below we will see how the 
topological operations on the surface code are translated into a measurement pattern of 
the MBQC on the 3D cluster state. 

5.2 Vacuum, defect, and singular qubit regions 

The cubic lattice is divided into three regions: the vacuum V, defect V, and singular qubits 
S (the detailed definitions are provided later). In the vacuum region, the topological 
quantum computation is protected through topological quantum error correction. The 
defect regions are utilized to implement topological quantum computation by braiding 
defects. We have two types of defects: the primal {D) and dual {D) defects. For simplicity, 
we only consider the primal defect. The extension to the dual case is straightforward by 
replacing primal by dual in the derivation. The primal defect D is defined as a set of 
primal cubes. The primal face qubits inside the primal defect (except for those on the 
boundary dD) are measured in the Z-basis to remove the corresponding bonds of the 
cluster state (or, equivalently, we can prepare the cluster state without those bonds from 
the beginning). On the boundary dD, the primal face qubits are measured in the X- 
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vacuum 


(b) 



Figure 5.3: (a) A defect pair logical qubit, (b) The logical operators at time 

step t. 


basis. The primal edge qubits belonging to the primal defect (including its boundary) are 
measured in the A-basis. The measurement pattern for the dual defect is defined similarly. 

Unfortunately, only Clifford circuits such as Pauli-basis preparations, measurements, 
and CNOT gates, are implemented in a topologically protected way. For universal quan¬ 
tum computation, magic states for the non-Clifford gates are injected on the singular 
qubits, which are always located in-between two defects. The injections are executed by 
measuring the singular qubits in the Y- and {X -|-y)/\/2-bases. These measurements cor¬ 
respond to injections of (|0)-|-e“*’^/^|l))/-\/2 and {\Qi) + e~™/^\1))/^2 (up to global phases), 
which are utilized to implement the S = and T = gates via gate tele¬ 

portation, respectively. The singular qubits are not topologically protected, because two 
defects are made close to each other resulting in shortening the code distance. However, 
we can obtain clean magic states with topologically protected Clifford gates through the 
magic state distillation protocols |BK05| . In this way, an arbitrary quantum computation 
is executed fault-tolerantly. Below, we will define these three regions more precisely and 
see how topological quantum computation is executed in a measurement-based way. 


5.3 Elementary operations in topological MBQC 

Definition of a logical qnbit 

The logical information is encoded by using a pair of two defects as shown in Fig. 15.31 
(a) and (b), where the measurements are done from left to right. The logical degree of 
information at time step t is described a primal 1-chain cisurrounding the defect and a 
dual 1-chain c\ connecting the two defects, as shown in Fig. l5.3l lab After measuring qubits 
up to the {t — l)th even and odd layers, according to the measurement patterns presented 
before, the following two operators may become logical operators, which commute with 
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Figure 5.4: A logical identity gate. The logical operators at time steps t and t' are related 
by the correlation surfaces K{c 2 ) and K{c 2 ) via the measurements. 

the stabilizer group of the remaining cluster state and are independent of it: 

= Z(ci), = A(^)Z(ci), (5.9) 

where ^ indicates the dual face qubits on the even layer at time step t that are the left 
neighbor of ci , as shown in Fig. 15.31 (b). These two operators anticommute with each 
other and represent a logical qubit. (If the cluster state ends at the even layer at the time 
step t, then the two logical operators are equivalent to the logical operators of the surface 
code. Because there are the time-like CZ gates for the Hadamard gates, the logical X 
operator in Eq. (15.9p accompanied by the Z operators.) 

Identity gate 

Next, we will see how these logical operators evolve with the measurements. We consider 
a correlation surface defined by a primal 2-chain C 2 and a dual 2-chain C 2 as shown in 
Fig. 15.41 A stabilizer operator A( 02 ) on the correlation surface C 2 surrounding the defect 
is obtained by multiplying the stabilizer generators Kf^ on the primal 2-chain C 2 : 

A(c2 ) = n^/: = Z{dc2)X{c2). (5.10) 

m 

Similarly, a stabilizer operator on the dual correlation surface C 2 is defined by multiplying 
Kf on the dual 2-chain C 2 : 

Jm 

A(c 2 ) = llKf2 = Z{dc 2 )X{c 2 ). (5.11) 
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Suppose measurements are done from the left to the right, except for those qnbits on the 
final even layer. Using the correlation surface, we obtain equivalence relations between 
the logical operators at time t and t': 

if ~ Z(cW)K(c2 ) = Z(cP)X(c2). (5.12) 

lJ) ~ X{^^)Z{c’i^)K{c2) = X{K)X{c2)Z{cP) (5.13) 

Here, A ^ B means that A and B are eqnivalent up to a multiplication of the stabilizer 
operator of the clnster state, meaning that both A and B act the same on the cluster 
state. When the qnbits on C 2 have been measnred in the X-basis, we obtain 

L« ~ (5.14) 

L« ~ 4+'), (5.15) 

where we assumed that all measurement outcomes are +1 for simplicity. This relation 
indicates that the logical information at time step t is propagated to time step t' without 
any operation, i.e., a logical identity operation. 

State preparation and measurement 

Next, we consider how the logical qubit is prepared from the vacuum. To prepare the 
eigenstate of 4\ we utilize the defect shown in Fig. [53] (a). By considering the correlation 
snrface C 2 , we obtain 

K{c2) = X{c2)Lf. (5.16) 

Because X(c 2 ) commutes with the measurements, the state at time step t is stabilized by 
4\ and hence a logical Z-basis state is prepared. Considering another snrface C 2 , shown 
in Fig. 15.51 (b), the state at time step t is also stabilized by Z{dc' 2 ) = L'4. Thus the 

pair of the defects is appropriately encoded into the code space. Both and L'4 a-ct 
equivalently as logical Z operators. 
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Figure 5.6: A logical A-basis state preparation. 

Next, we consider the defect shown in Fig. 15.61 Considering a correlation surface C 2 , 
we obtain 


K{c 2) = Z{dc2)X{c2). (5.17) 

After the measurements, the state at time step t is stabilized by = A(^)Z(ci), where 
Cl is a 1-chain on the tth even layer connecting the two defects. Thus, a logical A-basis 
state is prepared. Again, the logical state is stabilized by Z{ci)Z{di), with ci and di being 
a cycle surrounding each defect. Hence, we can choose either = Z{ci) or ^ = Z{c'i) 
to serve as the logical operator. The logical measurements of the defect pair qubits can 
be done with the same defects as the state preparations, but by reversing the time-like 
direction. 

CNOT gate by braiding 

Let us consider primal defects braiding around a dual defect as shown in Fig. 15.71 (a). 
Similarly to the previous case, we calculate the time evolution of logical operators by the 
measurements. The state at time step t is described by and 

corresponding to and respectively. We first consider a correlation 

surface C 2 with respect to as depicted in Fig. [521(b). Similarly to the identity gate, 
is transformed into An interesting thing happens when we consider the correlation 

surface C 2 with respect to as shown in Fig. 15.71 (c). The stabilizer operator on C 2 is 
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Figure 5.7: (a) A diagram for defect braiding for a logical CNOT gate, (b-e) Time 
evolutions of the logical operators and the corresponding correlation surfaces. 


given by 


A:(c') = Z(9c')A(c') 

= Z(c;(*))Z(ci(*'))Z(c'/*'^)Z(c'/)^(c'2), (5.18) 

where c" is a cycle in the dual defect of a loop as shown in Fig. 15.71 (c). Then after the 
measurements, we obtain an equivalence relation, 

Lf^L^pL'^p. (5.19) 

Note that inside the dual defect region, the dual face (primal edge) qubits are measured 

in the Z-basis, and hence we can obtain the eigenvalue of Z{c'(). After a similar argument 
using a defect surface C 2 and with respect to the dual 1 -chain ci and c'^, shown in 
Fig. 15.71 (di and (e), respectively, we obtain 

~ L^P (5.20) 

Lp ~ LplPl (5.21) 

These relations between the logical operators at time steps t and t' are equivalent to those 
for the CNOT gate. Thus, the defect braiding in Fig. 15.71 results in a logical CNOT gate. 
Now we realize that the correlation surface introduced in Chapter S] as a trajectory of the 
logical operator corresponds to the correlation surface defined by the stabilizer operator 
of the cluster state. 
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Figure 5.8: A state injection on the singular qubit. 

A singular qubit injection for magic state distillation 

So far, we have shown that the Clifford circuits, Pauli-basis preparation, measurements, 
and CNOT gate, can all be implemented in a topological way. Unfortunately, these opera¬ 
tions are not enough to generate universal quantum computation. To implement universal 
quantum computation, we inject V- and (X + y)/\/2-basis states by measuring singular 
qubits in the V- and (X + y)/\/2-bases, respectively, as shown in Fig. 15.81 Let us see how 
this measurement works. Similarly to the previous case, we have two correlation surfaces 
C 2 and C 2 : 

K(c 2 ) = Z(dc 2 )X(c 2 ) = Z(dc 2 )X(c 2 \s)X„ (5.22) 

K(c 2 ) = Z(dc 2 )X(c 2 ) = Z,Z(dc 2 \s)X(c 2 ), (5.23) 

where Ag is a Pauli operator on the singular qubit and [•]\s indicates a chain with a 
removal of an element corresponding to the singular qubit. Suppose the singular qubit is 
measured in the T-basis. After the measurements, the state at time step t is stabilized by 

K{c2)K{c2) (5.24) 

Thus, a logical T-basis state is prepared. When the singular qubit is measured in the 
{X + y)/\/2-basis, the state at time step t is stabilized by 

[K{c2) + K{c2)K{c2)]/V2 ~ (L^ + L?)/V2, (5.25) 

which means that a logical {X + y)/\/2-basis state is prepared. 

These states are utilized to implement S, T, HSH, and HTH gates using gate tele¬ 
portation with the CNOT gate. These gates form a universal set of gates. 
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5.4 Topological quantum error correction in 3D 


Next, we will see how topological quantum error correction is done in 3D. Indeed, in the 
3D case, the argument for the noisy syndrome measurements made in Chapter H] becomes 
more simple as follows. All measurements in the vacuum region are done in the A-basis. 
We consider a stabilizer operator on a unit primal cube g, 

K{dqn)= n (5-26) 

fm^d^n 

where there is no Z operator due to <9 o dqn = 0. This implies that, if there is no 
error, the parity of each six A-basis measurement outcomes on the primal cube is always 
even. The errors are described by using a dual 1-chain E = Z(ci). At a unit cube qn 
belonging to dci, we have \qn H dci\ = odd. (Recall that the primal 3-chain and the 
dual 0-chain are identified.) From a set of odd parity cubes dci, we estimate the actual 
location of errors E' = Z{c'i) such that dci = dc'^. If the total of the actual and estimated 
error chains ci -|- results in a trivial cycle, meaning that there is no defect inside the 
cycle, it can be contracted and removed by a continuous deformation. If the total of the 
actual and estimated error chains ci -|- results in a nontrivial cycle, meaning a cycle 
winding around a defect, EE' = Z{ci + c'^) may result in a logical operator. In such a 
case, the topological error correction has failed. This property is completely the same 
as the topological quantum error correction under faulty syndrome measurements argued 
in Sec. 14.61 If the error probability is smaller than a constant value (the threshold), the 
failure probability of the topological quantum error correction decreases exponentially in 
the characteristic size and distance of the defects. 

Inside the defect region, the face qubits are measured in the Z-basis. Especially, 
the Z-basis measurement outcomes near the defect boundary are employed to evaluate 
the correlation surface. Note that these Z-basis measurements and the removal of the 
corresponding bonds of the cluster state can instead be done by generating a cluster 
without connecting the corresponding bonds in advance. In such a case, the errors on the 
Z-basis measurements do not appear. We can obtain an additional parity and Kj ^ 
at the primal and dual faces on the boundary of the defects, respectively. If the errors 
on the face qubits fm and fm’ are suppressed, the errors on the boundary are reduced 
into errors on a toric code on a 2D surface, dD or dD. Again, if the error probability is 
sufficiently smaller than a constant threshold value, we can correct it faithfully. 

The A- and Z-basis state preparations and measurements, and the CNOT gate ob¬ 
tained by braiding, are topologically protected because we can execute these topological 
operations by keeping the defect size and distance larger than an arbitrarily large con¬ 
stant length. Unfortunately, through the state injection, we shrink the defect size into 
an elementary unit cell, where the defect size and distance become very small. Thus, the 
topological protection is broken down around the singular qubit (see Fig. 15.81 (a)). There 
will also be lower weight errors, which effectively increase the logical error probability on 
the injected logical states. However, noisy injected states can be purified by using the 
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topologically protected Clifford gates, the so-called magic state distillation. The V- and 
(X + y)/\/2-basis states are distilled by using the 7-qubit Steane and 15-qubit Reed- 
Muller codes, respectively as explained. The distilled states are utilized to implement 
non-Clifford gates via gate teleportation, as seen before. In this way, universal quantum 
computation is executed with arbitrary accuracy. 


5.5 Applications for MBQC on thermal states 


Topologically protected MBQC in 3D is useful to study the quantum computational ca¬ 
pacity of quantum many-body states at finite temperature. Consider the stabilizer Hamil¬ 
tonian of the 3D cluster state for topological MBQC |RBB03t IRBHr)5[ IBBD'*T)9] : 


Hfr = -J ^K(f) + '^K(f) 

. f f 

The thermal state at temperature T = l/(/3J) is given by 


(5.27) 


Pic = e ^'^f''/Tr[e 


(5.28) 


Using a unitary operator Ucz = W{fm fi) consisting of CZ gates on all nearest- 

neighbor two qubits, the thermal state can be mapped into the thermal state of an 
interaction-free spin model: 

UczPfcUlz = e-^^7Tr[e-^^f], (5.29) 

where 

i7f = - J ^ X, = UczHicUlz- (5.30) 

i 

Because is an interaction-free Hamiltonian, the stabilizer Hamiltonian, which we will 
call the free cluster Hamiltonian, does not undergo any thermodynamic phase transition. 

The thermal state of the free cluster Hamiltonian is given as a product state of the 
single-spin density matrix: 

Pi = e"^^VTr[e"^^f] = (5.31) 

i 

i 

where 


£i{p) = (1 - p)p pZipZi, 


(5.33) 


and ppj = e ^'^'^/(l -|- e ^'^'^). Because Si and Ucz commute, the thermal state of Hie is 
rewritten as 


Pic = UczPiUfjz = 




Ucz{\+){+\r^Ulz = 


W^iippj) 




(5.34) 
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Figure 5.9: A PEPS picture of the cluster state. 


where I'I'sd) is the ground state of FTfc, i.e., the 3D cluster state. This means that the 
thermal state is given as an ideal 3D cluster state, followed by an independent dephasing 
for each qubit with probability ppj = From the argument made in 

the previous section, if p < 2.9 — 3.3% and hence T = l/(/3J) < 0.57 — 0.59, we can 
then perform universal quantum computation reliably on the thermal state at a finite 
temperature, where the errors originating from the thermal excitations are corrected by 
the topological quantum error correction. On the other hand, in the high temperature 
limit T = l/(/3J) oo, the thermal state is given by a completely mixed state, and hence 
MBQC on it can be simulated classically. 

A projected-entangled-pair state (PEPS) picture [VC041IRBH051IBBD+09 allows us to 
obtain a lower bound for the possibility of a classical simulation. In the PEPS picture, the 
3D cluster state is described as follows. A maximally entangled pair |'0mes) = (|0)l+) + 
\l)\—))/y/2 is shared on each bond. On each site consisting of halves of the entangled pair, 
a projection 


| 0 )( 00 .. 0 | + | 1 )( 11 .. 1 | 


(5.35) 


is performed with an appropriate normalization. The resultant state is the 3D cluster 
state. Because the projection and the Z error £i commute, the effect of the thermal 
excitations on the shared entangled state can be determined beforehand: 

T-bond = £a{p)£b{p)\'>pMEs){'lpMEs\ (5.36) 

If p > (2 — y/2)l2, the decohered entangled pair pbond becomes a separable state. If two 
bonds per site are made separable, the 3D cluster state becomes a separable state. A 
sampling on such a resource state can be simulated efficiently classically |BBD'*~n9t IFT14j . 
In this case, the Z error probability per site has to be > \/2—l,i.e.,T = l/(/3) = 5.77J. 
The true critical temperature Tc between the classically simulatable and universal quantum 
computational phases is located in the range 0.59J < < 5.77J. Note that this model 

exhibits a transition of the computational capability, while there is no thermodynamic 
phase transition in the physical system |RBH05[ |BBD~*~09] . 

In classical information processing, a thermodynamic phase transition, or, more pre¬ 
cisely, an ordered phase below a critical temperature is utilized for robust information 
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storage in magnetic storage devises. While there is no such long range order in the pre¬ 
vious model, it is natural to ask whether or not a long range ordered phase is useful to 
enhance the measurement-based quantum computation on many-body thermal states for 
quantum information processing. To address this issue, Fujii et al. proposed an interacting 
cluster Hamiltonian [FNOMl^ . 

ifJ) 

Because interactions are introduced between the cluster stabilizers, this model is mapped 
by Ucz into an Ising model on a 3D lattice: 

Rising = UczH.MIz = Z ^/^/• 

ifJ) 

Thus, it undergoes a thermodynamic phase transition at a finite temperature. The degen¬ 
erate ground states 

Dcz|+)®” andHczl-)®’^ (5.37) 


are also the 3D cluster states, up to the simultaneous spin flipping due to the global 
symmetry. Because the eigenvalues of the cluster stabilizer have a long range order (they 
are likely to be aligned in the same direction) in a ferromagnetic ordered phase, the 
topologically protected MBQC on the symmetry-breaking thermal state has a special 
robustness against the thermal excitations. In Ref. [FNOMf^ . topological quantum error 
correction of this model is mapped to a correlated random plaquette Z 2 -gauge model in 
3D, where the disorder in the signs of the plaquettes has an Ising-type correlation. By 
using this property and the gauge transformation on the Nishimori line [NisOl] . Fujii et al. 
showed that the critical temperature of this model, and hence the threshold temperature 
for topological protection, is equal to the critical temperature of the 3D Ising model, 
which is the unitary equivalent model of the interacting cluster Hamiltonian. This means 
that the critical temperatures for the topological protection and the thermodynamic phase 
transition of the underlying physical system coincides exactly. Due to this fact, we can 
improve the threshold temperature for topological protection by one order of magnitude. 

While the above Hamiltonian employs multi-body interactions, the 3D cluster state 
can be generated from the thermal states of a nearest-neighbor two-body Hamiltonian for 
spin-3/2 and composite spin-1/2 particles via local filtering operations ILBK"*"!! IFM12] . 
Let us consider a system consisting of a spin-3/2 particle located at site r and a composite 
particle of two spin-1/2 particles located at the nearest-neighbor site r-|-i, with i = 1, 2, 3, 
as shown in Fig 15.101 (a). The Hamiltonian is given by 


H = 


^ Sr • (7r+l + Ir+2 + Ir+s) 


where Sr = {S^,Sr,S^) is the spin-3/2 operator of the center particle at the position 
r and Ir+a = ^r+a or Hr-i-a depending on the interaction types (line or dash). Here, 
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Ar+a = (^?+a) ^r+a> ^?+a) By+a = (-^r+a; -^r+a) -S^+a) are two independent spin- 
1/2 operators on the composite particle at the position r -|- a (a = 1,2,3). The above 
Hamiltonian H can be reformulated as 

H = ^l/r = A/2^(f2-S2-l7) 

r r 

where B = B+i + Ir +2 + .fr-i -3 and Tr = Sr + B- The ground state |G) = 0j. \gr) is given 
by Tr = 0, Sr = 3/2, and Ir = 3/2, where Lr{Lr + 1) {L = T, S, I) is the eigenvalue of the 
operator L^. Each center particle in the ground state |G) is filtered by using the POVM 
measurement; 

{F" = {Sf - l/4)/V6} {a = X, y, z). (5.38) 

If the measurement outcome is a = z, we obtain a four-qubit GHZ (Greenberger-Horne- 
Zeilinger) state |GHZ89| as the post-POVM measurement state: 

|GHZr) = ;^(|0 + ++) + |1-))) 

where |1) and |0) are eigenstates of S^ with eigenvalues -1-3/2 and —3/2, respectively, and 
|±) are the eigenstates of or with eigenvalues ±1, respectively. Even if we obtain 
other outcomes, we can transform the post-POVM measurement state to |GHZ^) by local 
operations. The four-qubit GHZ state is subsequently used to construct the 2D honeycomb 
cluster state, which is a universal resource for MBQG, by measuring the operators A^®B^ 
and A^ ® B^ on the bond particle as shown in Pig. 15.101 (a). 

In the case of finite temperature, we have the thermal state (^j. pr with pr = jZ 

instead of the ground state, where Z indicates the partition function and (3 = T~^ for 
a temperature T. Then, the GHZ state becomes a noisy, say thermal, GHZ state, Ur = 

In the low temperature case, the thermal GHZ state is calculated, 
in the leading order, to be £’ 4 (|GHZr)(GHZr|) with 

= (1 - - 3^2 - 3 ^ 3 ) [I] + qi [Zr] 

+Q2 [Zr+a] + qs [ZrZr+a\^ (5.39) 

a=l,2,3 a=l,2,3 

where gi, 52 ) and 53 are the error probabilities as functions of the temperature T, the 
Pauli Z operator Zb on the qubit at the position b, and [C]p = CpC\ respectively. The 
probability of other errors such as ZrZr+aZr+a' is several orders of magnitude smaller 

than gi, 2 , 3 - _ 

To obtain the 3D cluster state for topological MBQG, as done in Ref. jLBK~*~lI] . the 
five-qubit GHZ state IGHZ^) is generated in a similar way by using spin-2 particles and 
composite particles of spin-1/2, as shown in Eig. 15.101 (b). Instead of the spin-2 particles, 
spin-3/2 particles were employed in Ref. [FM12] to obtain the 3D cluster state shown in 
Fig. 15.101 (c). After the filtering operation and local operations, the two four-qubit GHZ 
states are connected to obtain the five-qubit GHZ state for building the 3D cluster state. 
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Figure 5.10: (a) A system consisting of spin-3/2 particles and composite particles of two 
spin-1/2. After the filtering operation on the ground state, we obtain cluster states, (b) 
A system consisting of spin-2 particles and composite particles of two spin-1/2 for the 3D 
cluster state, (c) A system consisting of spin-3/2 particles and composite particles of two 
spin-1/2 for the 3D cluster state. 


By using the threshold for topologically protected MBQC, we can calculate the thresh¬ 
old temperatures T = 0.21A and T = 0.18A for the cases of spin-2 and spin-3/2 center par¬ 
ticles, respectively [LBK"*"!! IFM12] . Accordingly, we can perform fault-tolerant universal 
measurement-based quantum computation even on the thermal states of local two-body 
Hamiltonians at finite temperature. 
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Appendix A 

Fault-tolerant quantum 
computation 


A.l Fault-tolerant syndrome measurements 

In the case of the stabilizer code, the syndrome measurements are done by simply mea¬ 
suring the stabilizer operators. Several QEC gadgets have been proposed to implement 
the stabilizer measurement fault-tolerantly |DS96[ [Knin51 I^e97] . 

DiVincenzo-Shor’s gadget — The first QEC gadget was proposed by David DiVincenzo 
and Peter Shor, who used cat states as ancillae for the syndrome measurement |DS96] . It 
is based on an indirect measurement of the observable A with the eigenvalues ±1: 


1 +) 

l^> 




I±A 


\' 4 >) 


(A.l) 


Eor example, the stabilizer S'! of the seven-qubit code can be measured as the transversal 
X measurements of the corresponding physical qubits in the code block; 


X error 


1+)-^ 

t -< 



f-'^X 

o 





o 

Q. 





0 r 










O 


)' ' 




i — 


two X errors 


(A.2) 
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Unfortunately, this measurement is not fault tolerant, because the errors in the CNOT 
gates [A = X in Eq. (jA.ip ] are spread by the following CNOT gates, as shown in the 
above circuit. To make it fault-tolerant, a cat state |cat) = (|000 • • • 0) -|- |111 • • • l))/^/2 is 
used as an ancilla for the measurement as follows: 


verification 



(A.3) 


where the cat state is verified before connecting with the code state. Because the qubits 
in the code block interact with different ancilla qubits, this measurement does not spread 
the errors in the CNOT gates. Similarly, other stabilizers S 2 , - ■ ■ ,Sq are measured fault- 
tolerantly to obtain the error syndrome. Instead of the verification, one can perform a 
suitable recovery operation by postprocessing the ancilla state after its interaction with 
the code state [DAnT] : 


cat state 



(A.4) 


The DiVincenzo-Shor’s QEC gadget and its improved version both require a lot of physical 
gate operations, which results in deterioration of the performance. 

Steane’s gadget — Subsequently, a relatively simple QEC gadget was proposed by An¬ 
drew Steane |Ste97] . where encoded ancilla states are used to extract the syndrome with 
transversal operations. In particular, for the case of the CSS code, the logical code states 
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can be used as ancilla states. The following circuit executes the Z and X error syndrome 
extractions by using the ancilla |0 l) states, 




|0l)- 


|0l)- 




|0l)- 






|0l)- 




X error detection 


Z error detection 


(A.5) 


Because the ancilla states are the logical code states, one can obtain the error syndrome 
by simply measuring the ancilla states. The syndrome extraction is repeated a couple of 
times to extract reliable error information. An optimized way to extract the syndrome 
information was proposed in Ref. [PVK97] . where the subsequent syndrome extraction 
is performed conditionally according to the preceding syndrome information. For these 
schemes to work fault-tolerantly, the encoded ancilla lO^,) states have to be prepared with 
high fidelity. This is achieved by using either verification or entanglement purification 
[KniOSl IDABOai IADBOK] . 

Knill’s gadget — Another interesting QEC gadget was proposed by Emanuel Knill 
[Kni05] . It is based on quantum teleportation as illustrated in the following circuit: 


I'fpL) 


logical Bell 
state 



logical Bell 
measurement 


(A.6) 


Here, the encoded data qubit {iPl) is teleported to the fresh encoded qubit of the ancilla 
Bell state. Thus, the encoded ancilla Bell state has to be prepared with high fidelity 
by using verification or entanglement purification, similarly to the Steane’s gadget. The 
logical Bell measurement completes the teleportation, namely error-correcting teleporta¬ 
tion. There is no need to identify the error syndrome, but it is sufficient to find the 
logical measurement outcomes of the logical Bell measurement. Thus, it is not neces¬ 
sary to repeat the syndrome extraction in this QEC gadget. The outcome of the Bell 
measurement is properly propagated to the subsequent computation as the Pauli frame 
|Kni05[ IPHNOHEI IDHNO^ . 
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A.2 Fault-tolerant gate operations 

Fault-tolerant computation is now executed by the logical gate operations that are followed 
by the QEC gadgets. This is illustrated for the fault-tolerant CNOT gate as follows: 



(A.7) 

where the code block is depicted as though it is a three-qubit code. A QEC gadget 
is attached to each logical output of the transversal CNOT gate. Because a single error 
never will be propagated as multiple errors in a fault-tolerant gate, a logical error is caused 
by two (or more) simultaneous physical errors. Denoting the number of such faulty pairs 
of error locations as C, the logical error probability is given by Cp‘^, with p being the 
physical error probability. If the physical error probability is sufficiently small, p < IjC, 
one can improve the accuracy of the gate and achieve a fault-tolerant computation. 

A.3 Concatenated quantum computation 

For a reliable computation of a large size, the logical error probability should be reduced 
arbitrarily. This is done by a concatenated fault-tolerant computation [KLZ98bl IKLZ98al 
IAB0971 IABO08j . In the concatenated computation, each physical gate operation is re¬ 
peatedly replaced by a logical gate operation followed by QEC gadgets. 

Suppose A is a quantum computation, which consists of some physical gate operations. 
Then the first level concatenated computation is defined by C(A), where the operation C 
indicates replacing each physical gate with a logical one, followed by QEC gadgets. For 
example, C(CNOT) is described in the diagram (IA.7p . The second level concatenated 
CNOT gate C oC(CNOT) is also described as follows: 



where each physical gate in the QEC gadgets is replaced by the logical one, followed by 
the QEC gadgets. By repeating this procedure, the Th level concatenated computation of 
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A is given by C^(A). Specifically, for a physical gate operation G (e.g., Hadamard, CNOT 
gate, etc.), C\G) is called the level-/ G gate. The /th level concatenated code state is 
called the level-/ qubit, denoted by|0(*^), 

As mentioned previously, the logical error probability of the level-1 gate is given 


by 


pa)=c(p(°))2. 


(A.9) 


where p^^'> = p and C denotes the number of all faulty pairs of error locations. The constant 
C differs between the logical gates. It is, however, sufficient to choose the maximum value. 
Due to the self-similarity of the concatenation, the logical error probability of the level-2 
gate is given in terms of by 

p(2) = C(p(b)2, (A.IO) 

Similarly, the logical error probability of the level-/ gate is obtained recursively as 


p(0 = c(p('-b)2 

= {Gp^^'>f'/G. 


(ATI) 


We conclude that if p^^^ < p^^ = 1/C, the logical error probability can be reduced super- 
exponentially with the concatenation level /. This is the so-called threshold condition. 

On the other hand, the resources usage, consumed for the level-/ gate is estimated 
roughly as 

/jW = N\ (A.12) 

where N indicates the total number of physical gates in the level-1 gate. Suppose that the 
size of the computation is 10”“^ = M. Then an accuracy of p^^ < 10“” is required for 
each logical gate at the highest level. The total resources to perform a reliable quantum 
computation of size M amount to 


( \ log2 N 

iog.„(/p(.))-0 " = 


(A.13) 


where I ~ log 2 [n/log^Q(Cp(*^))] is the number of levels necessary to achieve the required 
accuracy. This result clearly shows that if the physical error probability p^^^ is smaller 
than p^j^, one can execute quantum computation to an arbitrary accuracy with only poly- 
logarithmic overhead. This is the celebrated threshold theorem and the critical value pth is 
called the noise threshold [NCOOl IKit97[ IKLZ98bl IKLZ98al IAB097| . The noise thresholds 
have been calculated to be about 10“^ —10“^ for several fault-tolerant schemes under vary¬ 
ing degrees of assumption and rigor |KLZ98bl IKLZ98al IABQ971IABO081 ISte03l IKniOSl 
I AQpnhi lAcnTl iAOPn8i fAPnl Ist^ ICDTn9| . 
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Appendix B 

Decoding stabilizer codes 


Consider an n-qubit stabilizer code, whose stabilizer group is given hy Q = {Gi}. The 
group consisting of the logical operators is denoted by £ = {Li}. Suppose a Pauli product 
E G {I,X,Y, Z}^'^ acts as an error on the stabilizer code state. The error syndrome, a 
set of eigenvalues of the stabilizer generators, of an error E is denoted by S' = S{E). For 
each error syndrome S, we define a pure error operator R{S) G {/, X, T, Z}®” such that 
S = S[i?(S)]. The pure error operator R{S) is chosen arbitrarily as long as S = S[i?(S)]. 
The error E is decomposed uniquely into logical, stabilizer, and pure error operators; 

E = LiGjR[S{E)]. (B.l) 

We define a decoding map R, which computes the logical operator Li from the error E, 
i.e., Li = V{E). 

The decoding problem consists of finding an optimal logical operator Li that maximizes 
the posterior probability P{L\S): 


argmaxP(Li|S). (B.2) 

Li 

The posterior probability is given by 

P(L,|S) = 1^P(£)5[L,=P(£)]5[S = S(F;)] 

E 

= V E mG,R{s)]. 

Gj&g 

where 1/N is a normalization factor, 5(-• •) is an indicator function, and P{E) is the 
probability of the error E. In general, computing P{L\S) is hard because the summation 
over all stabilizer operators, Y^q ., takes an exponential time. However, if the code has a 
good structure, such as a concatenated code, we can manage it efficiently. 

A concatenated code is defined recursively using the logical qubits at the lower level 
as physical qubits at the higher level. Using the logical Pauli operators G at the 
/cth level {k = 1,2,..), we define a stabilizer group and logical operators 


(B.3) 

(B.4) 
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£(^+1) = at the {k + l)th level. The whole stabilizer group of the concatenated 

code is given by the union of all stabilizer groups where we consider the 1th 

concatenated code. A logical operator of the 1th concatenated code is given by a logical 
operator at the highest level, Lf'^ G At each level, we define a level-lc error . In 
the case of A; = 0, the level-0 error G {/, A, Y, Z}®'^ is a physical error. At a higher 
level, the level-lc error E^^'> is a level-lc logical operator, which will be dehned later. At 
each level, we define an error syndrome , i.e., a set of eigenvalues of 

the level-lc stabilizer operators under the level-(fc — 1) error The level-fc 

pure error operator is defined arbitrarily such that 

The level-lc error E^^'^ is now defined recursively as follows. Any physical error E^^'^ 
can be decomposed into the stabilizer, logical, and pure error operators of level-1: 

£(o) = (B.5) 

The obtained level-1 logical operator is further regarded as 

the level-1 error E^^'^ = Recursively, the level-lc error eW is decomposed into the 

stabilizer, logical, and pure error operators of level-(lc -|- 1): 

Then, we define the level-(lc -|- 1) error ( 5 (^+ 1 )), Let 

be such a level-lc decoding map which computes the level-lc 

logical operator from the level-(lc — 1) error At the highest level, we have a 

decomposition 


£(0) = rW Yl G^'i>R^'^^[S^'^\E^’^-^^)], (B.7) 

k=l 


where E^^'^ is defined by Eq. (IB.Gp . 

The decoding is performed by maximizing the posterior probability of the logical op¬ 
erator Lf\ conditioned on the union of the error syndrome up to the 1th 

concatenation level. 


P(Lf|5«) = l P(eW), 


(B.8) 


where E^^'^ is given by Eq. (IB.71) . Using the hierarchal structure, Eq. (|B.8p . can be rewritten 
as 


P{L? 11 ) 


(B.9) 
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P(5(0) 


H-1 


L?-»eCf-» • ' 


E 




^(i-i)g^(i-i) 

'-l — l 


P(5(z) 15^-1)) 


p(^(J-i)|^a-i)) 


(B.IO) 


If the physical errors occur independently for each qubit, we can factorize the posterior 
probability of the {I — l)th level into posterior probabilities for level-(Z — 1) code blocks 

= (B.ll) 

3 

where is a level-(Z — 1) logical operator acting on the jth code block, and is 

a level-(Z — 1) syndrome with respect to the stabilizer operator on the jth code block. By 
repeating this procedure, we can rewrite the posterior probability Eq. (IB.8h as a summation 
over logical operators defined at each concatenation level G where we have 

dehned the group of the level-Zc logical operators on the jih code block, 

The summation can be reformulated as a marginalization problem on a factor graph, 
which is a bipartite graph consisting of two types of nodes, circles, and boxes [PouOGj ■ 
The variable Xc and a function fb{dh) are assigned on each circle c and box b. Here, db is 
a set of circles neighboring b. Specifically, in the present case, the factor graph is a tree 
graph. The marginal distribution on a tree factor graph can be computed efficiently by 
using the brief-propagation method as follows. From circles to boxes, we pass a message, 

Me—>-fe(®c) — ^6'—>c(®c)- (^-12) 

b'GSc\b 

Then from boxes to circles, we pass another message 

lyb^cixc) = EH fbi^b) (B.13) 

Sb\xc b c'GSb\c 

By repeating these procedures alternatively, we can obtain P{Lf'^\S^^^) as a message 
Vb^eixe), where the fj,c^b{xc) and Vb^eixe) are updated from the bottom leaf nodes to 
the top root node, while using marginal distributions. By replacing ^c^b{xc) and fb{Sb) 
with and we obtain the 

posterior probability P{Lf'^\S^^'^) as Vb^c{xc) at the top root node. In this way, decoding 
by maximizing the posterior probability can be executed efficiently for the concatenated 
stabilizer codes |Pou06] . 
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•^-transformation, 122 
Z2 chain complex, 72 
7-qubit Steane code, 61 
9-qubit Shor code, 57 

Abelian anyon, 102 

Aharonov-Jones-Landau algorithm, 23 

Bloch sphere, 7 
braid group, 25 
byproduct operator, 47, 49 

Calderbank-Shor-Steane code, CSS code, 61 

Chernoff-Hoeffding bound, 20 

classical linear code, 59 

Clifford gate, 8 

Clifford operation, 35 

cluster states, 42 

code distance, 58 

completely-p ositive-t race-preserving (CP TP) 
30 

correlation surface, 120 

defect creation, annihilation, and movement, 
108 

defect pair qubit, 107 

deterministic quantum computation with one 
clean qubit, DQCl, 21 

feedforward, 49 
five-qubit code, 59 

gate teleportation, 48, 49 
Gottesman-Knill theorem, 39 
graph state, 41 

Hadamard test, 19 


homology class, 73 

Jones polynomial, 23, 24 
Jordan-Wigner transformation, 100 

Kauffman bracket, 25 
Kitaev’s compass model, 103 
Kitaev’s phase estimation algorithm, 22 
Kitaev’s toric code Hamiltonian, 101 
Kraus operators, 30 

logical operator, 34, 51, 58 
logical qubit, 58 

magic state, 63 
Majorana fermion, 100 

measurement-based quantum computation, 
46 

minimum-weight perfect matching, MWPM, 
84 

mutual duality relation, 94 

Nishimori line, 93 

one-bit teleportation, 47 

parity check matrix, 60 
path-model representation, 27 
Pauli group, 32 
Pauli operator, 7 
Pauli product, 32 
Pauli-basis measurement, 38 
plaquette and star operators, 77 
positive-operator-valued measure measure¬ 
ment, POVM measurement, 9 
projected-entangled-pair state, PEPS, 146 

quantum error correction, QEC, 55 
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quantum Fourier transform, 22 
quantum Gilbert-Varshamov bound, 93 
qubit, 6 

quotient group, 73 

random-bond Ising model, RBIM, 89 
Reed-Muller code, 64 
Reidemeister move, 121 

Shor’s factorization algorithm, 21 

stabilizer formalism, 32 

stabilizer generator, 33 

stabilizer group, 33 

stabilizer Hamiltonian, 99 

stabilizer QEC code, 58 

stabilizer state, 33 

stabilizer subspace, 33 

symmetry protected topological order, 100 

Temperley-Lieb (TL) algebra, 26 
three-qubit bit flip code, 56 
three-qubit phase flip code, 57 
topological calculus, 119 
topological color code, 96 
topological diagram, 119 
topological order, topological ordered sys¬ 
tem, 99 

topological stabilizer code, 72 
transversality, 57 
trivial cycle, 73 

weight, 58 


138 



References 


[AA06] 

[AAEL07] 

[ABD+09] 

[AB097] 

[ABO08] 

[AC07] 

[ADB05] 

[AGP06] 

[AGP08] 


Dorit Aharonov and Itai Arad. The bqp-hardness of approximating the 
jones polynomial. arXiv preprint quant-ph/0605181, 2006. 

Dorit Aharonov, Itai Arad, Elad Eban, and Zeph Landau. Polynomial 
quantum algorithms for additive approximations of the potts model 
and other points of the tutte plane. arXiv preprint quant-ph/0702008, 

2007. 

Panos Aliferis, Prederico Brito, David P DiVincenzo, John Preskill, 
Matthias Steffen, and Barbara M Terhal. Eault-tolerant computing 
with biased-noise superconducting qubits; a case study. New Journal 
of Physics, 11(1):013061, 2009. 

Dorit Aharonov and Michael Ben-Or. Fault tolerant quantum com¬ 
putation with constant error. In Proc. ACM STOC, pages I76-I88-, 
1997. 

Dorit Aharonov and Michael Ben-Or. Fault-tolerant quantum com¬ 
putation with constant error rate. SIAM Journal on Computing, 
38(4):1207-1282, 2008. 

Panos Aliferis and Andrew W Cross. Subsystem fault tolerance with 
the bacon-shor code. Physical review letters, 98(22):220502, 2007. 

Hans Aschauer, Wolfgang Diir, and H-J Briegel. Multiparticle entan¬ 
glement purification for two-colorable graph states. Physical Review A, 
71(1);012319, 2005. 

Panos Aliferis, Daniel Gottesman, and John Preskill. Quantum accu¬ 
racy threshold for concatenated distance-3 codes. Quant. Inf. Comput., 
6:97, 2006. 

Panos Aliferis, Daniel Gottesman, and John Preskill. Accuracy thresh¬ 
old for postselected quantum computation. Quant. Inf. Comput., 8:181, 

2008. 


139 



[AJL09] 

[AP09] 

[AS06] 

[ASWP13] 

[BAO+12] 

[Bar82] 

[BBC+93] 

[BBC+95] 

[BBD+09] 

[BDCP12] 

[BDSW96] 


Dorit Aharonov, Vaughan Jones, and Zeph Landau. A polynomial 
quantum algorithm for approximating the jones polynomial. Algorith- 
mica, 55(3):395-421, 2009. 

Panos Aliferis and John Preskill. Fibonacci scheme for fault-tolerant 
quantum computation. Physical Review A, 79(1):012332, 2009. 

Pablo Arrighi and Louis Salvail. Blind quantum computation. Inter¬ 
national Journal of Quantum Information, 4(05);883-898, 2006. 

Abbas Al-Shimary, James R Wootton, and Jiannis K Pachos. Lifetime 
of topological quantum memories in thermal environment. New Journal 
of Physics, 15(2):025027, 2013. 

H Bombin, Ruben S Andrist, Masayuki Ohzeki, Helmut G Katzgraber, 
and MA Martin-Delgado. Strong resilience of topological codes to de¬ 
polarization. Physical Review X, 2(2):021004, 2012. 

Francisco Barahona. On the computational complexity of ising spin 
glass models. Journal of Physics A: Mathematical and General, 
15(10):3241, 1982. 

Charles H Bennett, Gilles Brassard, Claude Crepeau, Richard Jozsa, 
Asher Peres, and William K Wootters. Teleporting an unknown quan¬ 
tum state via dual classical and einstein-podolsky-rosen channels. Phys. 
Rev. Lett, 70(13):1895, 1993. 

Adriano Barenco, Charles H Bennett, Richard Cleve, David P DiVin- 
cenzo, Norman Margolus, Peter Shor, Tycho Sleator, John A Smolin, 
and Harald Weinfurter. Elementary gates for quantum computation. 
Physical Review A, 52(5);3457, 1995. 

Sean D. Barrett, Stephen D. Bartlett, Andrew C. Doherty, David Jen¬ 
nings, and Terry Rudolph. Transitions in the computational power of 
thermal states for measurement-based quantum computation. Phys. 
Rev. A, 80:062328, Dec 2009. 

Hector Bombin, Guillaume Duclos-Cianci, and David Poulin. Universal 
topological phase of two-dimensional stabilizer codes. New Journal of 
Physics, 14(7);073048, 2012. 

Charles H Bennett, David P DiVincenzo, John A Smolin, and 
William K Wootters. Mixed-state entanglement and quantum error 
correction. Physical Review A, 54(5);3824, 1996. 


140 



[BE93] 

[BFK09] 

[BH12] 

[BJSll] 

[BK05] 

[BKM+14a] 

[BKM+14b] 

[BMD06] 

[BMD07] 

[BMD09] 

[BomlO] 

[Bomll] 


Hans-Jurgen Briegel and Berthold-Georg Englert. Quantum optical 
master equations; The use of damping bases. Physical Review A, 
47(4);3311, 1993. 

Anne Broadbent, Joseph Fitzsimons, and Elham Kashefi. Universal 
blind quantum computation. In Foundations of Computer Science, 
2009. FOCS’09. 50th Annual IEEE Symposium on, pages 517-526. 
IEEE, 2009. 

Sergey Bravyi and Jeongwan Haah. Magic-state distillation with low 
overhead. Physical Review A, 86(5):052329, 2012. 

Michael J Bremner, Richard Jozsa, and Dan J Shepherd. Classical 
simulation of commuting quantum computations implies collapse of 
the polynomial hierarchy. Proceedings of the Royal Society A: Mathe¬ 
matical, Physical and Engineering Science, 467(2126);459-472, 2011. 

Sergey Bravyi and Alexei Kitaev. Universal quantum computation with 
ideal Clifford gates and noisy ancillas. Phy. Rev. A, 71(2):022316, 2005. 

R Barends, J Kelly, A Megrant, A Veitia, D Sank, E Jeffrey, TC White, 
J Mutus, AG Fowler, B Campbell, et ah Superconducting quan¬ 
tum circuits at the surface code threshold for fault tolerance. Nature, 
508(7497);500-503, 2014. 

R Barends, J Kelly, A Megrant, A Veitia, D Sank, E Jeffrey, TC White, 
J Mutus, AG Fowler, B Campbell, et ah Superconducting quan¬ 
tum circuits at the surface code threshold for fault tolerance. Nature, 
508(7497);500-503, 2014. 

Hector Bombin and Miguel Angel Martin-Delgado. Topological quan¬ 
tum distillation. Physical review letters, 97(18): 180501, 2006. 

H Bombin and MA Martin-Delgado. Topological computation without 
braiding. Physical review letters, 98(16);160502, 2007. 

H Bombin and MA Martin-Delgado. Quantum measurements and gates 
by code deformation. Journal of Physics A: Mathematical and Theo¬ 
retical, 42(9);095302, 2009. 

Hctor Bombin. Topological order with a twist; Ising anyons from an 
abelian model. Physical review letters, 105(3):030403, 2010. 

H Bombin. Clifford gates by code deformation. New Journal of Physics, 
13(4);043005, 2011. 


141 



[BROl] 

[BROS] 

[BR07] 

[BT09] 

[BV93] 

[CDT09] 

[CGC+12] 

[CGM+14] 

[GhiOS] 

[GMS+14] 

[GNHM03] 


Hans J Briegel and Robert Raussendorf. Persistent entanglement in 
arrays of interacting particles. Physical Review Letters, 86(5) :910, 2001. 

Daniel E Browne and Terry Rudolph. Resource-efficient linear optical 
quantum computation. Physical Review Letters, 95(1):010501, 2005. 

Sergey Bravyi and Robert Raussendorf. Measurement-based quantum 
computation with the toric code states. Phy. Rev. A, 76(2):022304, 
2007. 

Sergey Bravyi and Barbara Terhal. A no-go theorem for a two- 
dimensional self-correcting quantum memory based on stabilizer codes. 
New Journal of Physics, 11(4):043029, 2009. 

Ethan Bernstein and Umesh Vazirani. Quantum complexity theory. In 
Proceedings of the twenty-fifth annual ACM symposium on Theory of 
computing, pages 11-20. AGM, 1993. 

Andrew W Gross, David P DiVincenzo, and Barbara M Terhal. A com¬ 
parative code study for quantum fault-tolerance. Quant. Inf. Comput., 
9:0541, 2009. 

Jerry M Ghow, Jay M Gambetta, AD Gorcoles, Seth T Merkel, John A 
Smolin, Ghad Rigetti, S Poletto, George A Keefe, Mary B Roth- 
well, JR Rozen, et al. Universal quantum gate set approaching fault- 
tolerant thresholds with superconducting qubits. Physical review let¬ 
ters, 109(6):060501, 2012. 

Jerry M Ghow, Jay M Gambetta, Easwar Magesan, Srikanth J Srini- 
vasan, Andrew W Gross, David W Abraham, Nicholas A Masluk, 
BR Johnson, Golm A Ryan, and M Steffen. Implementing a strand 
of a scalable fault-tolerant quantum computing fabric. Nat. Commun, 
5:4015, 2014. 

Andrew M Ghilds. Secure assisted quantum computation. Quantum 
Information & Computation, 5(6):456-466, 2005. 

AD Gorcoles, Easwar Magesan, Srikanth J Srinivasan, Andrew W 
Gross, M Steffen, Jay M Gambetta, and Jerry M Ghow. Detecting 
arbitrary quantum errors via stabilizer measurements on a sublattice 
of the surface code. arXiv preprint arXivNflO.hflO, 2014. 

I Ghiorescu, Y Nakamura, GJP Ma Harmans, and JE Mooij. Go- 
herent quantum dynamics of a superconducting flux qubit. Science, 
299(5614):1869-1871, 2003. 


142 



[CR99] 

[CS96] 

[DA07] 

[DAB03] 

[DCPlOa] 

[DCPlOb] 

[Dell4] 

[Deu85] 

[DG97] 

[DG98] 

[DHGB05] 

[DHN06a] 

[DHNOeb] 


William Cook and Andre Rohe. Computing minimum-weight perfect 
matchings. INFORMS Journal on Computing^ 11(2):138-148, 1999. 

A Robert Calderbank and Peter W Shor. Good quantum error- 
correcting codes exist. Physical Review A, 54(2); 1098, 1996. 

David P DiVincenzo and Panos Aliferis. Effective fault-tolerant quan¬ 
tum computation with slow measurements. Physical review letters, 
98(2):020501, 2007. 

Wolfgang Diir, Hans Aschauer, and H-J Briegel. Multiparticle en¬ 
tanglement purification for graph states. Physical review letters, 
91(10):107903, 2003. 

Guillaume Duclos-Cianci and David Poulin. Fast decoders for topolog¬ 
ical quantum codes. Physical review letters, 104(5) :050504, 2010. 

Guillaume Duclos-Cianci and David Poulin. A renormalization group 
decoding algorithm for topological quantum codes. In Information The¬ 
ory Workshop (ITW), 2010 IEEE, pages 1-5. IEEE, 2010. 

Nicolas Delfosse. Decoding color codes by projection onto surface codes. 
Physical Review A, 89(1):012317, 2014. 

David Deutsch. Quantum theory, the church-turing principle and the 
universal quantum computer. Proceedings of the Royal Society of Lon¬ 
don. A. Mathematical and Physical Sciences, 400(1818):97-117, 1985. 

Lu-Ming Duan and Guang-Can Guo. Preserving coherence in quan¬ 
tum computation by pairing quantum bits. Physical Review Letters, 
79(10):1953, 1997. 

Lu-Ming Duan and Guang-Can Guo. Prevention of dissipation with 
two particles. Physical Review A, 57(4):2399, 1998. 

W Diir, M Hein, JI Cirac, and H-J Briegel. Standard forms of 
noisy quantum operations via depolarization. Physical Review A, 
72(5):052326, 2005. 

Christopher M Dawson, Henry L Haselgrove, and Michael A Nielsen. 
Noise thresholds for optical cluster-state quantum computation. Phys¬ 
ical Review A, 73(5):052306, 2006. 

Christopher M Dawson, Henry L Haselgrove, and Michael A Nielsen. 
Noise thresholds for optical quantum computers. Physical review let¬ 
ters, 96(2):020501, 2006. 


143 



[DiV95] 

[DKLP02] 

[DlCDVdNMDll] 

[DN05] 

[DR05] 

[DS96] 

[Easl3] 

[Edm65] 

[EPR35] 

[Eey82] 

[PKLW03] 

[PM 12] 

[EMMC12] 


David P DiVincenzo. Two-bit gates are universal for quantum compu¬ 
tation. Physical Review A, 51(2):1015, 1995. 

Eric Dennis, Alexei Kitaev, Andrew Landahl, and John Preskill. 
Topological quantum memory. Journal of Mathematical Physics^ 
43(9):4452~-4505, 2002. 

G De las Cuevas, W Diir, M Van den Nest, and MA Martin-Delgado. 
Quantum algorithms for classical lattice models. New Journal of 
Physics, 13(9);093021, 2011. 

Christopher M Dawson and Michael A Nielsen. The solovay-kitaev 
algorithm. Quantum Information and Computation, 6:081-095, 2005. 

L-M Duan and R Raussendorf. Efficient quantum computation with 
probabilistic quantum gates. Physical review letters, 95(8):080503, 
2005. 

David P. DiVincenzo and Peter W. Shor. Eault-tolerant error correction 
with efficient quantum codes. Phys. Rev. Lett., 77:3260-3263, 1996. 

Bryan Eastin. Distilling one-qubit magic states into toffoli states. Phys¬ 
ical Review A, 87(3):032321, 2013. 

Jack Edmonds. Paths, trees, and flowers. Canadian Journal of math¬ 
ematics, 17(3):449-467, 1965. 

Albert Einstein, Boris Podolsky, and Nathan Rosen. Can quantum- 
mechanical description of physical reality be considered complete? 
Physical review, 47(10):777, 1935. 

Richard P Eeynman. Simulating physics with computers. International 
journal of theoretical physics, 21(6):467-488, 1982. 

Michael Preedman, Alexei Kitaev, Michael Larsen, and Zhenghan 
Wang. Topological quantum computation. Bulletin of the American 
Mathematical Society, 40(l):31-38, 2003. 

Keisuke Pujii and Tomoyuki Morimae. Topologically protected 
measurement-based quantum computation on the thermal state of a 
nearest-neighbor two-body hamiltonian with spin-3/2 particles. Phys¬ 
ical Review A, 85(1):010304, 2012. 

Austin G Eowler, Matteo Mariantoni, John M Martinis, and Andrew N 
Cleland. Surface codes: Towards practical large-scale quantum com¬ 
putation. Physical Review A, 86(3):032324, 2012. 


144 



[FNIK14] 


[FN0M13] 


[FSG09] 


[FTIO] 


[FT12] 


[FT 14] 


[FTP+05] 


[FWH12] 


[FYlOa] 


[FYlOb] 


[FYKI12] 


Keisuke Fujii, Makoto Negoro, Nobuyuki Imoto, and Masahiro Kita¬ 
gawa. Measurement-free topological protection using dissipative feed¬ 
back. arXiv preprint arXiv:1401.6350, 2014. 

Keisuke Fujii, Yoshifumi Nakata, Masayuki Ohzeki, and Mio Murao. 
Measurement-based quantum computation on symmetry breaking ther¬ 
mal states. Phys. Rev. Lett., 110:120502, Mar 2013. 

Austin G Fowler, Ashley M Stephens, and Peter Groszkowski. High- 
threshold universal quantum computation on the surface code. Physical 
Review A, 80(5):052312, 2009. 

Keisuke Fujii and Yuuki Tokunaga. Fault-tolerant topological one¬ 
way quantum computation with probabilistic two-qubit gates. Physical 
review letters, 105(25):250503, 2010. 

Keisuke Fujii and Yuuki Tokunaga. Error and loss tolerances of surface 
codes with general lattice structures. Physical Review A, 86(2):020303, 
2012. 

Keisuke Fujii and Shuhei Tamate. Computational quantum- 
classical boundary of commuting quantum circuits. arXiv preprint 
arXiv:U06.6932, 2014. 

Paolo Facchi, Shuichi Tasaki, Saverio Pascazio, Hiromichi Nakazato, 
Akira Tokuse, and DA Lidar. Control of decoherence; analysis 
and comparison of three different strategies. Physical Review A, 
71(2):022302, 2005. 

Austin G Fowler, Adam C Whiteside, and Lloyd CL Hollenberg. To¬ 
wards practical classical processing for the surface code. Physical review 
letters, 108(18):180501, 2012. 

Keisuke Fujii and Katsuji Yamamoto. Cluster-based architecture for 
fault-tolerant quantum computation. Phys. Rev. A, 81:042324, Apr 
2010. 

Keisuke Fujii and Katsuji Yamamoto. Topological one-way quantum 
computation on verified logical cluster states. Phys. Rev. A, 82:060301, 
Dec 2010. 

Keisuke Fujii, Takashi Yamamoto, Masato Koashi, and Nobuyuki 
Imoto. A distributed architecture for scalable quantum computation 
with realistically noisy devices. arXiv preprint arXiv: 1202.6588, 2012. 


145 



[GC99] 

[GFG12] 

[GHZ89] 

[Gibl4] 

[Got97] 

[Got98a] 

[Got98b] 

[Gro96] 

[GW09] 

[Har04] 

[HGG09] 

[HGEK14] 

[HDE+06] 


Daniel Gottesman and Isaac L Ghuang. Demonstrating the viability 
of universal quantum computation using teleportation and single-qubit 
operations. Nature, 402(6760):390-393, 1999. 

Joydip Ghosh, Austin G Eowler, and Michael R Geller. Surface code 
with decoherence: An analysis of three superconducting architectures. 
Physical Review A, 86(6);062318, 2012. 

Daniel M Greenberger, Michael A Horne, and Anton Zeilinger. Going 
beyond bells theorem. In Bells theorem, quantum theory and concep¬ 
tions of the universe, pages 69-72. Springer, 1989. 

Elizabeth Gibney. Quantum computer quest. Nature, 516(7529):24-26, 
2014. 

D Gottesman. Stabilizer codes and quantum error correetion. PhD 
thesis, Galifornia Institute of Technology, 1997. 

Daniel Gottesman. The heisenberg representation of quantum comput¬ 
ers. arXiv preprint quant-ph/9807006, 1998. 

Daniel Gottesman. The heisenberg representation of quantum comput¬ 
ers. arXiv preprint quant-ph/9807006, 1998. 

Lov K Grover. A fast quantum mechanical algorithm for database 
search. In Proceedings of the twenty-eighth annual ACM symposium on 
Theory of computing, pages 212-219. ACM, 1996. 

Zheng-Cheng Gu and Xiao-Gang Wen. Tensor-entanglement-filtering 
renormalization approach and symmetry-protected topological order. 
Physical Review B, 80(15) :155131, 2009. 

J. W. Harrington. Analysis of quantum error-correcting codes: sym- 
plectic lattiee codes and torie codes. PhD thesis, California Institute of 
Technology, 2004. 

Alioscia Hamma, Claudio Castelnovo, and Claudio Chamon. Toric- 
boson model; Toward a topological quantum memory at finite temper¬ 
ature. Physieal Review B, 79(24);245122, 2009. 

M Herold, ET Campbell, J Eisert, and MJ Kastoryano. Cellular- 
automaton decoders for topological quantum memories. arXiv preprint 
arXiv:1406.2338, 2014. 

M Hein, W Diir, J Eisert, R Raussendorf, M Van der Nest, and 
HJ Briegel. Quantum computers, algorithms and chaos. In Inter¬ 
national Sehool of Physics Enrico Fermi, volume 162, 2006. 


146 



[HFDVM12] 

[HLGll] 

[HPPOl] 

[HTPV08] 

[Huk99] 

[ICBB14] 

[IP 13] 
[JonlSa] 
[JonlSb] 
[JVMF+12] 

[JW28] 

[JW09] 

[Kas61] 


Clare Horsman, Austin G Fowler, Simon Devitt, and Rodney Van Me¬ 
ter. Surface code quantum computing by lattice surgery. New Journal 
of Physics, 14(12):123011, 2012. 

Min-Hsiu Hsieh and Franqois Le Gall. Np-hardness of decoding quan¬ 
tum error-correction codes. Physical Review A, 83(5);052331, 2011. 

Andreas Honecker, Marco Picco, and Pierre Pujol. Universality class 
of the nishimori point in the 2d±j random-bond ising model. Physical 
review letters, 87(4):047201, 2001. 

Martin Hasenbusch, Francesco Parisen Toldin, Andrea Pelissetto, and 
Ettore Vicari. Universal dependence on disorder of two-dimensional 
randomly diluted and random-bondztj ising models. Physical Review 
E, 78(1):011110, 2008. 

Koji Hukushima. Domain-wall free energy of spin-glass models: 
Numerical method and boundary conditions. Physical Review E, 
60(4):3606, 1999. 

S Iblisdir, M Cirio, O Boada, and GK Brennen. Low depth quantum 
circuits for ising models. Annals of Physics, 340(1):205-251, 2014. 

Pavithran Iyer and David Poulin. Hardness of decoding quantum sta¬ 
bilizer codes. arXiv preprint arXiv: 1310.3235, 2013. 

Cody Jones. Low-overhead constructions for the fault-tolerant toffoli 
gate. Physical Review A, 87(2):022328, 2013. 

Cody Jones. Multilevel distillation of magic states for quantum com¬ 
puting. Phys. Rev. A, 87:042305, Apr 2013. 

N Cody Jones, Rodney Van Meter, Austin G Fowler, Peter L McMa¬ 
hon, Jungsang Kim, Thaddeus D Ladd, and Yoshihisa Yamamoto. 
Layered architecture for quantum computing. Physical Review X, 
2(3):031007, 2012. 

P Jordan and Eugene P Wigner. About the pauli exclusion principle. 
Zeitschrift fur Physik, 47:631, 1928. 

Stephen P Jordan and Pawel Wocjan. Estimating Jones and homfly 
polynomials with one clean qubit. Quantum Information & Computa¬ 
tion, 9:264, 2009. 

Pieter W Kasteleyn. The statistics of dimers on a lattice: 1. the number 
of dimer arrangements on a quadratic lattice. Physica, 27(12):1209- 
1225, 1961. 


147 



[Kat95] 

[Kayl4] 

[KBF+14] 

[KCS14] 

[Kit95a] 

[Kit95b] 

[Kit97] 

[KitOl] 

[Kit03] 

[Kit06] 

[KL98] 

[KL12] 

[KLZ98a] 

[KLZ98b] 


Tosio Kato. Perturbation theory for linear operators, volume 132. 
springer, 1995. 

Alastair Kay. Implications of ignorance for quantum-error-correction 
thresholds. Physical Review A, 89(3):032328, 2014. 

J. Kelly, R. Barends, A. G. Fowler, A. Megrant, E. Jeffrey, T. C. White, 
et al. State preservation by repetitive error detection in a supercon¬ 
ducting quantum circuit. arXiv preprint arXiv:1408.0959, 2014. 

Eliot Kapit, John T Chalker, and Steven H Simon. Passive correction 
of quantum logical errors in a driven, dissipative system: a blueprint 
for an analog quantum code fabric. arXiv preprint arXiv: 1408.0959, 
2014. 

A Yu Kitaev. Quantum measurements and the abelian stabilizer prob¬ 
lem. arXiv preprint quant-ph/9511026, 1995. 

A Yu Kitaev. Quantum measurements and the abelian stabilizer prob¬ 
lem. arXiv preprint quant-ph/9511026, 1995. 

A Yu Kitaev. Quantum computations: algorithms and error correction. 
Russian Mathematical Surveys, 52(6):1191-1249, 1997. 

A Yu Kitaev. Unpaired majorana fermions in quantum wires. Physies- 
Uspekhi, 44(10S):131, 2001. 

A Yu Kitaev. Eault-tolerant quantum computation by anyons. Annals 
of Physics, 303(l):2-30, 2003. 

Alexei Kitaev. Anyons in an exactly solved model and beyond. Annals 
of Physics, 321(1):2~111, 2006. 

E. Knill and R. Laflamme. Power of one bit of quantum information. 
Phys. Rev. Lett., 81:5672-5675, Dec 1998. 

Kao-Yueh Kuo and Chung-Chin Lu. On the hardness of decoding quan¬ 
tum stabilizer codes under the depolarizing channel. In Information 
Theory and its Applications (ISITA), 2012 International Symposium 
on, pages 208-211. IEEE, 2012. 

Emanuel Knill, Raymond Laflamme, and Wojciech H Zurek. Resilient 
quantum computation. Science, 279(5349):342-345, 1998. 

Emanuel Knill, Raymond Laflamme, and Wojciech H Zurek. Resilient 
quantum computation: error models and thresholds. Proceedings of 
the Royal Society of London. Series A: Mathematical, Physical and 
Engineering Sciences, 454(1969):365-384, 1998. 


148 



[Kni05] 

[Kol09] 

[KR97] 

[KW41] 

[LB 12] 

[LBK+11] 

[LBSBIO] 

[LCW98a] 

[LCW98b] 

[LCYPP15] 

[Lic97] 

[Llo99] 

[LMPZ96] 


Emanuel Knill. Quantum computing with realistically noisy devices. 
Nature, 434(7029):39-44, 2005. 

Vladimir Kolmogorov. Blossom v: a new implementation of a min¬ 
imum cost perfect matching algorithm. Mathematical Programming 
Computation, l(l):43-67, 2009. 

Naoki Kawashima and Heiko Rieger. Finite-size scaling analysis of 
exact ground states for±j spin glass models in two dimensions. EPL 
(Europhysics Letters), 39(1):85, 1997. 

Hendrik A Kramers and Gregory H Wannier. Statistics of the two- 
dimensional ferromagnet. part i. Physical Review, 60(3):252, 1941. 

Ying Li and Simon C Benjamin. High threshold distributed quan¬ 
tum computing with three-qubit nodes. New Journal of Physics, 
14(9):093008, 2012. 

Ying Li, Daniel E Browne, Leong Chuan Kwek, Robert Raussendorf, 
and Tzu-Chieh Wei. Thermal states as universal resources for quan¬ 
tum computation with always-on interactions. Physical review letters, 
107(6):060501, 2011. 

Ying Li, Sean D Barrett, Thomas M Stace, and Simon C Ben¬ 
jamin. Fault tolerant quantum computation with nondeterministic 
gates. Physical review letters, 105(25) :250502, 2010. 

D. A. Lidar, 1. L. Chuang, and K. B. Whaley. Decoherence-free sub¬ 
spaces for quantum computation. Phys. Rev. Lett., 81:2594-2597, 1998. 

Daniel A Lidar, Isaac L Chuang, and K Birgitta Whaley. Decoherence- 
free subspaces for quantum computation. Physical Review Letters, 
81(12):2594, 1998. 

Olivier Landon-Cardinal, Beni Yoshida, David Poulin, and John 
Preskill. Perturbative instability of quantum memory based on effective 
long-range interactions. Phys. Rev. A, 91:032303, Mar 2015. 

WB Raymond Lickorish. An introduction to knot theory, volume 175. 
Springer Science & Business Media, 1997. 

Seth Lloyd. Obituary: Rolf landauer (1927-99). Aature, 400(6746):720- 
720, 1999. 

Raymond Laflamme, Cesar Miquel, Juan Pablo Paz, and Wojciech Hu¬ 
bert Zurek. Perfect quantum error correcting code. Physical Review 
Letters, 77(1):198, 1996. 


149 



[MBFIO] 

[MC02] 

[MEK12] 

[MF12] 

[MF13] 

[MFF14] 

[MFI14] 

[MNFT14] 

[MRR+14] 

[MS 77] 
[NCOO] 
[Nie04] 
[Nie06] 


Yuichiro Matsuzaki, Simon C Benjamin, and Joseph Fitzsimons. Prob¬ 
abilistic growth of large entangled states with low error accumulation. 
Physical review letters, 104(5) :050501, 2010. 

F Merz and JT Chalker. Two-dimensional random-bond ising model, 
free fermions, and the network model. Physieal Review B, 65(5):054425, 
2002. 

Adam M Meier, Bryan Eastin, and Emanuel Knill. Magic-state distil¬ 
lation with the four-qubit code. arXiv preprint arXiv:1204-4^21, 2012. 

Tomoyuki Morimae and Keisuke Fujii. Blind topological measurement- 
based quantum computation. Nature communications, 3:1036, 2012. 

Tomoyuki Morimae and Keisuke Fujii. Blind quantum computation 
protocol in which alice only makes measurements. Phy. Rev. A, 
87(5):050301, 2013. 

Tomoyuki Morimae, Keisuke Fujii, and Joseph F. Fitzsimons. Hardness 
of classically simulating the one-clean-qubit model. Phys. Rev. Lett., 
112:130502, Apr 2014. 

Akira Matsuo, Keisuke Fujii, and Nobuyuki Imoto. Quantum algorithm 
for an additive approximation of ising partition functions. Phys. Rev. 
A, 90:022304, Aug 2014. 

Tomoyuki Morimae, Harumichi Nishimura, Keisuke Fujii, and Shuhei 
Tamate. Classical simulation of dqcl _2 or dqc2 _1 implies collapse of 
the polynomial hierarchy. arXiv preprint arXiv:1409.6777, 2014. 

C Monroe, R Raussendorf, A Ruthven, KR Brown, P Maunz, L-M 
Duan, and J Kim. Large-scale modular quantum-computer architecture 
with atomic memory and photonic interconnects. Physical Review A, 
89(2):022317, 2014. 

Florence Jessie MacWilliams and Neil James Alexander Sloane. The 
theory of error correcting codes, volume 16. Elsevier, 1977. 

Michael A Nielsen and Isaac L Chuang. Quantum computation and 
quantum information. Cambridge university press, 2000. 

Michael A Nielsen. Optical quantum computation using cluster states. 
Physical review letters, 93(4):040503, 2004. 

M. A. Nielsen. Cluster-state quantum computation. Rep. on Math. 
Phys., 57:147, 2006. 


150 



[Nis79] 

[Nis81] 

[NisOl] 

[NLB13] 

[NN02] 

[NPT99] 

[OAIM04] 

[OF12] 

[Ohz09a] 

[Ohz09b] 

[Ohzl2] 

[Ohzl3] 

[ONR+14] 


Hidetoshi Nishimori. Conjecture on the exact transition point of the 
random ising ferromagnet. Journal of Physics C: Solid State Physics, 
12(23):L905, 1979. 

Hidetoshi Nishimori. Internal energy, specific heat and correlation func¬ 
tion of the bond-random ising model. Progress of Theoretical Physics, 
66(4):1169-1181, 1981. 

Hidetoshi Nishimori. Statistical physics of spin glasses and information 
processing, volume 187. Oxford University Press Oxford, 2001. 

Naomi H Nickerson, Ying Li, and Simon C Benjamin. Topological 
quantum computing with a very noisy network and local error rates 
approaching one percent. Nature communications, 4:1756, 2013. 

Hidetoshi Nishimori and Koji Nemoto. Duality and multicritical point 
of two-dimensional spin glasses. Journal of the Physical Society of 
Japan, 71(4):1198-1199, 2002. 

Yu Nakamura, Yu A Pashkin, and JS Tsai. Coherent control of 
macroscopic quantum states in a single-cooper-pair box. Nature, 
398(6730) :786-788, 1999. 

Takuya Ohno, Gaku Arakawa, Ikuo Ichinose, and Tetsuo Matsui. Phase 
structure of the random-plaquette z2 gauge model: accuracy threshold 
for a toric quantum memory. Nuclear Physics B, 697(3) :462~480, 2004. 

Masayuki Ohzeki and Keisuke Fujii. Duality analysis on random planar 
lattices. Physical Review E, 86(5):051121, 2012. 

Masayuki Ohzeki. Accuracy thresholds of topological color codes 
on the hexagonal and square-octagonal lattices. Physical Review E, 
80(1):011141, 2009. 

Masayuki Ohzeki. Locations of multicritical points for spin glasses on 
regular lattices. Physical Review E, 79(2) :021129, 2009. 

Masayuki Ohzeki. Error threshold estimates for surface code with loss 
of qubits. Physical Review A, 85(6):060301, 2012. 

Masayuki Ohzeki. Spin glass: A bridge between quantum comput¬ 
ing and statistical mechanics. In Mikio Nakahara and Shu Tanaka, 
editors. Lectures on Quantum Computing, Thermodynamics and Sta¬ 
tistical Physics, volume 8. World Scientific, 2013. 

Joe O’Gorman, Naomi H Nickerson, Philip Ross, John JL Morton, and 
Simon C Benjamin. A silicon-based surface code quantum computer. 
arXiv preprint arXiv.JfOd.Slfd, 2014. 


151 



[Pacl2] 

[PB02] 

[PBKLO04] 

[PHWL13] 

[Pou06] 

[PTBOlO] 

[PVK97] 

[Rau03] 

[RBOl] 

[RBB03] 

[RBH05] 

[RDQDS99] 

[Rei05] 


Jiannis K Pachos. Introduction to topological quantum computation. 
Cambridge University Press, 2012. 

Francesco Petruccione and Heinz-Peter Brener. The theory of open 
quantum systems. Oxford Univ. Press, 2002. 

David Poulin, Robin Blume-Kohout, Raymond Laflamme, and Harold 
Ollivier. Exponential speedup with a single bit of quantum informa¬ 
tion: Measuring the average fidelity decay. Physieal review letters, 
92(17):177906, 2004. 

Fabio L Pedrocchi, Adrian Hutter, James R Wootton, and Daniel Loss. 
Enhanced thermal stability of the toric code through coupling to a 
bosonic bath. Physical Review A, 88(6):062313, 2013. 

David Poulin. Optimal and efficient decoding of concatenated quantum 
block codes. Physieal Review A, 74(5);052333, 2006. 

Frank Pollmann, Ari M Turner, Erez Berg, and Masaki Oshikawa. En¬ 
tanglement spectrum of a topological phase in one dimension. Physieal 
Review B, 81(6):064439, 2010. 

MB Plenio, V Vedral, and PL Knight. Conditional generation of er¬ 
ror syndromes in fault-tolerant error correction. Physieal Review A, 
55(6);4593, 1997. 

R Raussendorf. Measurement-based quantum computation with cluster 
states. PhD thesis, Ludwig-Maximillians Universitat Miinchen, 2003. 

Robert Raussendorf and Hans J Briegel. A one-way quantum computer. 
Phys. Rev. Lett, 86(22):5188-5191, 2001. 

Robert Raussendorf, Daniel E Browne, and Hans J Briegel. 
Measurement-based quantum computation on cluster states. Physieal 
review A, 68(2):022312, 2003. 

Robert Raussendorf, Sergey Bravyi, and Jim Harrington. Long-range 
quantum entanglement in noisy cluster states. Phys. Rev. A, 71:062313, 
Jun 2005. 

FDA Aarao Reis, SLA De Queiroz, and Raimundo R Dos Santos. Uni¬ 
versality, frustration, and conformal invariance in two-dimensional ran¬ 
dom ising magnets. Physieal Review B, 60(9):6740, 1999. 

Ben W Reichardt. Quantum universality from magic states distillation 
applied to css codes. Quantum Information Proeessing, 4(3):251-264, 
2005. 


152 



[RH07] 

[RHG06] 

[RHG07] 

[RWH+12] 

[RZBB94] 

[SBIO] 

[SBD09] 

[Sho94] 

[Sho95] 

[Sho97] 

[SJ08] 

[STD07] 


Robert Raussendorf and Jim Harrington. Fault-tolerant quantum com¬ 
putation with high threshold in two dimensions. Phys. Rev. Lett., 
98:190504, May 2007. 

Robert Raussendorf, Jim Harrington, and Kovid Goyal. A fault- 
tolerant one-way quantum computer. Annals of physics, 321(9):2242- 
2270, 2006. 

Robert Raussendorf, Jim Harrington, and Kovid Goyal. Topological 
fault-tolerance in cluster state quantum computation. New Journal of 
Physics, 9(6):199, 2007. 

Beat Rothlisberger, James R Wootton, Robert M Heath, Jiannis K 
Pachos, and Daniel Loss. Incoherent dynamics in the toric code subject 
to disorder. Physical Review A, 85(2):022313, 2012. 

Michael Reck, Anton Zeilinger, Herbert J Bernstein, and Philip 
Bertani. Experimental realization of any discrete unitary operator. 
Physical Review Letters, 73(1):58, 1994. 

Thomas M Stace and Sean D Barrett. Error correction and degeneracy 
in surface codes suffering loss. Physical Review A, 81(2):022317, 2010. 

Thomas M Stace, Sean D Barrett, and Andrew C Doherty. Thresholds 
for topological codes in the presence of loss. Physical review letters, 
102(20):200501, 2009. 

P. W. Shor. Algorithms for quantum computation: discrete logarithms 
and factoring. In Proceedings of the 35th Annual Symposium on Foun¬ 
dations of Computer Science, SFGS ’94, pages 124-134, Washington, 
DC, USA, 1994. IEEE Computer Society. 

Peter W. Shor. Scheme for reducing decoherence in quantum computer 
memory. Phys. Rev. A, 52:R2493-R2496, 1995. 

Peter W Shor. Polynomial-time algorithms for prime factorization and 
discrete logarithms on a quantum computer. SIAM journal on com¬ 
puting, 26(5): 1484-1509, 1997. 

Peter W Shor and Stephen P Jordan. Estimating Jones polynomials 
is a complete problem for one clean qubit. Quantum Information & 
Computation, 8(8):681-714, 2008. 

Krysta M. Svore, Barbara M. Terhal, and David P. DiVincenzo. 
Noise threshold for a fault-tolerant two-dimensional lattice architec¬ 
ture. Quant. Inf. Comp., 7:297-318, 2007. 


153 



[Ste96] 

[Ste97] 

[Ste99] 

[Ste03] 

[TD04] 

[TL71] 

[VC04] 

[VdNDDM04] 

[VGW96] 

[VKL99a] 

[VKL99b] 

[VL98] 

[VMLFYIO] 


Andrew M Steane. Error correcting codes in quantum theory. Physical 
Review Letters, 77(5) :793, 1996. 

Andrew M Steane. Active stabilization, quantum computation, and 
quantum state synthesis. Physical Review Letters, 78(11):2252, 1997. 

Andrew M Steane. Efficient fault-tolerant quantum computing. Nature, 
399(6732):124-126, 1999. 

Andrew M Steane. Overhead and noise threshold of fault-tolerant quan¬ 
tum error correction. Physical Review A, 68(4):042322, 2003. 

B. Terhal and D DiVincenzo. Adaptive quantum computation, constant 
depth quantum circuits and arthur-merlin games. Quantum Informa¬ 
tion & Computation, 4(2): 134-145, 2004. 

Harold NV Temper ley and Elliott H Lieb. Relations between 
the’percolation’and’colouring’problem and other graph-theoretical 
problems associated with regular planar lattices: some exact results 
for the’percolation’problem. Proceedings of the Royal Society of Lon¬ 
don. A. Mathematical and Physical Sciences, 322(1549):251-280, 1971. 

Frank Verstraete and J Ignacio Cirac. Valence-bond states for quantum 
computation. Phy. Rev. A, 70(6):060302, 2004. 

Maarten Van den Nest, Jeroen Dehaene, and Bart De Moor. Graphi¬ 
cal description of the action of local Clifford transformations on graph 
states. Physical Review A, 69(2):022316, 2004. 

Lev Vaidman, Lior Goldenberg, and Stephen Wiesner. Error prevention 
scheme with four particles. Physical Review A, 54(3):R1745, 1996. 

Lorenza Viola, Emanuel Knill, and Seth Lloyd. Dynamical decoupling 
of open quantum systems. Physical Review Letters, 82(12) :2417, 1999. 

Lorenza Viola, Emanuel Knill, and Seth Lloyd. Dynamical decoupling 
of open quantum systems. Phys. Rev. Lett., 82:2417-2421, 1999. 

Lorenza Viola and Seth Lloyd. Dynamical suppression of decoherence 
in two-state quantum systems. Physical Review A, 58(4):2733, 1998. 

Rodney Van Meter, Thaddeus D Ladd, Austin G Fowler, and Yoshi- 
hisa Yamamoto. Distributed quantum computation architecture using 
semiconductor nanophotonics. International Journal of Quantum In¬ 
formation, 8(01n02):295-323, 2010. 


154 



[WFHll] 

[WHP03] 

[Wit 89] 
[WSB+04] 

[WY06] 

[Yoslla] 

[Yosllb] 

[YPA+03] 

[YR03] 

[ZLCOO] 

[ZR97] 

[Zur84] 


David S Wang, Austin G Fowler, and Lloyd CL Hollenberg. Surface 
code quantum computing with error rates over 1%. Physical Review A, 
83(2):020302, 2011. 

Chenyang Wang, Jim Harrington, and John Preskill. Confinement- 
higgs transition in a disordered gauge theory and the accuracy thresh¬ 
old for quantum memory. Annals of Physics, 303(l):31-58, 2003. 

Edward Witten. Quantum field theory and the jones polynomial. Com¬ 
munications in Mathematical Physics, 121 (3):351-399, 1989. 

Andreas Wallraff, David I Schuster, Alexandre Blais, L Frunzio, R-S 
Huang, J Majer, S Kumar, Steven M Girvin, and Robert J Schoelkopf. 
Strong coupling of a single photon to a superconducting qubit using 
circuit quantum electrodynamics. Nature, 431 (7005): 162-167, 2004. 

Pawel Wocjan and Jon Yard. The jones polynomial: quantum algo¬ 
rithms and applications in quantum complexity theory. arXiv preprint 
quant-ph/0603069, 2006. 

Beni Yoshida. Classification of quantum phases and topology of logical 
operators in an exactly solved model of quantum codes. Annals of 
Physics, 326(1): 15-95, 2011. 

Beni Yoshida. Feasibility of self-correcting quantum memory and ther¬ 
mal stability of topological order. Annals of Physics, 326(10) :2566- 
2633, 2011. 

Tsuyoshi Yamamoto, Yu A Pashkin, Oleg Astafiev, Yasunobu Naka¬ 
mura, and Jaw-Shen Tsai. Demonstration of conditional gate operation 
using superconducting charge qubits. Nature, 425(6961):941-944, 2003. 

Nadav Yoran and Benni Reznik. Deterministic linear optics quan¬ 
tum computation with single photon qubits. Physical review letters, 
91(3):037903, 2003. 

Xinlan Zhou, Debbie W Leung, and Isaac L Chuang. Methodology for 
quantum logic gate construction. Phy. Rev. A, 62(5):052316, 2000. 

Paolo Zanardi and Mario Rasetti. Noiseless quantum codes. Physical 
Review Letters, 79(17):3306, 1997. 

Wojciech H Zurek. Reversibility and stability of information processing 
systems. Physical Review Letters, 53(4) :391, 1984. 


155 



